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Introduction 



The title quantum optics covers a large range of possible courses, and so 
this introduction intends to explain what this course does and does not aim 
to provide. Regarding the negatives, there are several things this course 
deliberately avoids: 

• It is not a course on quantum information theory. Some basic notions 
of coherent states and entanglement will be assumed, but will not be 
the focus. 

• It is not a course on relativistic gauge field theories; the majority 
of solid state physics does not require covariant descriptions, and so 
it is generally not worth paying the price in complexity of using a 
manifestly covariant formulation. 

• As far as possible, it is not a course on semiclassical electromagnetism. 
While at times radiation will be treated classically, this will generally 
be for comparison to a full quantum treatment, or where such an 
approximation is valid (for at least part of the radiation) . 

Regarding the positive aims of this course, they are: to discuss how to 
model the quantum behaviour of coupled light and matter; to introduce 
some simple models that can be used to describe such systems; to dis- 
cuss methods for open quantum systems that arise naturally in the context 
of coupled light and matter; and to discuss some of the more interesting 
phenomena which may arise for matter coupled to light. Semiclassical be- 
haviour will be discussed in some sections, both because an understanding 
of semiclassical behaviour (i.e. classical radiation coupled to quantum me- 
chanical matter) is useful to motivate what phenomena might be expected; 
and also as comparison to the semiclassical case is important to see what 
new physics arises from quantised radiation. 

The kind of quantum optical systems discussed will generally consist 
of one or many few-level atoms coupled to one quantised radiation fields. 
Realisations of such systems need not involve excitations of real atoms, but 
can instead be artificial atoms, i.e. well defined quantum systems with dis- 
crete level spectra which couple to the electromagnetic field. Such concepts 
therefore apply to a wide variety of systems, and a variety of character- 
istic energies of electromagnetic radiation. Systems currently studied ex- 
perimentally include: real atomic transitions coupled to optical cavitiesp]; 



v 



vi 



INTRODUCTION 



Josephson junctions in microwave cavities (waveguides terminated by re- 
flecting boundaries) [21 |3J; Rydberg atoms (atoms with very high princi- 
ple quantum numbers, hence small differences of energy levels) in GHz 
cavities 0]; and solid state excitations, li.e. excitons or trions localised in 
quantum dots, coupled to a variety of optical frequency cavities, includ- 
ing simple dielectric contrast cavities, photonic band gap materials, and 
whispering gallery modes in disks [5]. 

These different systems provide different opportunities for control and 
measurement; in some cases one can probe the atomic state, in some cases 
the radiation state. To describe experimental behaviour, one is in gen- 
eral interested in calculating a response function, relating the expected 
outcome to the applied input. However, to understand the predicted be- 
haviour, it is often clearer to consider the evolution of quantum mechanical 
state; thus, both response functions and wavefunctions will be discussed. 
As such, the lectures will switch between Heisenberg and Schrodinger pic- 
tures frequently according to which is most appropriate. When considering 
open quantum systems, a variety of different approaches; density matrix 
equations, Heisenberg-Langevin equations and their semi classical approxi- 
mations, again corresponding to both Schrodinger and Heisenberg pictures. 

The main part of this course will start with the simplest case of a single 
two-level atom, and discuss this in the context of one or many quantised 
radiation modes. The techniques developed in this will then be applied 
to the problem of many two- level atoms, leading to collective effects. The 
techniques of open quantum systems will also be applied to describing las- 
ing, focussing on the "more quantum" examples of micromasers and single 
atom lasers. The end of the course will consider atoms beyond the two-level 
approximation, illustrating what new physics may arise. Separate to this 
main discussion, the first two lectures stand alone in discusing where the 
simple models of coupled light and matter used in the rest of the course 
come from, in terms of the quantised theory of electromagnetism. 



Lecture 1 



Quantisation of 
electromagnetism in the 
Coulomb gauge 

Our aim is to write a theory of quantised radiation interacting with quan- 
tised matter fields. Such a theory, e.g. the Jaynes-Cummings model (see 
next lecture) has an intuitive form: 



Hj.c. = Uka\a k + Vi°t + 9i,k°t a k + H - c -] • 



The operator al creates "a photon" in the mode with wavevector k, and so 
this Hamiltonian describes a process where a two-level system can change 
its state with the associated emission or absorption of a photon. The term 
0Jki\a k then gives the total energy associated with occupation of the mode 
with energy u^. While the rest of the course is dedicated to studying such 
models of coupled light-matter system, this (and in part the next) lecture 
will show the relation between such models and the classical electromag- 
netism of Maxwell's equations. 

To reach this destination, we will follow the route of canonical quan- 
tisation; our first aim is therefore to write a Lagrangian in terms of only 
relevant variables. Relevant variables are those where both the variable and 
its time derivative appear in the Lagrangian; if the time derivative does not 
appear, then we cannot define the canonically conjugate momentum, and 
so cannot enforce canonical commutation relations. The simplest way of 
writing the Lagrangian for electromagnetism contains irrelevant variables 
- i.e. the electric scalar potential (j) and gauge of the vector potential A; 
that irrelevant variables exist is due to the gauge invariance of the theory. 
Since we are not worried about preserving manifest Lorentz covariance, 
we are free to solve this problem in the simplest way — eliminating the 
unnecessary variables. 
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1.1 Revision: Lagrangian for electromagnet ism 



To describe matter interacting with radiation, we wish to write a La- 
grangian whose equations of motion will reproduce Maxwell's and Lorentz's 
equations: 



V x B = /x J + /U £oE 
V-B = 

m a r a = q a [E(r a ) + r a x B(r c 



V • E = p/e 
V x E = -B 



(1.2) 
(1.3) 
(1.4) 



Equations (1.3) determine the structure of the fields, not their dynamics, 
and are immediately satisfied by defining B = Vx A and E = — V^— A. Let 
us suggest the form of Lagrangian £ that leads to Eq. (1.2) and Eq. (1.4): 

C = Y1 \ m ^ 2 + § [dV[E 2 - c 2 B 2 ] + ^ Qa [r a • A(r a ) - 0(r a )] . 

a a 

(1.5) 

Here, the fields E and B should be regarded as functionals of <f> an d A. 
Note also that in order to be able to extract the Lorentz force acting on 
individual charges, the currents and charge densities have been written as: 



p(r) = ^q a 5(r - r a ), J(r) = } q a r a S(r - r a ). 



(1.6) 



The identification of the Lorentz equation is simple: 
d dC d 



d 



[m a r a + q a A(r a )] = m a r a + q a (r a ■ V)A(r a ) + q a — A(r Q ) 



dt dr a dt 
dC 

= t\ — = <?a V [r a ■ A(r a ) - 0(r Q )] 

= q a [(r Q • V) A(r a ) + r a x (V x A(r a )) - V0(r a )] , 



thus one recovers the Lorentz equation, 



r a x [V x A(r a )] - V<£(r a ) 



at 



A(r a ) 



(1.7) 



Similarly, the equation that results from can be easily extracted; since 
dC/d(j) = 0, this becomes 



dC 



e V ■ E - ^2q a 5(r - r c 



0. 



(1.8) 



Finally, the equations for A are more complicated, requiring the identity 



OA 



(V x A) 2 = 2V x (V x A), 



(1.9) 



which then gives: 
d dC d 



dC 



dtdA dt OA 



-e c 2 V x (V x A) + ^2q a r a 5(r - r a ), 

a 

(1.10) 
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which recovers the required Maxwell equation 

- eo^E = --^V x B + ^2q a r a 5(r - r a ). 



Mo 



(1.11) 



Thus, the Lagrangian in Eq. (1.5), along with the definitions of E and B 
in terms of A and eft produce the required equations. 



1.2 Eliminating redundant variables 

As mentioned in the introduction, we must remove any variable whose 
time derivative does not appear in the Lagrangian, as one cannot write 
the required canonical commutation relations for such a variable. It is 



clear from Eq. (1.5) that the electric scalar potential <j> is such a variable. 



Since <p does not appear, it is also possible to eliminate <$> directly from 



the equation dC/dcf); using Eq. (1.8) and the definition of E, this equation 
gives: 

-e V- A - e o V 2 - p(r) = 0. (1.12) 
Rewriting this in Fourier space, one has: 

0(k) = ^f^ + ik-A(k)V (1.13) 



= 



We can now try to insert this definition into the Lagrangian, to eliminate 
(f). To do this, we wish to write E 2 and B 2 in terms of cfi and A; it is therefore 
useful to start by writing 

- Ej(k) = ik^ik) + Aj(k) = i^p(k) + (s jk - i fe (k). (1.14) 

This means that the electric field depends on the charge density, and on 
the transverse part of the vector potential, which will be written: 

4 k = - ^Sr) ^( k )- (i-iB) 

The transverse 1 part of the vector potential is by definition orthogonal to 
the wavevector k, and so the electric field is the sum of two orthogonal 
vectors, and so: 

|E(k)| 2 = |A^(k)| 2 + ^lp(k)| 2 . (1.18) 



The combination: 

^(k)=fe- fc ^) (j.i.i, 



is the reciprocal space representation of the transverse delta function; with appropriate 
regularisation[6] Complement Ai], it can be written in real space as: 

Sf k (r) = 2 -6(r)S 3k + ( 3 ^ - S jk ) . (1.17) 
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Similarly, the squared magnetic field in reciprocal space is given by: 
|B(k)| 2 = (k x A(k)) • (k x A*(k)) 



= k 2 5< 



kjkk 



A j (k)A* k (k) = k 2 \A ± (k)\ 2 . 



(1.19) 



2 



Thus, the field part of the Lagrangian becomes: 
dV [E 2 -c 2 B 2 ] = ^-J d 3 k^\p(k)\ 2 

+ £ ^d 3 fc(|A ± (k)| 2 -c 2 fc 2 |A ± (k)| 2 ) . (1.20) 

The notation f has been introduced to mean integration over reciprocal 
half-space; since A(r) is real, A*(k) = A(— k), thus the two half spaces 
are equivalent. This rewriting is important to avoid introducing redundant 
fields in the Lagrangian; the field is either specified by one real variable at 
all points in real space, or by two real variables at all points in reciprocal 
half-space. Similar substitution into the coupling between fields and matter, 
written in momentum space gives: 



-cm— matter 



= 2^kj d 3 k[J(k) ■ A*(k) - p(k)<f (k)] 



i 



2m d 3 k 



J(k)-A*(k)-p(k) 



\e k 2 



ik 

¥ 



A*(k) 



:i.2ii 



This can be simplified by adding a total time derivative, C — > C + dF/dt; 
such transformations do not affect the equations of motion, since they add 
only boundary terms to the action. If: 

-tk- A*(k)' 



2'Stj d 3 kp(k) 



k 2 



then 
C 



cm— matter 



dF 
~dt 



2^kjd 3 k ^ 



J(k)-p(k) 



ik 



A*(k) 



(1.22) 

ML 21 

e k 2 



(1.23) 

Then, using conservation of current, p(k) + ik ■ J(k) = 0, one finally has: 



£ cm— matter 



dF 
~dt 



2)tj d : ''k 



J^k) • A ±; "(k) 



e k 2 



(1.24) 



Note that this set of manipulations, adding dF/dt has eliminated the longi- 
tudinal part of the vector potential from the Lagrangian. The form chosen 
for F is such that this procedure is equivalent to a gauge transformation; 
the chosen gauge is the Coulomb gauge. Putting everything together, one 
has: 



-^coulomb 

+ £o 



,j d 3 k (jA^k)! 2 - c 2 A; 2 | A ± (k)\ 2 + 2K [j x (k) ■ A x *(k)] ) . (1.25) 
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Thus, the final form has divided the interaction into a part mediated 
by transverse fields, described by A ± , and a static (and non-retarded) 
Coulomb interaction. Importantly, there are no irrelevant variables left 



in Eq. (1.25) The Coulomb term can also be rewritten: 



Vcoul = -/^l^k)! 2 = £ - ^ (1.26) 

Note that since the Coulomb interaction is non-retarded, both the Coulomb 
and transverse parts of interaction must be included to have retarded in- 
teractions between separated charges. 



1.3 Canonical quantisation; photon modes 



We have in Eq. ( 1.25| ) a Lagrangian which can now be treated via canon- 
ical quantisation. Since only the transverse part of the field A appears in 



Eq. (1.25), we can drop the superscript label in A from here on. To pro- 



ceed, we should first identify the canonical momenta and the Hamiltonian 
and then impose canonical commutation relations. Thus, 



n(k) 



dr a 

9£ coulomb 

<9A(k)* 



mr a + q a A(r a ) 
e A(k). 



(1.27) 
(1.28) 



Then, constructing the Hamiltonian by H 

1 r 



H 



y 

^ 2m a 

a 

+ £ 7 d 3 k 



Pc 



^2iPiRi — £, one finds: 
q a A(r a )} 2 + V COVL i 



|n(k) 



+ c 2 k 2 \A(k)\- 



[1.29) 



In order to quantise, it remains only to introduce commutation relations 
for the canonically conjugate operators. Noting that A(r) has only two 
independent components, because it is transverse, it is easiest to write 
its commutation relations in reciprocal space, introducing directions ek, n 
orthogonal to k with n = 1, 2 ; then: 



[ri,a,Pj,p] = ihSapbij 

A e . „(k),n ek , .(k')l =ih5(k-k')5 r , 



;i.3o) 
;i.3i) 



This concludes the quantisation of matter with electromagnetic inter- 
actions in the Coulomb gauge. It is however instructive to rewrite the 
transverse part of the fields in terms of their normal modes. The second 



line of Eq. (1.29) has a clear similarity to the harmonic oscillator, with 



a separate oscillator for each polarisation and momentum. Rewriting in 
normal modes thus means introducing the ladder operators: 



2hck 



c^ ek , ra (k) + -n eki „(k) 



(1.32) 
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or, inverted one has: 



A(r) 



EE 

k n=l,2 



(1.33) 



Inserting this into Eq. (1.29) gives the final form: 



H = ^2 ^ a ~ qa ^ Ya ^ 2 + ^2 ^2 ^kaka k + v con \. 



(1.34) 



k n=l,2 



1.4 Dipole, two- level, and rotating wave 
appr oximat ions 



Equation (1.34) applies to any set of point charges interacting with the 
electromagnetic field. In many common cases, one is interested in dipoles, 
with pairs of opposite charges closely spaced, and much larger distances 
between the dipoles. In this case, there are a number of approximations 
one can make to simplify calculations. This section will briefly illustrate 
these approximations: neglect of the A 2 terms, the dipole approximation, 
projection to two-level systems, and the rotating wave approximation. The 
study of the approximate models that result will be the subject of the rest 
of this course. 

Neglect of the A 2 terms in expanding [p a — q a A(r a )] 2 can be justified 
in the limit of low density of dipoles; considering only a single radiation 



mode, the contribution of the A 2 term can be rewritten using Eq. (1.33) 
as: 

SH A » = ?--±^(a+tf) . (1.35) 
A V 4me w V J v ' 

Thus, this term adds a self energy to the photon field, which scales like the 
density of dipoles. The relative importance of this term can be estimated 
by comparing it to the other term in the Hamiltonian which is quadratic 
in the photon operators, fkJka\.a k . Their ratio is given by: 

N q 2 N 3 (IZydV 

thus, assuming particles are more dilute than their Bohr radius, neglect 
of 5H A 2 is valid for frequencies of the order of the Rydberg for the given 
bound system of charges. 

Turning next to the dipole approximation, consider a system with just 
two charges: charge +q mass m\ at R + r/2 and — q mass m2 at R — r/2. If 
r <C A where A is a characteristic wavelength of light, then one may assume 
that A(R + r/2) ~ A(R — r/2) ~ A(R), and so the remaining coupling 
between radiation and matter is of the form 

<W A . P -'/(—- — ) -A(R). 
1 mi m-2 



1.5. FURTHER READING 
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Then, in the case that one can write Hq, a Hamiltonian for the dipole 
without its coupling to radiation one can use p/m = x = i[Ho,x]/h, thus 
giving: 

H = H -i^[H ,r}-A(K) + ^ fo^on + W (1.37) 

k,n=l,2 

To further simplify, one can now reduce the number of states of the 
dipoles that are considered; currently, there will be a spectrum of eigen- 
states of Ho, and transitions are induced between these states according to 
(ipf\r\ipi). Finally, restricting to only the two lowest atomic levels and to 
a single radiation mode, one has a model of two-level systems (describing 
matter) coupled to bosonic modes (describing radiation). This model is 
known as the Jaynes-Cummings model: 

u 1 ( 6 g{a + a)) \ . k + n iq\ 

H = - , + . )+nu)na ] a. (1.38) 

Here we have introduced the energy splitting e between the lowest two 
atomic levels. In terms of the upper and lower states, \a) and \b), the 
atom-photon coupling strength g can be written: 

^=q(b\[H ,r)\a) - 1 = e Gk ' n ' dba , (1.39) 
2 y/2sQf}WkV ^2eohwkV 

where dj, a = g (6|r|a) is the dipole matrix element. The matrix notation 
represents the two levels of the dipole. We assumed Hq commutes with the 
parity operator x — > —x, and so the coupling to radiation appears only in 
the off-diagonal terms in the two-level basis. The final approximation to 
be discussed here, the rotating wave approximation, is appropriate when 
e ~ wo- Considering g as a perturbation, one can identify two terms: 

9 ( a \ . g ( at \ 

A - = ^ „t n . A cross = ^ n n , (1.40) 



2\at J ' cross ~ 2 V a , 

where A co , the co-rotating terms, "conserve energy"; and A cross do not. 
More formally, the effects of A croS s give second order perturbation terms 
like g 2 /(u>o + e), while A co give the much larger g 2 /(cjq — e)- 



1.5 Further reading 

The discussion of quantisation in the Coulomb gauge in this chapter draws 
heavily on the book by Cohen- Tannoudji et al. [6] . 



Questions 



Question 1.1: Transverse delta function 



Prove that Eq. (1.17) is the Fourier transform of Eq. (1.16). It is helpful 
to consider a modified version of Eq. (1.17), with a factor exp(— mk 2 ) 
take m — > only at the end of the calculation. 



and 
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Question 1.2: Thomas-Reiche-Kuhn sum rule 

Prove the equality: 

[r a ,[H ,r p }} = ^, (1.41) 

stating the most general form of Hq for which this is true. 

From the expectation of this commutator, show that the dipole matrix 
elements between any state a and all other states b obeys the relation: 

£ I4 b | 2 (£ 6 - £ a ) = £ (1-42) 

b a 



Lecture 2 



Gauge transformations, and 
quantum electrodynamics in 
other gauges 



In the previous section, we choose to work in the Coulomb gauge, adding 
a total time-derivative to the Lagrangian which had the effect of removing 
all dependence of the Lagrangian on An. This choice, which is equivalent 
to choosing to impose A» = hugely simplified the subsequent algebra, but 
is not strictly necessary. This chapter describes the consequences of mak- 
ing other choices for Air, such choices turn out to correspond to making a 
unitary transformation in the quantum problem. In the special case of one 
or a few localised systems of charge — where charges within a system are 
separated by far less than the wavelength of light — a change of gauge to 
the electric dipole gauge can simplify calculations, and provides further un- 
derstanding of the relation between the instantaneous Coulomb interaction 
and the photon-mediated terms in Eq. (1.34). 



2.1 Freedom of choice of gauge and classical 
equations 

The classical Lagrangian in the previous section depended on the longitu- 
dinal part of the vector potential only in the coupling between matter and 



electromagnetic fields, Eq. (1.21) 



£em-matter = 2%tj d 3 k [J(k) • A*(k) - p(k)0*(k)] 



2${jd 3 k 



J(k).A-(k)- /)( k)(£^-^.A*(k) 



(pi) 
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LECTURE 2. QED IN OTHER GAUGES 



Using the continuity equation p + ik ■ J = to eliminate Jy = ip/k and 
breaking A into transverse and longitudinal parts yields 



23U d 3 k 



j ± (k).Ai(k) + ip(k)Aj(k)-J^ 



+ ^P(k)Af (k) 



It is therefore clear that taking any functional form for 
to a change to the Lagrangian which may be written as: 



5L 



gauge 



d A k 



k 



99* (Aj 



(2.1) 
<?(Aj_) leads 



(2.2) 



The choice of A» is the choice of gauge; by choosing A» as a function of 
Aj_, one can ensure, for example, that n • A(r) = for some fixed n. Since 
such a condition corresponds only to an additional total derivative, the 
action changes only by boundary terms, and so the equations of motion are 
not affected. The definition of canonical momentum however does change. 
Rewriting 

d_ 

at 



5L 



gauge 



w 



drk 



k 



pg* (A ± ) 



+ 2%ljd 3 k 



i dg* (A±) A 

k P QA A ± (kj 



the canonical momentum thus becomes: 

dC 



n(k) 



3A ± (k)* 



n 



Coulomb 



(k) + 



ip(k) dg* (A. 



k 



dA 



(2.3) 
(2.4) 



Quantum formalisms resulting from different gauges 

Once we have found the canonical momentum in a given gauge, we can 
quantise in this gauge by promoting the dynamical variables and their 
canonical conjugates to operators, and imposing canonical commutation re- 
lations. This means [Aj_ n (k), n n /(k')] = ih5(k — k')8 nn ', with n(k) being 
the new canonical momentum in the new gauge: The canonical momentum 
corresponds to a different combination of physical fields in different gauges. 
If we write 



F{A ± ,p) = 2^ftj d 3 k ] 

then we may compare the two formalisms: 
A_L >new (k) =A_L j0ld (k) 

n new (k) =n oM (k) + 



pg* (A. 



dF 



dA 



(2.5) 

(2.6) 
(2.7) 



..new (k)' 

If we consider a state \ir) which is an eigenstate of the old momentum, 
n o id(k)|7r) = \n\ir), then in the new formalism, this state obeys: 

dF 



n new (k) 



<9A_L incw (k) 



Itt) 



-ih- 



d 

dA±^ new 

iF(A ± ,p)/h 



(k) 

ill 



dF 
5A_L incw (k) 



Itt) 



d 



5A 



,new 



(k) 



e -iF(A ± ,p)/h>y 



(2.8) 

(2.9) 
(2.10) 
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Hence, the eigenstates of the momenta of the old and new formalism are re- 
lated by the unitary transform exp[—iF(A.±,p)/K] — one can immediately 
see that this relation of eigenstates of operators works for all canonical mo- 
menta, as well as working for the position operators which commute with 
the unitary transforms. It is worth stressing that the statement proved here 
is that all gauge transforms correspond to unitary transforms of the quan- 
tum formalism, and not the other way around. It is also worth highlighting 
(as was implicit in the calculation), that this occurred because the gauge 
transforms corresponded classically to canonical transformations of the La- 
grangian (addition of a total time derivative), and that it is more generally 
true that canonical transformations of classical dynamical systems corre- 
spond to unitary transformations of the associated quantum problem. 



Formal equivalence of gauges 

Formally, since a gauge transformation corresponds to a unitary transfor- 
mation, calculations of physical quantities in different gauges should exactly 
match, i.e. 



(i) 



(D 



(i) 



(1) 



U j UO (1) U^U 



(1) 



(2) 



G (2) 



(2) 



(2.11) 

Despite this formal equivalence, there have been long and convoluted 
arguments about gauge invariance in approximate methods. One particu- 
larly dangerous feature of such transformations is that the "natural" iden- 
tification of a bare and perturbation Hamiltonian can be different between 
different gauges (this is true of the discussion below), and so eigenstates 
of the bare Hamiltonian are not a gauge invariant quantity. Two general 
solutions to such problems are 

• To use a basis of states that correspond to eigenstates of a physical 
(i.e. gauge invariant) operator, such as the mechanical momentum, 
physical fields, and the total energy. 

• To use only the formal equivalence, and demand a transformation of 
states as well as of operators when switching between gauges. 



2.2 Transformation to the electric dipole gauge 

The preceeding working allows us to directly transform between two quan- 
tum pictures, rather than needing to changes gauges classically and re- 
quantize. As an illustration of this in a case where it is particularly clear, 
this section describes the transformation to the electric dipole gauge for a 
problem of two systems of localised charges, where each system is overall 
neutral. We consider a hierarchy of scales, such that within a given system 
all charges are within a distance I and I< A, |R^4 — Rb|, where Ra_b are 



the centres of mass of the two systems. This is illustrated in Fig. 2.1 

Since ! < A, we can approximate the electromagnetic fields as being 
those at the centre of mass of each system. In the Coulomb gauge, this 
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B 



Figure 2.1: Two localised systems, consisting of charges which 
are closely spaced in each system. 



then can be written following the notation in Eq. (1.34) as: 

+ ^ Wk /4a k + Kll + + Kil, (2-12) 

k,n 

in which the Coulomb interaction has been divided into parts limited to each 
subsystem as well as a cross term. This cross term is instantaneous (as all 
the Coulomb terms are), while the real physical interaction is retarded. To 
recover this retardation it is however necessary to include the transverse 
parts of the fields. From this point onwards we will set K = 1 

Transformation of Hamiltonian 

The unitary transformation that will simplify the results is of the form: 



U = exp 



[id A • A_l(Ra) + id B ■ Aj_(Rb)] , = ^ q a (r a - R A ), 

(2.13) 

in which we have introduced the dipole moments d Aj B of the two sys- 
tems. This transformation clearly commutes with the position operator, 
and ought to commute with the transverse vector potential (this is shown 
to be true explicitly below). However, it does not commute with p a nor 
a iL' a k- O ne ms t ea d finds: 

U ] p a U = p a - i-^— ( i ^ <I»{r a - R A ) • A_l(R a ) ) 

dr « \ a J 
= p Q + q a A ± (R A ) - Qa^j ■ A±(Ra) + ° (x) • 

(2.14) 



In this, the term ^2 a q a = due to neutrality of each subsystem, and the 
higher order terms in l/X are ignored due to the long wavelength approxi- 
mation. Hence all that remains is Wp a U ~ p a + q a A±(R A ). 
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For the photon annihilation operators no approximations are needed. 
The unitary transformation can be rewritten as 



U = exp < 



d A e 



d B e 



-ik-Ri 



4 + H.c.L (2.15) 



which has the form of a shift operator for the annihilation operators: 



U^a\JJ = a k + i 



d A e 



+ d B e- ik ' Rs 



(2.16) 



From this we can reconfirm that the unitary matrix commutes with the 
transverse vector potential, 



k,n 



y/2e uj k V L 



a ik r al + e ik r a k 



(2.17) 



hence, 



U*A ± (r)U = A ± (r) 



+ E^V e n -d Ae jk -( r - R -)+e n .d se lk -( r - R -)-H.c. . (2.18) 



k,n 



2e ^kV 



It is clear that due to the antisymmetry of the sum under k <-> — k the 
sum vanishes, so Aj_(r) commutes with the transformation. Hence, the 
transformed form of the Hamiltonian is given by: 



AB 
coul 



+ 



+ E 

k,n 
1 

2^V 



UJ k 



v / 2^ k 17 



(e n - d A ? + (e n - d B ) 2 + (e n - d A ) (e n - d B ) 2 e ik '( R ^- R ^ 



(2.19) 



where the last line has relabelled the sum via k <-> — k. 

Transformation of physical fields 

To simplify this expression, we need two connected things: 

• To identify the combination of photon operators that appears in the 
second line as the light-matter coupling. 



To simplify the cross terms between the A and B systems. 
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It is tempting to use the definition of the transverse electric field in the 
old gauge, 



2e V 



ik-r t ik-r 



(2.20) 



to identify the combination of operators, however this relation between 
the physical operators and the annihilation operators is true only in the 
Coulomb gauge, not the electric dipole gauge. Transforming this electric 
field using Eq. (2.16) gives: 



tfE ± U = E ± -Y 



k,n 



e 

2e V L 



dA£ ,ik.(r-R A ) + dse ik.(r-R B ) 



+ H.C. 



(2.21) 



This sum is symmetric under k <-> — k and so does not vanish. Instead 
it can be identified in terms of the polarisation. Since the dipole matrix 
element of system A is d^, we can write PA( r ) = d^<5(r — Ra); this places 
the polarisation at the centre of charges, and is measured w.r.t a reference 
state in which all the charges are at the centre of mass, i.e. Vr Q = R A . If 
we then use the transverse delta function to find the transverse part of this 
expression, we find that: 



_(r) = ^J>„ [e n • d A e* k -( r - R ^ 



(2.22) 



and hence: 



k,ra 



2e V 



k,n 



E±(r) + EiW = E±M (2 . 23 ) 



= 



This solves our first question; the second is also considerably simplified 
when we note that the term: 



E 47^- d ^(^- d ^ eJk ' (RA " Rfl) 

k,n 

d 3 r 




k,n 

f -PAl(r)P B ,l(r). (2.24) 
J eo 



Thus, we can write: 

H = y vL + y 



.->,.,, /__> ^coul ^coul Kiioul 



+ Yl ^ a W + D ±(Rvl) • d A + D ± (R B ) • d B 

k,n 

+ Sf ipolar + Sf ipolar + J ^P A ±(r)P B ,±(r). (2.25) 



2.3. ELECTRIC DIPOLE GAUGE FOR SEMICLASSICAL 
PROBLEMS 
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Instantaneous nature of interaction 

To complete the simplification, it is necessary to write the cross-term in 
the Coulomb interaction similarly to Eq. (2.24). The Coulomb interaction 
arose by writing the longitudinal part of the electric field in terms of the 
charge density. Inverting this procedure easily gives: 

V C Z = jjd 3 r [(E M (r) + E M (r)) 2 - E M (r) 2 - E M (r) 2 ] (2.26) 
= e J d 3 rE M (r) • E M (r) = j- J d 3 rP lA (r) ■ P M (r). (2.27) 

The last step made use of the reference distribution of charges; since this 
distribution has no net charges (each system is neutral, and all charges 
coincide), then V • D = 0, i.e. Dm = 0, so En = Pii/eo- This allows us to 
combine: 

+ [ -PA,l(r)PRl(r) = / — P,(r)P B (r) = (2.28) 
J eo J e 

where the last line follows from the definition of the polarisation, which is 
zero outside each region of charge. 
Hence, the final Hamiltonian is: 

'2 n 2 

fy _ V" Pa , \" y P , t/AA , V BB , V A , V B 

9^T + 9^T + coul + coul + ^dipolar + ^dipolar 
a ^ nl a a z,n f3 

+ ^Wk4a k + D ± (R A ) -d A + D ± (R B ) d B . (2.29) 



k,n 



This transformation has therefore removed the instantaneous interac- 
tion between the two systems; all interaction is now described by the quan- 
tised D field. Since Dm = 0, the transverse and total fields are equivalent, 
hence the transverse D field is physical and retarded. 



2.3 Electric dipole gauge for semiclassical 
problems 

The previous section showed the unitary transformation to the electric 
dipole gauge for quantised radiation coupled to quantised matter. For 
the simpler case of a single localised system, the transformation would be 
U = expfzd • Aj_(0)], yielding: 

„2 

i* = 12 T 9 - + Fcoul + Sdi P° lar + w k4 a k + D ± (0) • d. (2.30) 

a lrUa k,n 

For comparison, this section describes the even simpler transformation 
that applies when the radiation is assumed to be classical; i.e. it has no 
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free dynamics, but rather the operator Aj_ is replaced by a time depen- 
dent field corresponding to some incident radiation. In this case the initial 
Hamiltonian is: 

H s .c. = J2 r [Pa ~ g«A ± (0, t)f + Koui. (2.31) 

a 

The gauge transformation is formally the same as before, but as now Aj_ 
is not an operator but a time-dependent field we must use the result: 

WJW) = Uidtty) + (id t U)\iP) = HU\i>) (2.32) 

idt\i/>) + (uUd t U)\iP) = U^HU\iP) (2.33) 

and so: 

H = U ] HU - UUd t U. (2.34) 



The transformation of p a is exact as in Eq. (2.14). Since Aj_ is a c- 
number, it commutes with the unitary transformation, and so the final 
result is: 

Hs.c. = V + Koui - idt [id • A_l(0)] 

a 

III (y 



Relating ^-matrix elements in different gauges 

In the semiclassical well as the formal equivalence between gauges, 

there are "less formal" relations that can be shown to hold. One example 
arises if we consider a field Aj_(t) = AoA(t) cos(a>t), where \{t) is a function 
that goes to at ±oo, but \(t) = 1 for t\ < t < t2, with smooth interpola- 
tion between: Since this means U(t) = exp[«A(t)d • Aocos(wi)], it is clear 

Figure 2.2: Form of X(t) for S'-matrix 




that at long times U{t) = 1, and hence eigenstates in the two gauges are 
coincident. 

This allows us to consider the S'-matrix: 

Sab = (b\ e iH ^U{hM)e- iHotl \a) . (2.36) 

Since U = 1 at long times the bare Hamiltonian Hq is also the same between 
the two gauges, so this S'-matrix should be gauge invariant. Let us consider 
its calculation perturbatively in the two gauges 



2.4. PITFALLS OF PERTURBATION 
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Calculation in A p gauge 

In the Coulomb gauge, the perturbation Hamiltonian (to first order in Ao) 
is given by: 

H int (t) = - • Aocos(arf). (2.37) 

z — ' m a 

a 

Using the result of time-dependent perturbation theory: 

U(t 2 ,h) = l-i f 2 dre-^ o(t2 - r) i7 int (r)e- l//o(r -' l) + . . . (2.38) 
Jti 

we have: 

rt 2 



S ab = Sab~i / 2 dre l ^ T (b\ H int (r) \a) (2.39) 
Jti 

= 5 ab + i f 2 dre iuJbaT cos(^) A • (b\ 



PaQa i v 

a 



x , o ^(^ba + w) + 5(w 6a - w) /Tj-V^rrr II \ 

= <5 a6 - 7r [<5(cj ba + w) + <5(w fea - a;)] uj ba A ■ (b\ d \a) . (2.40) 

In the last line we have used the fact that the bare Hamiltonian depends 
only on p a via p^/2m a , allowing one to rewrite the matrix element of 
momentum in terms of the dipole matrix element. 

Calculation in E • r gauge 

In the electric dipole gauge, the perturbation is instead 

fl int (t) = -d-E(t). (2.41) 



where E(i) = — A(t) ~ Aquj sin(u>i), and so starting from Eq. (2.39) we 
have: 

rti 

Sab = Sab + i / dTe iWbaT sin(a;t)a;Ao • (b\ d |a) 
Jti 

= 5ab-n [S(uj ba -u)- 5(uj ba + oj)] wA • (b\ d |a) , (2.42) 

where we have made the same resonant approximation as above. Although 
functionally different, the two expressions do match, because uj ba = w is 
ensured by the delta function. For a further example of this kind of equiv- 



alence, see question 2.1 



2.4 Pitfalls of perturbation 

The idea of writing the free Hamiltonian Hq is in fact rather subtle to define 
correctly: Atoms are held together by electromagnetic forces, and as men- 
tioned above, separation into transverse, photon mediated, and Coulomb 
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terms leads to apparently non-retarded interactions. Further, Hq differs 
between different gauges. 

The question of most practical relevance is how the parameters in a 
semi-phenomenological model of a coupled light-matter system should be 
related to experimentally measured quantities. 

2.5 Further reading 

The majority of this chapter is also based on discussions in the book by 
Cohen- Tannoudji et al. [61. The subtlety mentioned briefly, about how 
one can define a "free Hamiltonian" for atoms when the binding within an 
atom is due to the electromagnetic field has been discussed in the context 
of which gauge should be used in calculating off-resonant transitions. See 
for example [7H2]- 

Questions 

Question 2.1: Two-photon transitions 

If we choose states \b), \a) such that (6| d \a) = 0, then the leading order 
contribution to the transition comes from two photon processes. Assuming 
a b, the leading order contribution to the S'-matrix can be written: 



Sab = - dr dr'e^e^' (b\ H int (r) \c) (c\ H mt (r') \a) (2.43) 



2.1. (a) Making the resonant approximation, show that the expressions 
for S a b m the two gauges can be written as: 



A = 2 / H -z^ ^ ^l Ao . (6 | d )c) Ao . (c) d , o) (2.44) 

4 ^ LO ca - LO 

S^ = 2,^M £ ^^Ao • (6| d |c) Ao • (c\ d \a) (2.45) 

4 — UJ ca — U 



2.1.(b) Prove the equivalence of these expressions by showing that: 




c 



2 



■ (b\ di\c) ■ (c\ dj | a) = 



(2.46) 



c 



UJ ca - UJ 



if 2u = u ba . 

[Hint: use the two-photon resonance condition to show that 
D oc ^(wfe c - uj ca ) • (b\ di\c) • (c| dj \a) , 



(2.47) 



c 



then rewrite this sum in terms of the commutator [[Hq, di], dj].] 



Lecture 3 



One two-level system 
coupled to photons: the 
Jaynes-Cummings model 



The aim of this first section is to discuss the interaction between a single 
two- level system and various different configurations of photon fields. We 
will begin by discussing the interaction with a semiclassical light field, then 
consider a single quantised mode of light in an initially coherent state, and 
finally consider a continuum of quantised modes. In future lectures we will 
then consider combinations of these cases, to describe realistic cavity QED 
systems, in which the cavity picks out a particular mode, but coupling to a 
continuum of background modes also exists. Even in the case of a classical 
applied field, the behaviour of a two-level system is not entirely trivial, 
as the two-level system is non-linear, i.e. it does not have an harmonic 
spectrum. 

3.1 Semiclassical limit 

On applying a classical light field to a two-level system, we can adapt 
the Jaynes-Cummings model by neglecting the dynamics of the photon 
field, and replacing the photon creation and annihilation operators with 
the amplitude of the time-dependent classical field. This gives the effective 
Hamiltonian: 



We have made the rotating wave approximation, and hence only included 
the positive/negative frequency components of the classical field in the rais- 
ing/lowering operators for the two-level system. To solve this problem, and 
find the time-dependent state of the two-level atom, it is convenient to 
make a transformation to a rotating frame. In general, if 




(3.1) 




(3.2) 



19 



20 



LECTURE 3. JAYNES CUMMINGS MODEL 



then 



„ ( Hn + f(t) H 12 e~ 2 ^ , 

H ^ H ~i ff 2ie W H 22 -f(t) )' (3 - 3) 



thus, in this case we use / = —tot/2, giving: 

H= 1 ( e ~ UJ ^ A (3-4) 

Since this transformation affects only the phases of the wavefunction, we can 
then find the absorption probability by considering the time averaged prob- 
ability of being in the upper state. The eigenvalues of Eq. ( |3.2[ ) are ±0/2 
where O = y(e — uj) 2 + g 2 X 2 . Writing the eigenstates as (cos 9, sin 9) T , 
and (— sm9, cos 6) T corresponding respectively to the ± roots (assuming 
tan(#) > 0). Substituting these forms, one finds the condition: tan(2#) = 
g\/(e — uj). This form is sufficient to find the probability of excitation. If 
the two-level system begins in its ground state, then the time-dependent 
state can be written as: 

t T ) = sin e( cos M e -^/ 2 + cos B ( - sin e ) (3 5) 

ipi J \ smff J \ cos 9 J 



i — cos(rjf) 



Pex = |^T| 2 = |2sinecos0sin(m/2)|2- ,,,,,, 

(3.6) 

Thus, the probability of excitation oscillates at the Rabi frequency £1 = 
— uj) 2 + OA) 2 , and the amplitude of oscillation depends on detuning. 
On resonance (e = uj), one can engineer a state with = = 
or a state with |^| = 1 by applying a pulse with a duration (gA)i = t 1 "/ 2 
or 7r respectively. 

3.2 Single mode quantum model 

Let us now repeat this analysis in the case of a quantised mode of radiation, 
starting with the Jaynes-Cummings model in the rotating wave approxima- 
tion. 



H=-[ t + uJ a)a. (3.7) 



2 V 9a j -e 

Rabi oscillations in the Jaynes-Cummings model 

In the rotating wave approximation the total number of excitations iV ox = 
a z + a) a is a conserved quantity, thus if we start in a number state of the 
light field: 

\n,l) = ^ 7 L\0,l), (3.8) 
\Jn\ 



then the general state at later times can always be written: 

\m)=a(t)\n-l,^)+p(t)\n,i). (3.9) 



3.2. SINGLE MODE QUANTUM MODEL 
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Inserting this ansatz in the equation of motion gives 

1 



id t 



a{t) 



ii 



ul + 



[e-uj) 

gy/n 



9V n 
(e-uj) 



a(t) 
Pit) 



(3.10) 

After removing the common phase variation exp[— i(n— 1 /2)ut[, the general 
solution can be written in terms of the eigenvectors and values of the matrix 
in Eq. (3.10), and reusing the results for the semiclassical case, we have the 
state 



a{t) 



-i(n-l/2)u)t 



sin 



e -int/2 



COS I 

sin i 



+ cos I 



JQt/2 



- sin I 
cos I 



and the excitation probability is thus: 

1 — cos(f2t) 
cx = [ 2 

where we now have tan(2#) = g\fn/{e 



g 2 n 



g 2 n + (e — uj) 2 
uj) and VL = y^(e 



(3.11) 
(3.12) 



uj) 2 + g l n. 



Collapse and revival of Rabi oscillations 

Restricting to the resonant case e = uj, let us discuss what happens if the 
initial state does not have a well defined excitation number. Since the oscil- 
lation frequency depends on n, the various components of the initial state 
with different numbers of excitations will oscillate at different frequencies. 
We can therefore expect interference, washing out the signal. This is indeed 
seen, but in addition one sees revivals; the signal reappears at much later 
times, when the phase difference between the contributions of succesive 
number state components is an integer multiple of 2ir. 

Let us consider explicitly how the probability of being in the excited 
state evolves for an initially coherent state 



-l A l 2 / 2 exp(Aat)|0,i) = e-l A l 2 / 2 ^ 



A' 



(3.13) 



If resonant, the general result in Eq. (3.11 ) simplifies, as 9 = tt/4. Following 
the previous analysis, the state at any subsequent time is given by: 

lo/A - P HA| 2 +«^)/2 



E 



(Ae- 



cos l — - — I |n, 1) + sin — - — ) \n — 1, J) 



(3.14) 



Thus, the probability of being in the excited atomic state, traced over all 
possible photon states, is given by 



P.- 



n 

1 1 



2 = e -|A|- 



V- |A| 2n . 2 (gyfUt 



111 



111 



[gy/nt) 



(3.15) 
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We consider the case A>1. In this case, one may see that the coefficients 
|A| 2n /n! are sharply peaked around n ~ |A| 2 by writing: 

^— = -=exp[-/(n)], f(n) =nln(n) - n - nln(|A| 2 ), (3.16) 
n\ V27rn 

By differentiating / twice, one finds f'(n) = ln(n/|A| 2 ), /" = 1/n, and 
hence one may expand f(n) near its minimum at n = |A| 2 to give: 

/(n = |A| 2 + x)^-|A| 2 + ^ (3.17) 



Using this in Eq. (3.15) gives the approximate result: 

^2 



P ex ~ jft < exp 

2 2y ^iw itr 



m 

m 



2|A|- 



(3.18) 



where it is assumed that | A| 2 is large enough that the limits of the sum may 
be taken to ±oo. 

Three different timescales affect the behaviour of this sum. Concentrat- 
ing on the peak, near m = 0, there is a fast oscillation frequency gX, de- 
scribing Rabi oscillations. Then, considering whether the terms in the sum 
add in phase or out of phase, there are two time scales: collapse of oscilla- 
tions occurs when there is a phase difference of 2tt across the range of terms 
\m\ < a m contributing to this sum, i.e. when gT co n apse ( y/\ A| 2 + a m — |A|) = 
it. Taking a m = |A| from the Gaussian factor, this condition becomes 
(7T co ii a psc — 27T. On the other hand, if the phase difference between each suc- 
cessive term in the sum is 27r, then they will rephase, and a revival occurs, 
with the condition gT vcVnai (^J \\\ 2 + 1 - |A|) = 2ir, giving #T re vivai = 47r|A|. 
Thus, the characteristic timescales are given by: 

^oscillation — TTT' -^decay — — i -^revival — • (3.19) 

g\M g g 



and the associated behaviour is illustrated in Fig. 3.1 

Let us now consider how this behaviour may be approximated ana- 
lytically. The simplest approach might be to replace the sum over number 
states by integration, however this would inevitably lose the revivals, which 
are associated with the discreteness of the sum allowing rephasing. Thus, 
it is necessary to take account of the discreteness of the sum, and hence the 
difference between a sum and an integral. For this, we make use of Poisson 
summation which is based on the result: 

Y J 5(x-m) = Y J e 2 ™" 1 => E = E / dxe 2mrx f(x). (3.20) 

m r m r 

This allows us to replace the summation by integration, at the cost of 
having a sum of final results, however in the current context this is very 
helpful, as the final sum turns out to sum over different revivals. Applying 
this formula, we may write P cx ^ \ [1 - £ r 3fi{A(r, |A|,t)}] where 

A(r, |A|,i)= f — . ^ X exp 2irirx ttth + itjt ( \\\ -\ — ttt ttt? I 

1,1 M J ^/2tt\W Y 2|A| 2 y V 2|A| 8|A|V 

(3.21) 



3.3. MANY MODE QUANTUM MODEL 



- IRREVERSIBLE DECAY3 




Figure 3.1: Collapse and revival of Rabi oscillations, plotted for 
A = V200. 



Completing the square yields: 



A(r,\X\,t) = 7 / dxexn 

11 1 ! 



l L igt 

Wv + W\ 



x 2 + i I 2irr + ) .r 



Jgt\X\ 



: exp 



x ^ dx exp 

e igt|A| 

Vl + i5*/4|A| 



|A| 2 (27rr + gt/2|A|)- 
2 l + ist/4|A| 



2|A| : 



1 + 



4|A| 



x — i\X\ 



exp 



lAl 2 (27rr + gt/2|A|r 
2 l + i 9 t/4|A| 



2|A| 



2 27rr + gt/2|A| V 
l + ^/4|A| J 



(3.22) 



From this final expression one may directly read off: the oscillation fre- 
quency <?|A|; the characteristic time for the first collapse 2y/2/g (by taking 
r = and assuming gt/\X\ is small); the delay between successive revivals 
47r|A|/g; and finally one may also see that the revival at gt = 4-7rr|A| has 
an increased width and decreased height compar ed t o the initial collapse, 



given by a factor y/l + 7r 2 r 2 as is visible in figure 3.1 



3.3 Many mode quantum model — irreversible 
decay 

So far we have discussed a two level system coupled to a single mode, with 
light described either classically or quantum mechanically. One important 
difference between these situations is that when described classically, there 
was stimulated absorption and stimulated emission (hence the periodic os- 
cillation), but if the classical field had vanished, then there would have been 
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no transition. For the quantum description, spontaneous emission also ex- 
isted, in that the state |0, f) hybridises with the state |1, j), due to the +1 
in the \Jn + 1 matrix elements. We will now consider the effects of this +1 
factor when we couple a single two level system to a continuum of radiation 
modes. 

Our Hamiltonian will remain in the rotating wave approximation, but 
no longer restricted to a single mode of light, and so we have: 



H = eS z + ^4,n«k,n + 

n,k k 



gk,T. 
~2~ 



ats + H.c. 



where u>k = ck and 



gk,n _ ee k;?1 • dfc a 

2 " y/2eou k V ' 



(3.23) 



(3.24) 



Wigner— Weisskopf approach 

To study spontaneous emission, we start from the state |0,|), in which 
there are no photons. The subsequent state can be written as: 



a(t)|OT) + X)Ac,n(t)|lk,ni), 



(3.25) 



where |lk,n) contains a single photon in state k,n. Then, writing the 
Schrodinger equations as: 

id t a = |a + ^ ^pAc, n , id t Pk, n = (-| + ^fc) /?k,n + ^pa, (3.26) 

k,n 

one may solve these equations by first writing a = ae" !ei / 2 , (3 k = /3j,e*( e / 2_a;fc )* , 
and then formally solving the equation for with the initial condition 
p k (0) = 0, to give: 



i 9 ^ J' dt'e^-^'a(t'). (3.27) 



Substituting this into the equation for a gives: 



d t a = - 



k,n 



9k,n 



J(oj k -e)(t'-t)^ t l 



1 dt > [ ^r(u,) e ^- e >c-«>a(f) 
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where we have defined: 



I» = 2tt^ 

k,n 



I ffk,n 
l~2~ 



<5(w - w kj „) 



(3.28) 



(3.29) 



Before evaluating T, let us find the resultant behaviour, under the assump- 
tion that T is a smooth function, meaning that T does not significantly vary 
over the range e ± V, i.e. dT/duj <C 1. In this case (which corresponds to 



3.3. MANY MODE QUANTUM MODEL 



- IRREVERSIBLE DECAY5 



the Markov approximation, meaning a flat effective density of states) the 



integral over u in Eq. (3.28) becomes a delta function, so that: 



d t a = - J* dt'T(e)5(t' - t)a(t') = -^a(t) (3.30) 

Thus, the probability of remaining in the excited state decays exponentially 
due to spontaneous decay, with a decay rate T(e) for the probability. The 
approach used in this section is known as the Wigner-Weisskopf approach. 

Effective decay rate 

We have not yet calculated the effective decay rate appearing in the previous 



expression. Let us now use the form of in Eq. (3.24) to evaluate T(e). 
Writing the wavevector in polar coordinates, with respect to the dipole 
matrix element d a j> pointing along the z axis (see Fig. |3.2[), we then have: 



r(e) = 2vr^3 J!! d<pde S m9k 2 dk6(ck-e)^^Y,^-^ 2 - ( 3 - 31 ) 




Figure 3.2: Directions of polarisation vectors in polar coordinates 



From Fig. 3.2 it is clear that z • = and one may see that z • 
eki = s'mO. For concreteness, the k direction and polarisations can be 
parametrised by the orthogonal set: 



sin o cos ( 
sin 8 sin ( 
cos 6 



eki 






f — sin^ \ 


, e k2 = 


cos <fi 




V o ) 



Putting all the terms together, one finds: 



e 2 \d ab \" 
(2^) 2 2e 



dO sin 3 



ckcdk 



5(ck 



2 J, 



1 4|tU| 2 e 

47T£o 3c 3 



(3.32) 



(3.33) 



26 



LECTURE 3. JAYNES CUMMINGS MODEL 



This formula therefore gives the decay rate of an atom in free space, asso- 
ciated with the strength of its dipole moment and its characteristic energy. 



Further reading 

The contents of most of this chapter is discussed in most quantum optics 
books. For example, it is discussed in Scully and Zubairy [10], or Yamamoto 
and Imamoglu 



Questions 

Question 3.1: Collapse and revival 

Consider the "linear" model: 

H = uja)a + ec^c + g{a)c + c ] a). (3.34) 

What is its spectrum? If one starts from a coherent state of photons, 
exp(Aat)|0, 0), what is the expectation of occupation of the "atomic" mode 
c after time t. Explain why there is no collapse. 
Consider now the model [101 Sec. 6.1]: 

7F7 v 



H= L -[ 9aVa>a Wa. (3.35) 

2 y gva'aa' —e J 

Once again, find the time evolution of the probability of being in an excited 
state, and thus find the timescales for oscillation, collapse and revival in 
this case. [Note that this problem is exactly solvable, no approximations 
should be required] . 

Question 3.2: Noisy classical driving 



Let us consider the problem in Eq. (3.1) perturbatively; i.e. to leading 
order in g\; for non-degenerate perturbation theory to be valid here, we 
must for the moment assume we are away from resonance. Show that in 
first order time-dependent perturbation theory, the probability of being in 
the excited state initially grows as: 

P CX (T) = T J drX(t + r)X*(t), X{t) = (3.36) 

[Hint: Find the standard time-dependent perturbation result for the am- 
plitude in the excited state, then rewrite the amplitude squared to bring 
out a factor of time delay T.] 

Now, consider the case where in addition to the desired driving there is 
some random field, e JW °* — > e tu) o t + r i{t) _ jf j s a Gaussian random variable 
with (rj(t)) = 0, and (rj(t)rj(t')) = T2<5(t — t'), show that the transition rate 
now becomes 

,^^>^_EW (3 ,7) 



2 ( e -u, ) 2 + r 



2 
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3. A Further properties of collapse and revival 

When considering collapse and revivals of the Jaynes-Cummings oscilla- 
tions, we discussed only the probability of finding the atom in an excited 
state. We can however understand more of where this collapse and re- 
vival originates from by considering the entanglement between the two- 
level system and the photon field. The entanglement is defined in terms of 
S = —Tr(pTLS^ n PTLs)i where ptls is the reduced density matrix of the 
two level system: 

ms = X>h/>><#i>. (3-38) 



summing over e.g. number states of the photon field. Using our previous 
results for Rabi oscillations, the full state of the photon field can be written 
as: 



W)) 



n 



(i)|n-l,T) + 5>n(*)M). 



with (restricted to the resonant case): 



a n (t) 



A" 



-|A| 2 /2 



-isin {g^/nt/2) 
cos (gy/nt/2) 



In terms of these, the reduced density matrix will have the form: 
Kl 2 a* 1 (3 n \_fl/2 + Z X + iY 



Ptls 



+ 1 



X-iY 1/2 -Z 



(3.39) 



(3.40) 



(3.41) 



For a two dimensional density matrix, since the trace must be unity, by 
normalisation, the eigenvalues (and hence the entanglement entropy) are 
entirely determined by the determinant. In terms of the above parametrisa- 
tion, the determinant is 1/4-R 2 where R 2 = X 2 +Y 2 +Z 2 . If \R\ = 1/2, the 
determinant is zero, and so the eigenvalues are 1,0, i.e. this is a pure state, 
with no entanglement. If \R\ = 0, then the determinant is 1/4, and the 
eigenvalues are 1/2,1/2, with maximum entanglement. Thus, by plotting 
1/4 — R 2 we have a non-linear parametrisation of the entropy. The vector 
R is the Bloch vector parametrisation of the density matrix, something we 
will make much use of in subsequent lectures. 



From the form of the diagonal matrix elements, and the sum in Eq. (3.15 ) 
we have that 1/2 + Z = P cx and so 



Z 



2 ^ 



y/l+igt/4\\\ 



exp 



lAl 2 (27rr + gt/2|A|r 
2 l+igt/4\\\ 



(3.42) 



This is non-zero only near the revivals, so if we consider what happens away 
from the revivals (which is the majority of possible times), then we have 
|i?| 2 = X 2 + Y 2 , which can be written as: 



\R\ 



2n n\e^ Z 



A 



2^n + 1 



sin | ^(a/ti + Vn + 1) j +sin (^-(y/n- \Jn + 1) 



(3.43) 
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In this expression, the sum and difference sine terms can be treated sepa- 
rately. The first term, depending approximately on s'm(gt^/n) will behave 
very similarly to the excitation probability, with the only difference being 
this is the imaginary part, rather than the real part, of the complex expres- 
sion. Thus, once again, away from revivals, this term is small. The last 
term is however different. Expanding for large n, one has: 

\R\ ~ ] 
1 

~~ 2 

where the last line replace the sum by an integral (not worrying about re- 
vivals of this discrete sum). This looks superficially similar to the excitation 
probability, but the timescales are different — the oscillation timescale here 
is t ~ 4| A| /g; this is the same timescale as the revivals of the excitation prob- 
ability. Such a result makes sense — collapse of the excitation probability 
occurred because of dephasing, i.e. decoherence, arising from entanglement 
between the two level system and the photon field. The revival timescale 
is thus connected to the period of the entanglement oscillations. 
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Figure 3.3: Oscillations of the entanglement (top panel) and its 
connection to collapse and revivals of the excitation probability 
(bottom panel) in the case of A = 30. 
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Lecture 4 



Two level systems coupled 
to many modes: Density 
matrix description 

The aim of this lecture is to repeat the physics discussed at the end of 
the last lecture — i.e. the relaxation of a two-level system coupled to a 
continuum of radiation modes — but to present it in the context of the 
density matrix of the two-level system. This framework will then allow 
us to consider other relaxation and decoherence processes, such as those 
arising due to noise sources, or inhomogeneous broadening in an ensemble 
of atoms. We will then illustrate the use of this approach by discussing the 
behaviour of a two-level system illuminated by a coherent state of light, but 
able to radiate into a continuum of modes. In this lecture we will study the 
steady state of this system, and in subsequent lectures we will look at the 
spectrum of emitted light, which will require knowing more about the full 
time-evolution of the system. The formalism we develop here will be useful 
also when we later consider the laser, and its more quantum mechanical 
variants of single atom lasers and micromasers. 



4.1 Density matrix equation for relaxation of 
two-level system 

We begin by considering again the problem of a two-level system coupled to 
a continuum of radiation modes, which in the previous lecture we solved via 
the Wigner-Weisskopf approach, of eliminating the Schrodinger equations 
for the state of the radiation field. Let us now apply the same idea in a 
more formal scheme, by eliminating the the behaviour of the continuum 
modes in order to write a closed equation for the density matrix of the 
system. If we define the system density matrix as the result of tracing over 
all "reservoir" degrees of freedom, we can write an equation of motion for 
this quantity as: 

^-ps = -iTv R {[H SR ,p SR }}. (4.1) 
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in which S means state of the system, R the reservoir (in this case, the 
continuum of photon modes), and Ttr traces over the state of the reservoir. 
After removing the free time evolution of the reservoir and system fields, the 
remaining system-reservoir Hamiltonian can be written in the interaction 
picture as: 

H SR = E| (<r-ale**-* + * + a„jt~**) . (4.2) 



In general Eq. (4.1) will be complicated to solve, as interactions between 
the system and reservoir lead to entanglement, and so the full density ma- 
trix develops correlations. This means that the time evolution of the system 
would depend on such correlations, and thus on the history of system reser- 
voir interactions. However, the same assumptions as in the previous lecture 
will allow us to make a Markov approximation, meaning that the evolution 
of the system depends only on its current state — i.e. neglecting memory 
effects of the interaction. 

For a general interaction, there are a number of approximations which 
are required in order to reduce the above equation to a form that is straight- 
forwardly soluble. 

• We assume the interaction is weak, so psr = ps <8> PR + SpsR, where 
8p S R ^ 0(g k ). Since p s = Ttr(psr), it is clear that Tv R (Sps R ) = 0, 
however the result of Tr R (Hs R Sps R ) may be non-zero. Such corre- 
lations will be assumed to be small, but must be non-zero in order 
for there to be any influence of the bath on the system. In order 
to take account of these small correlation terms, generated by the 
coupling, we use the formal solution of the density matrix equation 
PSR{t) = -ifdt'[H SR (f),ps R (t!)} to write: 



d_ 
dt 



Ps (t) = -Ttr {J*dt' [H SR (t), [H S R(t'),p S R(t')]] | . (4.3) 



We may now make the assumption of small correlations by assuming 
that ps R (t') = ps(t')®PR(t')} an d that only the correlations generated 
at linear order need be considered. 

We further assume the reservoir to be large compared to the system, 
and so unaffected by the evolution, thus p R (t') = p R (0). Along with 
the previous step, this is the Born approximation. 

Finally, if the spectrum of the bath is dense (i.e. the spacing of energy 
levels is small), then the trace over bath modes of the factors e lUkt 
will lead to delta functions in times. This in effect means that we 
may replace ps(t') = Ps(t)- This is the Markov approximation. 

One thus has the final equation: 

j t p s (t) = -Tr R Udt' [H SR (t), [H SR (t'),p s (t)®p R (0)))\. (4.4) 
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To evaluate the trace over the reservoir states, we need the trace of 
the various possible combinations of reservoir operators coming from Hsr. 
These will involve pairs of operators a k ,a\. Using cyclic permutations, 
these trace terms can be written as thermal average {X) = Tr r(prX). For 
bulk photon modes, the relevant expectations are: 

[ a k a k>) = (44') = °' ( a l a k>) = Skk'nk, (afe4') = 5 kk'(n k + 1). 

(4.5) 

In these expressions we have assumed nk / 0, i.e. we have allowed for 
a thermal occupation of the photon modes. The results of the previous 
lecture assumed an initial vacuum state, and thus nk = 0. 

With these expectations, we can then write the equation of motion for 
the density matrix explicitly, writing = exp[i(ujk — e)(i — t')], to give: 

J t Ps(t) = |y J q dt' [{cr~a + n k £ k + <J + <J~{n k + l)C k ) ps 

- (a~pa + (n k + 1) + a + pa'n k ) (£ fc + Ck) 

+p s (a-a+rikCk + + 1)&)] • (4.6) 

In this expression, the first term in parentheses comes from the expres- 
sion Tr^ [HsR(t)Hsii(t')psR], the second comes from the conjugate pair 
TrR[H S R(t)psRH S R{t')} and t <-»■ t' [hence (£ fe + Q)], and the final term 
comes from Tr^ [psRH S R(f)H S R(t)]. 

At this point, we can now simplify the above expression by using the 
assumed fiat density of states of the photon modes, so that just as in the 
Wigner-Weisskopf approach, we may write: 

£ I f 2 a = m - *'), El? 2 = r ^(* - o- (4-7) 



Here we have assumed that not only is the density of states flat, but that 
also n k is sufficiently flat, that we might write an averaged occupation n. 
Such an approximation is only valid at high enough temperatures — we 
will discuss this in some detail in a few lectures time. 

d r 

-J t Ps{t) = - 2 n { a ~ a+ Ps ~ 2v + Ps°~ + Ps<r~° + ) 

- ^(n + 1) {a + a-ps - 2a~ p s cr + + ps^^) (4.8) 

The factor of 1/2 here comes from the regularised half integral of the delta 
function. The two lines in this expression correspond to stimulated absorp- 
tion (which exists only if n > 0), and leads to excitation of the two-level 
system, and to emission (both stimulated and spontaneous). If n = 0, only 
emission survives, and the results of the previous lecture are recovered for 
the relaxation to equilibrium. 

One can see the behaviour from this equation most clearly by writing 
the equations for the diagonal components: 

r_ r 

dtPn = -2™( _2 ^U) - 2 -(^ + 1 )( 2 PTt) = -dtPU (4-9) 
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Then, using p^ = 1 — p^, one may find the steady state probability of 
excitation is given by: 

= r[n(l- m )-(n + l) m ] - m = __. (4.10) 

This formula is as expected, such that if one uses the Bose-Einstein occupa- 
tion function for the photons, n = [exp(/3e) — l] -1 , the excitation probability 
of the two-level system is the thermal equilibrium expression. 



Equations of motion with coherent driving 

As a slightly less trivial example of the effect of the bath, let us now turn 
to the combination of coherent driving and decay. We will assume from 
hereon that the reservoir photon modes are empty 1 . The evolution of the 
density matrix is then controlled by: 

d t p = -i [H, p]-^ {a-a + Ps - 2a + p s a~ + p s a~ a + ) (4.11) 

where H describes the coherent, Hamiltonian dynamics, which we will take 
to be given by: 

H = ea z + ^e~ iu)t a + + H.c (4.12) 

One may then write out the equations of motion for each component ex- 
plicitly, giving 

re got _ iu}f go. ju}f el 

o t pn = -i [-ptt + Y e ~ puy c -Pn2\~ T Pih I 4 - 13 ) 



d tPU = ~ 



^Pn + Y Pu~PnY e ~ Pn \2 



r 

9 

(4.14) 



2 PU- 



These two equations are sufficient, since the unit trace property of a density 
matrix implies pn = 1 — p^, and the off diagonal elements are complex 
conjugate of one another by the Hermitian structure of the density matrix. 
To solve these equations, it is convenient to first go to a rotating frame, 
so that one writes = Pe~ tujt , and secondly to introduce the inversion, 
Z = (/9ff — Pn)/2. In terms of these variables, one has: 

dt Z = -i 9 -^(P*-P)-r(z+ 1 -) (4.15) 



2 v ' \ 2 , 

d t P = i(oo-e)P + i 9 ^2Z - -P. (4.16) 

To further simplify, we may separate the real and imaginary parts of P = 
X + iY, and write A = to — e for the detuning of the optical pump, so that 
one has the three coupled equations: 



1 This limit is frequently relevant in cavity quantum Electrodynamics experiments, 
since the confined photon modes start at energies of the order of feV~ 10 4 K, so the 
thermal population of such photon modes is negligible. 
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8 t X = -AY - (4.17) 
d t Y = +AX + gaZ - (4.18) 
d t Z = -gaY-v(z+^\. (4.19) 

These are Bloch equations for the Bloch vector parametrisation of the den- 
sity matrix (in the rotating frame): 

The three terms, A, ga, T correspond to: rotation about the Z axis (i.e. 
phase evolution); rotation about the X axis (i.e. excitation — nb the ob- 
servation that the rotation rate is ga, so a duration gat = n would lead to 
the inverted state); and T causes the length of the Bloch vector to shrink. 

As noted in the appendix to lecture |3j for a two level system, the eigen- 
values are directly given by the determinant of the density matrix. Hence, 
the entanglement entropy, describing the decoherence of the system due to 
its entanglement with the reservoir, depends monotonically on the length 
of the Bloch vector, X 2 + Y 2 + Z 2 . 



Before solving Eq. (4.17-4.19), we will first consider other kinds of de- 
phasing that may affect the density matrix evolution, leading to a gen- 
eralised version of these equations, allowing for both relaxation and pure 
dephasing. 



4.2 Dephasing in addition to relaxation 

In addition to relaxation, where excitations are transfered from the sys- 
tem to the modes of the environment, pure dephasing is also possible, in 
which the populations are unaffected but their coherence is reduced. In 
considering this kind of dephasing, it is useful to also broaden our view 
from considering a single two level system, to describing measurements on 
an ensemble of many two level systems. For the moment we will assume 
that even if there are many two level systems, they all act independently. 
As such, the expectation of any measurement performed on this ensemble 
can be given by: 

(X) = Y J Tr( Pl X) =Tr(pX), where p=^Pi. (4.21) 

i i 

We thus have two types of decoherence one may consider: 

• Broadening within a two-level system, coming from time dependent 
shifts of energies (and possibly also coupling strengths). A variety of 
sources of such terms exist; examples include collisional broadening 
of gaseous atoms, where the rare events in which atoms approach 
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each other lead to shifts of atomic energies. In solid state contexts, 
motions of charges nearby the artificial atom can also lead to time 
dependent noise. 

• Inhomogeneous broadening, where the parameters of each two-level 
system vary, and so off-diagonal matrix elements, involving time de- 
pendent coherences, are washed out in p. This is particularly an issue 
in solid state experiments, where the two-level systems may involves 
the dimensions of a fabricated or self- assembled artificial atom, and 
these dimensions may vary between different systems. 

To describe both such effects together, we consider adding noise, 

e^e + r](t). (4.22) 
This noise will then lead to a decay of the off diagonal correlations, via: 

P(t) -> P(t) (exp (i 7itf)dA \ (4.23) 

We will take (r](t)) = 0, and two distinct limits for the correlations of r/. 



Fast noise limit 



In order to extract a tractable model of time-dependent noise, the simplest 
limit to consider is that r/(i) has a white noise spectrum, i.e. it is Gaussian 
correlated with: 

(V(t)v(t'))=2 l0 5(t-t'). (4.24) 

This assumption is equivalent to assuming that the noise is "fast" , i.e. that 
any correlation timescale it does possess is sufficiently short compared to 
the dynamics of the system that it may be neglected. 

With this assumption of instantaneous Gaussian correlations, one may 



then expand the exponential in Eq. (4.23) and write: 



Pit) -P(i) (l + if* dt'vV) ~ \ ff dt'dt"r,(t f ) V (t f ) + ...\ 
=P(t) 



270 



2! ./o 4! 6 



t \ 2 

dt' | + 







P{t) e -~< 0t 

(4.25) 



with the combinatoric factors coming from the number of possible pairings 
in the Wick decomposition of {rf 1 ). 

This noise thus leads to an enhanced decoherence rate for the off di- 
agonal terms, so that T/2 should be replaced by I/T2 = (r/2) + 70 in 
the density matrix equations. One may equivalently define the relaxation 
time \ jT\ = T. In general T2 < T±, as dephasing is faster than relaxation, 
however the factors of two in the above mean that the strict inequality is 
T2 < 2Ti, with equality holding when there is only relaxation. 
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The Bloch equations are thus given by: 



d t X 



-AY-—X 
T 2 



8 t Y = +AX + gaZ - —Y 
dtZ = -gaY - i-fz +\ 



(4.26) 
(4.27) 
(4.28) 



Static limit (inhomogeneous broadening) 

To describe inhomogeneous broadening, rj should represent the system- 
dependent variation of energies, and should have no time dependence. The 
same behaviour is also relevant if time-dependent noise terms vary slowly 
on the timescale of an experiment. (For example, slow noise may arise 
from a nearby impurity have multiple stable charging states, so that the 
charge environment can take multiple different values; if the charging and 
discharging of this impurity is slow compared to the experiment, this energy 
shift will effectively give inhomogeneous broadening when averaged over 
multiple experimental shots.) 



For this case rj is static, but randomly distributed, and Eq. (4.23) be- 
comes: 



P(t) -> P(t) / dr]p(r))e ir,t . 



(4.29) 



If we consider a Lorentzian probability distribution for rj with width 71, i.e. 
p(rf) = 271/(7! + rj 2 ), then one finds P(t) — > P{t) exp(— 7i|i|). 

If one has both noise terms as in the previous section, and inhomoge- 
neous broadening, one may then distinguish T 2 as defined above, and T|, 
where 1/T 2 = (r/2) + 70 + 7i- 



Distinguishing T2 and T| 

Since both 70 and 71 lead to decay of the coherence functions, a simple 
experiment measuring coherence will see the lifetime T£. However, the 
dynamics of each individual two level system remains coherent for the longer 
timescale T 2 . Such extra coherence can easily be seen in any sufficiently 
non-linear measurement. An example of how the difference of T 2 and T| 



can be measured is given by photon echo, illustrated in Fig. 4.1 




Figure 4.1: Cartoon of photon-echo experiment, which removes 
dcphasing due to inhomogeneous broadening, leaving only true 
dcphasing rate to reduce intensity of revival. 
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From the Bloch equations, Eq. (4.26-4.28), it is clear that a resonant 
pulse with duration gat = ir/2 brings the Bloch vector to the equator. 
From there, inhomogeneous broadening means that (in a rotating frame at 
the mean frequency) the Bloch vectors for each individual two-level system 
start to spread out. However, by applying a second pulse with gat = it, 
a 7r rotation about the X axis means that whichever Bloch vector spread 
at the fastest rate is now furthest from the Y axis, and the subsequent 
evolution sees the vectors re-converge. At the final time, the coherence of 
the resulting pulse will have been reduced by 70 and T, but 71 will have 
had no effect. 



4.3 Power broadening of absorption 

Having developed this formalism, we may now use this to re-examine the 
question we started the previous lecture with — what is the effect of shining 
classical light on a two-level atom, but now accounting for decay. Whereas 
in the previous lecture, the absence of decay meant that the result was 
Rabi oscillations, with decay a steady state is eventually found. Looking 



for steady state solutions of Eq. (4.26-4.28) one finds from the equations 
for X and Z that: 

Substituting these into the equation for Y, and recalling P ex = Z + 1/2, 
gives: 

ga 



= P O x 

After re- arranging, one then finds: 



a AT 2 1 1 

A — 7f + 7f — 
gali ±2 gali 



2 



(4.31) 



1 TiT 2 {gaf = 1 (Ti/T 2 ){gaf 

ex 2(T 2 A) 2 + l + T 1 T 2 (<?a) 2 2 A 2 + [1 + TiT 2 (ga) 2 ] /Tf 1 ' 

The addition of damping to the equations of motion has thus had a 
number of important consequences: 

• It is responsible for a steady state existing at all (but note that this 
steady state is in a frame rotating at the pump frequency u, so the 
physical coherence is time dependent.) 

• It gives a finite width to the resonance, even at weak pump powers 

- in the Hamiltonian case, as ga — ► 0, the width of the resonance 
peak vanished. 

• It modifies the overall amplitude at resonance, by a factor depending 
on the power, such that at ga — * 0, there is no response. 

One feature of this absorption probability formula that can be remarked 
on is that the linewidth depends on the intensity of radiation. This is 
a consequence of the nonlinearity implicit in a two-level system. Such 
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dependence is not so surprising given that in the Hamiltonian dynamics, 
the resonance width already depended on the field strength. One may note 
that if one considered an harmonic atomic spectrum, rather than a two-level 



atom, no such broadening would be seen (see question 4.2). 




Figure 4.2: Power broadening of absorption — probability of 
excitation vs atom-photon detuning for increasing field strength. 



4.4 Further reading 

Once again, the contents of this chapter can be found in most standard 
quantum optics books, e.g-JTUl H2] . The discussion of the resonance fluo- 
rescence problem in particular can be found in chapter 4 of Meystre and 
Sargent III [T2|. 

Questions 

Question 4.1: Fast noise via density matrix equations 

As another way to consider fast noise terms, consider interaction be- 
tween the two-level system and a bath of radiation modes given by the 
interaction picture coupling: 

H SR = ^Cn (B n e- iuJnt + 4e^*) (4.33) 

n 

where B n are the quantised modes of a field whose field strength shifts the 
energies of the two-level system. 

4.1. (a) By making the Born-Markov approximation, show that this leads 
to a density matrix equation of the form: 

^Incise = "70 (f-2aV z ), (4-34) 

and determine how 70 depends on the temperature of the reservoir if the 
reservoir follows a thermal Bose-Einstein distribution. 
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4.1. (b) Show that this describes pure dephasing, i.e. it has no effect on 
the diagonal elements of the density matrix, and leads to a decay of the 
off-diagonal terms dtPn = — 7oPtJ.- 

Question 4.2: Absence of power broadening for an harmonic os- 
cillator 

Consider an harmonic oscillator mode H sys = etfb, coupled to a contin- 
uum of radiation modes via the interaction picture Hamiltonian: 

H SR = y (bale 1 ^-^ + b^a^'^ . (4.35) 
k 

4. 2. (a) Show that (for an empty bath) this leads to the density matrix 
equation of motion: 

d tP = - E (V V - 2bptf + ptfb) (4.36) 

4.2.(b) Now consider adding the Hamiltonian dynamics describing pump- 
ing by a coherent field, with H = (ga/2)(b'e~ iu ' t + H.c.) + etfb. Determine 
the density matrix equation equation of motion in the Fock basis (i.e. num- 
ber state basis). 

4.2.(c) Show that the ansatz: 

1 ( -iga/2 \ n f iga/2 \ m 

satisfies this equation of motion. 

4.2.(d) From the normalised density matrix, find the average excitation 
probability, and show that there is no power broadening. 



Lecture 5 

Resonance Fluorescence 



At the end of the last lecture, we found the steady state excitation proba- 
bility for a two-level system pumped by a coherent light field. This steady 
state emerged due to the competition between coherent pumping, and in- 
coherent decay into the continuum of radiation modes. The aim of this 
lecture is to discuss the spectrum of this fluorescence, when the atom is 
driven near resonance. This is intended both as an illustration of applying 
the density matrix equation of motion approach to more complex problems, 
and also to reveal further aspects of the behaviour of a two-level system. 



5.1 Spectrum of emission into a reservoir 

Our aim is to calculate the spectrum of the emission into the continuum 
of photon modes. Formally, the spectrum of radiation is given (via the 
Wiener-Khint chine theorem, applicable to a stationary process, which ours 
is) by: 



/oo 
e iut {E + (t)E-(0)) =2K 
-oo 



e tut (E + {t)E-{0)) 



(5.1) 



where we have divided: 



E(t) = E + (t) + E-(t), E-(t) =Y, £ ka k , n (t)e n , k , (5.2) 

k 

with £}, = y^LOk/^eoV . The first challenge is to write the time dependent 
reservoir operators in terms of system operators. If we work in the Heisen- 
berg picture, then this is straightforward, the interaction Hamiltonian was: 

Hsn = E 9 i {^^^ + <r + « fc e <(e ~ h,fc) ') • (5.3) 
k 

and so the solution to the Heisenberg equation for a^t) can be written as: 
M*) = f fdt'e^-W-Va'it') + o fc (-oo). (5.4) 
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We will neglect second term dfc(— oo), since we assume the state at t = — oo 
to be a vacuum, i.e. the only photons are those coming from the atom. In 
this case we may write: 

E~(t) = fdt'^^^~^'-%-(t')e n ^ (5.5) 

k 

Because (at least near resonance) £ ^ oc the sum in the above expression 
is Markovian as in the previous lecture, and so we may write E~ (t) oc a~ (t) , 
hence up to an overall constant, we have: 



I(v) oc 23? 



/•oo 

Jo 



(5.6) 



This reduces the problem to that of finding a two-time correlation func- 
tion of the open quantum system. This is however something we have not 
yet calculated — the density matrix completely describes single time cor- 
relation functions, but in general more knowledge is required to find the 
correlation of operators at two different times. However, for Markovian 
baths it turns out that it is possible to relate this two-time correlation to 
the density matrix evolution via the quantum regression theorem. 



5.2 Quantum regression "theorem" 

The quantum regression theorem can be stated in terms of the time evolu- 
tion of single-time correlation functions, governed by: 

(B(t + T)) = Y,<Xi(T)(B i (t)), (5-7) 

i 

where Bi(t) is some complete set of operators, and the functions a^r) solve 
the averaged equations of motion for Then, the theorem states that 

one can write: 

(A(t)B(t + T)C(t)} = J2^(r)(A(t)B t (t)C(t)), (5.8) 

i 

Note that ordering is important, as A, Bi, C are non-commuting operators. 

In order to prove the theorem, it is first necessary to state more explore 
further what is implied by the Markovian approximation. 



Significance of the Markovian approximation 

Formally stated, a Markovian system is one where the future evolution is 
governed by the current state, and not by any history of the system. In the 
current context, current state means the state of the system, not including 
the baths. Thus, the assumption is that one can write an equation of 
motion for the system density matrix that depends only on the current 
value of the system density matrix. Therefore, the baths must have no 
memory of previous states of the system. 
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Generically, without the Markovian approximation, the equation of mo- 
tion of the system density matrix could be written as 



dp sit) 
dt 



i[H s ,p s (t)]+ f dt'T{t-t')p s {t'). (5.9) 

J —oo 



The term T{t — t') describes how the state of the bath influences the system 
at time t, and depends on the state of the system at the earlier time t'. The 
Markovian approximation is that T(t — t') = T5{t — t'). If one describes 
the bath by a sum over many independent modes, these delta correlated 
response functions imply a dense spectrum of bath modes; this is what was 



used in Eq. (4.7) 



There is another consequence of the Markovian limit, which was also 



used earlier in Eq. (4.4). The full evolution really depends on psn; validity 
of the Markovian approximation then requires that the state of the system 
is sufficiently described by ps = Ti R ps R . This means that correlations 
between the system and the reservoir must be unimportant, and so it is 
sufficient to write psn(t) = Ps(t) ® Pn(t)- It is this statement, as we see 
next, that implies the quantum regression theorem. 

Under what physical conditions is the Markovian approximation a good 
description of a real system? Firstly, to have truly delta correlated noise 
requires a fiat spectrum of the baths, as we have already assumed. However, 
a fiat spectrum is only required for the bath modes that couple almost 
resonantly to the system (i.e. \ujj — u\ < T). If one considers a reservoir 
which is not in its vacuum state — i.e. with some thermal occupation - 
then a Markovian description of the system requires that not only is the 
density of states fiat, but that we may also approximate: 



\9k 
I 2 



V^-^)* = nTS(t). 



i.e. we require that varies slowly for over the range of k for which \u>j — 
u\ < V; e.g. for a thermal distribution, we require T>I. This is clearly 
always an approximation; the consequences of this approximation will be 
investigated further in the next lecture. 

Proof of regression theorem 

The Markovian approximation implies that the evolution of single-time 
expectations can be written as follows, in terms of the unitary evolution 
U (r) of system and reservoir: 

(B(t + r)> = Tt s Tt r (u\T)B(t)U(T) Ps (t) p fl (t) ' 
= Tr 5 [B(t)Tv R (u(T)ps(t) ® PR {t)U\r) 
= Tr 5 \B{t)Tv R ( Ps {t + t)® p R (t + t))]. (5.10) 



Here we made use of the cyclic property of the trace, and the relation 
between evolution of density matrices and that of Heisenberg operators. In 
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an analogous way, the two-time expectation evolves according to: 
(A(t)B(t + r)C(t)) = Tt s Tt r (A(t)rf (T)B(t)U(r)C(t)p s (t) p R (tfj 



B(t)Tv R (U(t) [C(t)p s (t)A(t)]® PR (t)rt(T] 
B(t)Tv R ( [C Ps A] (t + t) ® p fl (t + r; 



(5.11) 



Thus, the only difference between the right hand sides of Eq. (5.101 and 



Eq. (5.11 ) is to replace the initial system density matrix ps with the product 



CpsA. However, in the Markovian approximation, we know that Eq. (5.10) 



can be solved using Eq. (5.7), for any initial density matrix. Thus, one can 
use this solution to give 



{A{t)B{t + r)C{t)) = Y J ^{r){B l {t)) 



Ps(t)^C(t) Ps (t)A(t) 

6 

Y^cxii^TrslB^C^ps^Ait)] 

i 

£ ai(r)Ti s [A(t)B t (t)C(t) Ps (t)] (5.12) 



which is the desired result in Eq. (5.12). 



Thus, to find a two-time correlation in practise, the scheme is this: find 
the general time-dependent solution of one of the operators, and how it 



depends on initial expectation values — find the explicit form of Eq. (5.7) 
and then insert the second operator(s) in these initial expectations. 



5.3 Resonance fluorescence spectrum 

From the quantum regression theory, we now have that: 



I{y) oc 23? 



e iut Tr [a + p(t)} 



p(t=0)=a-p 



(5.13) 



i.e. we should solve the density matrix equation of motion to find p(t) given 
the initial condition a~ po, where po is the equilibrium density matrix. 

Since the equation of motion for the density matrix is a first order 
differential equation, dtp = M[p], where M is a linear super-operator, we 
can easily solve by Laplace transforming to give: 



M[p(s)] - sp(s) 



-v PO, 



(5.14) 



where the right hand side encodes the initial conditions. From the Laplace 
transform p(s), the Fourier transform corresponding to the intensity spec- 
trum is then given by: 



I(u) oc 25ft {Tr [a + p(i 



-iv 



+ 0)]}. 



(5.15) 
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Let us now specify the various parts of this for the case of resonance 
fluorescence, i.e. A = uj — e = 0. WE will work in the rotating frame as in 
the previous lecture. The steady state condition then becomes: 



1 



Po 



2[1 + T^ga) 2 } 



T^ga) 2 
-iT 2 (ga) 



iT 2 {ga) 
2 + T 1 T 2 (ga) 2 



(5.16) 



(making use of X = and Y, Z as calculated in the previous lecture.) 
The important part of the spectrum however comes from the super- 



operator in Eq. (5.14). If we write density matrices as a vector: 



PtT PIT 

Pn pu 



( y°n \ 

PIT 

pu 
V pu J 



(5.17) 



then in this four dimensional space, the super-operator can be written as a 



4x4 matrix, following Eq. (4.15) 



M 



( 



\ 



-1/Ti 

-iga/2 
iga/2 
1/Tx 



-iga/2 
-l/T 2 


iga/2 



iga/2 


-l/T 2 
-iga/2 



One can thus find the full density matrix as: 



I(v) cx -2K 



/ 



(0 1 0)[M-(-«/ + + )l] 



-i 



\ 

iga/2 

-iga/2 
/ 








(5.18) 



\ 



TiT 2 {ga) 2 
V +iT 2 (ga) J 



. (5.19) 



When inverting the 4x4 matrix, there will in general be four poles. How- 
ever, because the density matrix evolution is trace-preserving, one of these 
poles will necessarily be at zero, corresponding to the conserved quantity. 
The corresponding eigenvalue will be the steady state (the trace-preserving 
property ensures that the steady state is non-trivial). Since we know that 
in the steady state there is a non-zero polarisation, the weight of this zero 
pole will not in general vanish. This means there is a delta-function peak in 
the spectrum. Such a peak describes the elastic scattering of pump photons 
into the radiation modes. 

The remainder of the poles will describe inelastic scattering. To find 
these poles, we should solve: 



= Det [M + iv] =iv[iv 



1 

T 2 



(5.20) 



This has the zero mode as discussed above, and three other poles: 
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The first of these will describe a broadened resonance near elastic scattering 
(remembering that this is in the rotating frame) . The second term describes 
(at large enough a), a resonance shifted by ±<x When the weights of these 



poles are calculated, as shown in Fig. 5.1, the central peak has twice the 
weight of the other two peaks. This characteristic shape is known as the 
Mollow triplet, and has a simple interpretation. 




Figure 5.1: Spectrum of resonance fluorescence, from Eq. (5.191. 



Interpretation: of the Mollow triplet 

The origin of the Mollow triplet is the effective energy levels of the two- 
level system in the presence of a coherent field. If one considers in general 
a two-level system coupled to photons, there are Rabi oscillations. In the 
frequency domain, these oscillations correspond to a splitting of the eigen- 
states resulting from hybridising the states with n excitations; i.e. the states 
|n — 1,|), \n, I) mix and lead to eigenstates split by g\fn (if on resonance). 

If the system can now decay, it will emit an excitation into the contin- 
uum, and undergo a transition to a state with one fewer excitation in the 
system. However, starting from one of the hybridised eigenstates with n 
excitations, it is possible to undergo a transition to either of the eigenstates 
with n — 1 excitations — non-zero matrix elements connect them all. Thus, 
in general there are four transition frequencies. 

In the limit of a strong classical field, all four of these matrix elements 



have the same strength, and moreover, as indicated in Fig. 5.2 the Rabi 
splitting becomes almost n independent when n > 1, so two of the fre- 
quencies will match. This then leads to the characteristic 1:2:1 ratio seen 



in Figure 5.1 



5.4 Further reading 

The discussion of the resonance fluorescence problem in particular can be 
found in chapter 16 of Meystre and Sargent III [12 . The discussion of the 
quantum regression theorem in this lecture is however more closely related 
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Figure 5.2: Characteristic behaviour at high excitation number, 



ith 



\a\ , demonstrating a triplet of transition frequencies, 



with the intensity ratio 1:2:1. 



to the density matrix based approach given in Scully and Zubairy [TO] ; The 
Onsager-Lax (quantum regression) theorem was introduced by Lax |13j . 



Lecture 6 



Quantum stochastic methods 
and limitations of the 
Markovian approximation. 



The previous two lectures, have introduced the density matrix equations of 
motion to describe a system coupled to a bath, and have shown how the 
Markovian approximation implicit in these equations of motion leads one 
to the quantum regression theorem. This lecture reviews the foundations 
of these methods for open systems. Two further formalisms are presented, 
which can be seen as the stochastic equations of motion for which the prob- 
ability evolution is that of the density matrix. In the latter of these — the 
Heisenberg-Langevin equations — the consequences of the approximations 
required for a Markovian equation of motion are more explicit; these are 
discussed in terms of the fluctuation dissipation theorem. 

We will consider just the problem of a system coupled to a single decay 
bath, which gave us the equation of motion: 

d , s k 



dr 



p(t) = — — n \ aa) p — 2a} pa + paxvj 



+ (n + 1) (a) ap - 2apa t + pa) a) . (6.1) 

6.1 Quantum jump formalism 



Consider evolution under Eq. (6.1) when n = 0, i.e. when only decay can 



occur. Then, after a short time St, one has: 
ps(t + St) = p s (t) - ^ (a)ap s + p s a^a^ + n5tap s a) + 0(5t 2 ) 

1 2~ a ^ a ^ P- s 2~ a ^ a ) ^ K ^ap s -j KSta^ + 0(5t 2 ). 

(6.2) 

These two contributions to the density matrix can be interpreted as the 
conditional density matrices that arise under "no photon loss" and "one 
photon loss" ; they are added together because of the trace over states of 
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the bath. The probability of these two final states can be found from the 
trace of the conditional density matrices, hence: 



-Pone photon loss = Tr \^nStap s a' 1 J = nSt {a" a) . (6.3) 

This is as one intuitively expects; the probability of losing a single photon 
in [t, t + St] is given by the rate of loss of each photon multiplied by the 
total number of photons. Thus, the evolution can be written in terms of 
the probabilities of transition to the (suitable normalised) 

(1 - nStaU/2) |V) , . , [bA) 

\ ' =U=L P = 1 - kSM\o) a\ib) 

k yJ{il)\l-K6tcta\il)) 

This quantum jump formalism can be understood as describing the loss 
process in terms of the environment repeatedly measuring the number of 
photons in the cavity, and hence decohering the states with different num- 
bers of photons. 

Example: Photon loss from a single-photon state 

The above procedure can be illustrated for a particularly simple case: if the 
initial state contains a single photon then from the density matrix equation: 

dt(n + k\p\n) = — ^(2n + k)(n + k\p\n) 



+ Ky/(n + l)(n + k + l)(n + k + l\p\n + 1) (6.5) 

it is clear that the initial condition Vra : (n + k\p\n) = is preserved except 
for k = 0, and so the only equations to solve are 

d t (l\p\l) = -d t (0\p\0) = -k(1\ P \1). (6.6) 

so the full solution is: 

/9 = e - Kt |l)(l| + (l- e -«*)|0)(0|. (6.7) 

In the quantum jump formalism it is easy to reproduce this result; at 
each time step we have that: 

I in\ e> — 

(6.8) 




A more complicated example is given in Question 6.1 for decay of a 
coherent state. In that case, since a coherent state is unchanged under 
the loss of a photon, one has the counterintuitive result that it is failure 
to observe a photon that leads to the density matrix evolving. This is 
less surprising if one considers that the state has an uncertain number of 
photons to start with, and failure to observe any photons escape indicates 
that the wavefunction has been projected by observation onto a state with 
fewer photons. 
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6.2 Heisenberg-Langevin equations 



Just as one can solve a closed quantum-mechanical problem in either the 
Schrddinger or Heisenberg pictures, quantum stochastic approaches can be 
constructed in either picture. The previous sections generally described 
a stochastic evolution of the wavefunction, this section instead describes 
stochastic evolution of the system operators. Starting once again from 
the same coupled bath and system, and working in the rotating frame as 



normal, one has the Hamiltonian of Eq. (4.2), 



(4.2) 



From this Hamiltonian, the Heisenberg operator equations of motion are: 
id = [a, H SR ] = £ C k b k e-^-^ 1 (6.9) 



ib k = [bk,Hsn] = (kae 



(6.10) 



To eliminate the bath degrees of freedom, one can integrate Eq. (6.10), and 



substitute it in Eq. (6.9). Thus, 



bk(t) = -Kk / dt'e^-^'ait') + 6 fc (0) 
Jo 



which, substituted into Eq. (6.9), gives: 
ft 



[ dt' 
Jo 



a{t') 



(6.11) 



. (6.12) 



Then, defining the second term in brackets as F a (t), and making the same 
Markovian approximation: 



du 
2^ 



kF(u) 



(6.13) 



as in Eq. (4.7), so that the first term in brackets becomes nS(t — t'), this 

1 



equation becomes: 



;K a + F a (t). 



(6.14) 



Here F a (t) is a stochastic operator; it has quantum mechanical commu- 
tation relations related to its definition in terms of b k ] but b k is a bath 
operator, and the state of the bath is random — drawn from a thermal 
ensemble, so b k has different expectations for each realisation of the bath. 

Preservation of commutation relations 

The operator nature of F a (t) is apparent if one considers the commutation 
relations: 



= Cfce- i ( w - w *K*-*'> = nS(t - t'). (6.15) 
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This non-commuting nature is essential to ensure the preservation of com- 
mutation relations in the time-evolution of the operators a(t). The solution 



of Eq. (6.14) can be written 



a(t) = a(0)e" K */ 2 + / dt' F a {t')e- K{t -^l 2 . 
Jo 



(6.16) 



Using this, one can write: 



a(t),at{t)] = \a(0),a\0) 



— Kt 



+ f dt' f dt" 
Jo Jo 



F a (t'),F^t")] e -«*+«(f+*")/2 j ( 6 . 17 ) 



in which we have assumed no correlation between the initial state of the 
system and that of the bath. Then, using [a, at] = 1 in the initial state, 



and the commutator in Eq. (6.15) this becomes: 
a(t),a f (i) 



e~ Kt + 



-Kt 



f dt' 

Jo 

+ [1 



Kt 



■Kt+Kt' 



(6.18) 



This shows that the existence of the non-commuting stochastic terms was 
essential to preserve the operator commutation relations. 



Finite bath occupation; other correlations 

To reproduce all the results of the density-matrix formalism, it is necessary 
also to have expressions for the correlation functions of the stochastic opera- 
tors F a (t) as well as their commutator. Given the form for the commutator, 
all correlation functions can be found from this and the anticommutator, 



k,k' 



-i(io—ijJk)t+i(LO—ujy)t' 



-i(u-uj k )(t-t r ) 



2vr \ 2 



//roth ( ^ ) 6(t- t'). 



(6.19) 
(6.20) 



The last line here makes the additional assumption of T > k, i.e. that the 
temperature is large compared to the bandwidth, with the bandwidth taken 
approximately equal to the imaginary width. This assumption is clearly 
necessary if the entire formalism is to be Markovian (i.e. to depend only on 
current state, rather than history), and is identical to the approximation 



made in Eq. (4.7). The validity and limitations of this last approximation 



are the subject of the next section. 
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6.3 Fluctuation dissipation theorem and the 
Markovian approximation 

The fluctuation-dissipation theorem relates the correlation function of an 
operator with its response function. For the photon operators considered 
here, the correlation function can be defined by: 

C(fi) = f°° drC(r)e mT , C(r) = \ ({a(t), a\t + r)}) . (6.21) 

The response function can be found by considering the response to a 
perturbation: 

5H = Y^ fne int a j + H.c. (6.22) 
n 

The response function a(O) is given by 

(a(t)) = Y,amfoe im . (6.23) 



The discussion below uses the Heisenberg-Langevin formalism to eval- 
uate the relation between these operators with and without the second 
Markovian approximation above. 

The response function is the same in either case, since it is averaged 
and so does not depend on the stochastic operator F a (t). In the following 
it will be necessary to distinguish between the bare operators a(t) and the 
gauge transformed versions a(t) = a{t)e~ lut . Thus, 

(ft + f) <*(*)> ^E^ -"* (6-24) 
n 

where we have used the Heisenberg equation following from the gauge trans- 
formed perturbation Hamiltonian. This yields 

hence a(tt) = - oS) - in/2). 

Next, consider the fluctuation correlation function at long times, so 
that the initial conditions are not involved. In this case, Eq. (6.16) should 
become: 

a (t) = e iut f dt'F a (t')e- K{t - tl)l2 (6.26) 

and so: 

„ — 'iuJT ft ft+T 

C(r) = — J ^ J J^{{FjfiX{f)})fT<^ 1f -^ (6.27) 

We now have two choices; we may either use the exact relation in 
Eq. (6.19), or the Markov approximation in Eq. (6.20). 
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Without second Markovian approximation 



Considering first Eq. (6.19), we have: 

CM = — /"* dt' -^coth ( ^) W'-t')-«Vt + r-t>-t»)/2 

2 J-oc J-oc J-oo 27r V 2 / 

' Bv\ e i(u-u)(t+r-t) 



2 7_oo 2vr \ 2 J [i(u> - v) + K/2][-i(uJ - v) + k/2] 

1 f°° dv (p v \ e~ iuT 

= 2 L UJ ( W -,)2 + («/2) 2 (6 - 28) 

hence it is immediately clear that: 

C <°> " (?) ( M - n) " + (./2)» " » th (?) 9[a(!1)l (6 ' 29) 
which is the required fluctuation dissipation theorem. 

With second Markovian approximation 



In contrast, if we use Eq. (6.20) we have that: 



C( T ) = — [ dt' [ T dt'k C oth(^]5{t' -t'^e^ 21 ^- 1 '-^/ 2 

2 J-oo J-oc V 2 / 



2 V 2 



KCOth ^ ^ e -«(2 t+ r-2t')/2 



coth ( e -i.|r[/a (6 . 30) 



2 V 2 

from which it follows straightforwardly that 

C(O) = coth f ^ 3f[a(fi)]. (6.31) 

Thus, use of the Markovian approximation means that the fluctuation 
dissipation theory fails, unless T ^> (O — u) > k. Such a result can also be 
seen to be a limit for the use of the quantum regression theorem. In that 
case, the high temperature (classical) answer is known as the Onsager theo- 
rem. The quantum extension of the Onsager theorem is an approximation, 
that derives from the Markov approximation in the dynamics 



6.4 Further reading 

The topic of Heisenberg-Langevin equations, and stochastic operator equa- 
tions is discussed briefly in Scully and Zubairy [10] , and more fully in Gar- 
diner and Zoller [14J The implicit dependence of the quantum regression 
theorem on the use of the Markovian approximation is highlighted in the 
papers of Ford and O'Connel pJil HB], and the response by Lax [T7]. This 



topic is addressed again in question 6.2 
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Questions 



Question 6.1: Comparison of quantum jump and density matrix 
equations of motion. Consider an initial coherent state, 

p = exp(-|a| 2 )e aat |0)(0|e a * a (6.32) 

following the equation of motion of the density matrix with n = 0, i.e. 



K 

2 L 



a l ap + pa 1 a — 2apa 



(6.33) 



6.1. (a) Show that this equation of motion can be satisfied by the ansatz 
p(t) =exp(-|a(t)| 2 )e a( * )at |0)(0|e Q * Wa (6.34) 

and find the equation satisfied by a(t). 

Now consider this question in the quantum jump formalism. If the 
problem can be written as: 



p(t) = J2 P n(t)\Mt))(Mt)\ 



(6.35) 



where \ip n ) corresponds to the state after having lost exactly n photons. 

6.1. (b) Show that \ip n (t)) = |?/>o(i)); (i.e. that losing a photon does not 
change the state). 

6.1. (c) Find P n (t), and verify that ^2 n P n (t) = 1, and hence show that 

P (t) = \Mt))(Mt)\- 



6.1. (d) Show that at time t, the non-Hermitian evolution leads to 

|Vo(i)) = ■^exp(- K ta t a/2)|^o(0)) 
where M is an appropriate normalisation. 



(6.36) 



6.1. (e) Prove that this result matches that in Eq. (6.34). 

Question 6.2: Fluctuation dissipation theorem, and quantum re- 
gression 

This extended question repeats the analysis in the lecture of the limit 
of validity of the Markov approximation, but for a damped quantum har- 
monic oscillator. It shows how the above results do not depend on other 
approximations that were involved, such as rotating wave approximation. 
Starting from the Hamiltonian, 



H 



muJ 2 , J?_ 

2 + 2m 



(Qj ~x) 2 + 



1±_ 
2ra,- 



(6.37) 
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One can easily derive the equations of motion: 

x = p/m p = -mwpX + m j u; 'j( ( lj ~ x ) 



Qj = Pj/ m J 



Pj = -mjuj^qj-x) 



Eliminating pj and p, one can write: 

Qj = ~Uj {Qj -x); x = -uj x + ^ — uij (qj 



x) 



(6.38) 



6. 2. (a) Show that the equation for qj(t) is satisfied by: 



Qj(t) ~ x{t) 



q°j cos(uji) H — sin(u;jt) 



/d£ cos[ujj(t — t')]x(t') 
-oo 

(6.39) 



Substituting this in the equation for x one finds: 

ft 



X + OJqX 



-2 f dt'j(t - t')x(t') + F(t) 



(6.40) 



where 



27(* - t') = cos M* " 0] 



w = £ 



J , ,2 



lJ.: 



m 



Pi 



q® cos(ujjt) + — - — sm(ujjt) 



Let us assume a flat spectrum, so j(t — t') = j5(t — t'). This implies 

-► / (6.41) 
m J it 



The correlations of F require knowledge of correlations of p° and 5°, the 
initial momenta and coordinates. Using the reasonable assumption of mi- 
correlated, diagonal (thermal) states, one has (jjPq®) = 0, and: 



2m j to j (q»q k >) = 2-^- = S jk (2n( Wj ) + 1) = S jk coth j . (6.42) 

6.2.(b) Using the above correlations and definition of F(t), show that: 

\ (F(t)F(t>) + F(t')F(t)) = £ ^ cos M* " H)\ ^th (^) 

j 

= / duj^- cos\uj(t — t')l coth ( — — — ) . 
Jo 7rm \2k B Tj 

(6.43) 
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6.2. (c) Show that the general solution to: 

x + jx + ujqX = F(t), 
can be written (using u)\ = — 7 2 /4) as: 



x(t) 



dt'e~^~^ Sin[uJl{t - t,)] F(t'). 



(6.44) 
(6.45) 



6.2. (d) Thus, show that: 



C(r) = - (x(t + r)x(t) + x(t)x(t + t)) , 



is given by: 
C(t) = 



du cos(cjt)— — ^-r — — 

Tim (ijj z — lOq) z + 7 z u> z 



Fourier transforming, this gives the correlation function 



/oo 
dte- iujt C{t) 
-oo 



7 



(6.46) 



coth U^J- (6 - 47) 



> o.o , > o ''-"li I .,. .. ) • (6.48) 
(oj z — uj z ) z + j z lo z m \2k B T J 

Let us now relate this fluctuation correlation function to the response 
function. To define the response function a(t), consider the averaged equa- 
tion motion, responding to a driving force f(t) = e luJt fo: 

fit) 



(x) +j(x) +u (x) 



m 



(6.49) 



6.2. (e) Show that the response function, (x) u = a(u;)/ w is given by: 

1 1 



a(u>) 



m luq — (J 1 + i"fuj 



(6.50) 



Hence, show one can write the fluctuation-dissipation theorem as: 

C(u,) = 3[a(u;)]cothf^Y 

Now, let us repeat the above calculations with the Markovian approxima- 
tion. 



(6.51) 



6.2. (f) Show that the general time-dependent behaviour of (x(t)), fol- 



lowing Eq. (6.49) is given by: 

(X(t)) = 



-jt/2 



(xq) [ cos(uii) — sinfu^i) ) + — — sin(cJit) 

2u>i J mu>i 



(6.52) 



Then, using the Markovian approximation for quantum regression, one can 



find C(t) as defined in Eq. (6.47), by replacing the initial density matrix 



with pit) — > x(t)p(t). Assuming equilibrium single-time correlations, one 
has (x(t)p(t) + p(t)x(t)) = 0, and {x(t) 2 ) = (2n + l)/2mu Q ). Thus, 

2n+l 



C(r) 



7T/2 cos(c;t) 



7 
2wi 



sin(o;ir) 



(6.53) 
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6.2.(g) Taking the equilibrium population 2n + 1 = coth(u;o/2A;.BT), 
show that this yields 



6.2.(h) By considering the fluctuation-dissipation theorem, explain why 
in the classical limit ujq 3> ksT, an exact regression theorem holds, i.e. 
in this limit, the equations governing two-time and single-time correlation 
functions are exactly related by the quantum regression theorem, even with- 
out Markovian or weak coupling approximations. 




(6.54) 



Lecture 7 



Two-level atom in a cavity: 



In the previous few lectures we have discussed a single two-level atom cou- 
pled either to a continuum of radiation modes, or considering a cavity, 
which picks out a particular mode. In this lecture we will consider more 
carefully the situations in which a cavity will pick out a single mode, by 
describing how the system changes as one goes from no cavity, via a bad 
cavity (which modifies the field strength but does not pick out a single 
mode) to a good cavity, which can pick out a single mode. 

To study this problem throughout this crossover, we will consider a two 
level system coupled both to a cavity pseudo-mode (which itself decays), 
and also coupled to non-cavity modes, providing incoherent decay. We will 
show how in this model, a bad cavity can lead to either enhanced decay, 
and a good cavity can lead to periodic exchange of energy between the 
cavity and the two-level system, i.e. the Rabi oscillations we previously 
discussed for a perfect cavity. We will then consider some examples of 
currently studied experimental cavity QED systems, discussing the values 
of the relevant parameters describing coupling and decay. 

Before discussing the crossover between bad and good cavities, we first 
investigate a toy model of a ID cavity, to put into context the meaning of 
the pseudo-mode description of the cavity mode. 

7.1 The Purcell effect in a ID model cavity 

The aim of this section is to describe how a cavity modifies the rate of 
decay of a two-level system, coupled to one-dimensional radiation modes. 
Following the previous discussion of system-reservoir coupling, the decay 
is characterised by the combination of reservoirs density of states and its 
coupling to the atom, given by: 



k 

Recalling that gu oc the field strength associated with a single photon, 
we can find the effect of the cavity on T(lj) by determining the spatial 




QED 




(7.1) 
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profiles of the modes of the full system, and inserting these into the above 
sum. 



L/2 



a/2 



Figure 7.1: Toy model of cavity, length a, in quantisation volume 
of length L. 



We will consider a cavity of size o, embedded in a quantisation volume 
of length L, where we will take L — ► oo later on. (See figure 7.1 ). Our aim 



is to find the spatial profile of the electric field modes, and thus to calculate 
T(u) oc ^2S{lo -ck)\£ k (x = 0)| 2 . (7.2) 

k 

Both in the cavity, and outside, the modes will be appropriate combina- 
tions of plane waves. Restricting to symmetric solutions (as antisymmetric 
solutions will vanish at x = 0), and matching the boundary at x = ±L/2, 
we may write: 



£k{x) 




\x\ < 



\x\ > 



(7.3) 



A k cos(kx) 

-B fc sin [k (\x\ - |)] m ^ 2 

As L — > oo, the normalisation condition will approach the simple result 
B k = y/2/L, since the normalisation integral will be dominated by the parts 
outside the cavity. We now need to find the effect of the imperfect barriers 
on this mode in order to relate A k ,B k , and thus find the quantisation 
condition specifying k, and the mode amplitude, £ k (x = 0) = A k . 

As a simple model, let us consider a varying dielectric index, e = Eq[\ + 
T)5{\x\ — a/2)], giving the equation: 



dx 2 



£ k {x) 



£k 



Let us focus on x > 0; then, given the form of £k(x) = 0(a/2 - 
0(x — a/2)g(x) with f" = —k 2 f and g" = —k 2 g, we may write: 

d 2 



dx 2 



£k 



-k 2 £ k + 25 



[g'(x) - f'(x)] + 5' (x 



(7.4) 
x)f(x) + 

~ f(x)} 
(7.5) 



Thus, to solve Eq. (|7.4|) we require that: 

f 



f 



f 



+ 9 



(7.6) 



'a\ /a\ ,/a\ ,/a\ r]k 2 

^) =9 {2)> 9 (2)- f {2) = ~ 

Substituting the forms of f(x),g(x) into this equation, one may eliminate 
A k , B k to give the eigenvalue condition: 
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For finite L, this provides a restriction on the permissible values of k. As 
L —* oo, the permissible values of k become more dense, such that, in this 
limit, one may instead consider this equation as depending on two separate 
values kL and ka; if L 3> a, it is possible to significantly change kL without 
modifying ka. Thus, we will instead rewrite this equation as a condition 
that specifies the value of kL given a fixed ka: 



cos (t) 



kr/ 

1 H sin 

2 



ka\ 

- jcos 



~2) 



sin 



2 ) 



COS 



A: a 
Y 



(7-8) 

Using this condition, one may then eliminate kL, and write A^ in terms of 
Bk = \/2/L and functions of ka, i.e.: 

, 1 (7.9) 

L y/1 + kr)sm(ka/2)cas(ka/2) + {k-q/2) 2 cos 2 {ka/ 2) 

If we define A = krj/2, then we may use this formula to give the effective 
decay rate: 



rffc 



10 







1 + Asin(fea) + ^A 2 (l + cos(fca)) 

r 



l + #) +Aa/1 + t cos(fca-flo) 



A 2 



tan 0q 



A' 



(7-10) 



u 

3 




C0/O3 n 



Figure 7.2: Effective decay rate vs frequency for the toy model of 



a ID imperfect cavity. Inset compares the exact result, Eq. (7.10) 



to the Lorentzian approximation of Eq. (7.131 



The form of this function in general is shown in Fig. 7.2 For small k, the 
effect of the cavity is weak (i.e. kr] <C 1), and so there is little modification 
compared to the result without the cavity. For larger k, there are sharp 
peaks, which can be described by an approximately Lorentzian form (see 



the inset of Fig. 7.2). This Lorentzian form can be found by expanding 



Eq. (7.10) near its peaks, which are at: 

k a = 6 + (2n + l)vr. 



(7.11) 
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We can then expand T(uj) near ojq = cko- If A 3> 1, and if we may neglect 
its k dependence, then we find: 



- -i ^ • (7-12) 

Let us introduce the full width half maximum, k, such that this expression 
becomes: 

L = PM^l f713) 

which leads to the definitions: 

k 2c „ 2 4c 



2 AV 



P = A Z = — . (7.14) 



Here To is the result for a ID system in the absence of a cavity. The 
factor P describes the maximum value of T/Tq at resonance. At resonance, 
the presence of the cavity enhances the decay rate, and P is the Pureed 
enhancement factor. It is helpful for comparison to later results to rewrite 
this as: 

J»_ 4 !5JL _£<,(»), Q^Jl (7.15) 

where Q is the quality factor of the cavity mode. The enhancement on 
resonance therefore depends both on the quality factor, and on the mode 
volume. 

Away from resonance, the decay rate decreases instead of increasing, 
since there are no modes with significant weight inside the cavity; the 
minimum value of T/Tq occurs when k^a = 9q + 2nir, and is given by 
T/Tq = 1/A 2 = 1/P, hence for this ID case, P describes both the enhance- 
ment on resonance, and the reduction off resonance. 

As r] increases further, the Lorentzian density of states becomes sharply 
peaked, and so the Markovian approximation for decay of the atom no 
longer holds, instead one can use: 

d /•*,_,,, f duo r P(K/2) 2 e^-^( t '^) 



-a~ = - J dt'<j-{t') 



dt J w J 2ir (to - lo ) 2 + (k/2) 2 

= _iv^ f^ ym ^ (716) 

This form describes the behaviour one would see if the two-level system 
were coupled to another degree of freedom, with frequency oo and decay 
rate k/2. Hence, we are led in this limit to introduce the cavity pseudo- 
mode as an extra dynamical degree of freedom. The density of states implies 
that the cavity mode cannot be treated as a Markovian system, it has a 
non-negligible memory time 1/k, and so we must consider its dynamics. 
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7.2 Weak to strong coupling via density matrices 

This section considers the case in which a single cavity mode must be 
treated beyond a Markovian approximation, by considering the Jaynes- 
Cummings model, along with relaxation of both the two-level system and 
the cavity mode. The model is thus given by: 

-p = —i [H, p] — — (a^ap — 2apa) + pa) a 



dt r 2 

r 

- — (a + a~p-2a-pa + + pa + a-) (7.17) 
H= | (a + a + a~a^ +ea z +uj Q a)a. (7.18) 

One should note that in this equation, a) is the creation operator for a 
pseudo-mode of the system, i.e. it does not describe a true eigenmode; the 
true eigenmodes are instead superpositions of modes inside and outside the 
cavity, just as in the previous section. Because the pseudo-mode overlaps 
with a range of true eigenmodes, the probability of remaining in the pseudo- 



mode will decay (see question 7.1 for the ID case). This coupling between 
cavity modes and modes outside the cavity can be described in a Markovian 
approximation, leading to the decay rate k in the above density matrix 
equation. 

In addition to the decay of the pseudo-mode, we also include a rate V 
describing decay of the two level system into modes other than the cavity 
mode. While in the one dimensional description considered previously, no 
such other channel exists, in three dimensions, if the cavity is not spherical, 
then decay into non-cavity directions is possible. In this case, one may 
consider V = TO,/4tt, depending on solid angle. For the experimental 
systems discussed below, this is generally the case, due to the need to 
have access to insert atoms into the cavity. More generally T' describes 
the possibility of relaxation into any mode other than the cavity mode; in 
solid state contexts, other possible reservoirs often exist (e.g phonons, other 
quasiparticle excitations etc). 



Let us solve the equation of motion, Eq. (7.17), starting in an initially 
excited state | 1,0). From this initial state, the only possible other sub- 
sequent states are | j,l) and | |,0). As the last of these states cannot 
evolve into anything else it may be ignored, and we can write a closed set 
of equations for three elements of the density matrix: 

Pa = Pto,to, Pb = Cab = Pto,|i ( 7 - 19 ) 

By taking appropriate matrix elements, one finds: 

jP A = -i 9 -{C AB - C AB ) - T'P A (7.20) 
jPb = ~i\{C AB - C* AB ) - K P B (7.21) 
jCab = -%(e - co )C AB + i 9 -(P A - P B ) - ^~C AB (7.22) 
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One may further simplify these equations by noting that PaPb~ \Cab\ 2 = 
is conserved, and thus writing Pa = \a\ 2 ,Ps = \P\ 2 ,Cab = ct/3* . Substi- 
tuting this leads to the simpler equations: 

dt ( P ) = ( - J/2 zA -1/2 ) ( /? ) (7 ' 23) 

where we have written A = ujo — e. In general this leads to two frequencies 
for decaying oscillations of the excitation probability: 

2 



(<"-y) (W*A-|) + ^ = 0. (7.24) 



Bad cavity limit — Purcell effect in 3D 

In the limit of a bad cavity, i.e. k <C g, T, the two frequencies corresponding 
to the above equation correspond separately to decay of the excited two- 
level system and decay of photons. Because these are on very different 
timescales, one may treat the coupling perturbatively writing v = —iT'/2 + 
S, where at 0{5) the determinant equation becomes: 

iS (iA - \ + Q + ^ = 0, (7.25) 
and since k 3> T this bad cavity limit gives: 

This describes a cavity-enhanced decay rate, T e g = r'-|-r cav , where the cav- 
ity decay rate describes the Purcell effect, as discussed above. On resonance 
we have: 

9 
2 



2 1 

K 



r — 4 

1 cav — ^ 

and recalling the earlier 3D results: 

_ u)l\d ab \ 2 g 2 _ too\d ab \' 



(7.27) 



2e V 



(7.28) 



(where we have assumed resonance, so u>o = e), one may write: 

E^-p-4^I- J-Q^ (729) 

This closely resembles the one dimensional result before; the Purcell en- 
hancement depends both on the quality factor, and also on the mode vol- 
ume. An estimate of the minimum decay rate off resonance in this case 

yields r caVimin = r caVjmax = (^j. However, if the back- 

ground decay rate T' also exists, this off- resonant reduction is not as rele- 
vant as the on-resonant enhancement. 
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Strong coupling — Rabi oscillations 

If the cavity is sufficiently good, then rather than decay, the excitation 
probability oscillates. In the resonant case of A = 0, one can write the 
general solution for the oscillation frequency as: 




(7.30) 



In order that Rabi oscillations exist, it is thus necessary that g > (T — k)/2. 
In order that they should be visible before decay has significantly reduced 
the amplitude, it is necessary that one has g > k, T. In this limit, one has 
strong coupling, and decay is indeed very strongly modified. 



7.3 Examples of Cavity QED systems 



Having discussed the role of the parameters g,K,T, and their relation to 
whether a system is weak or strong coupling, we now consider a variety of 
experimentally studied cavity QED systems, discussing the values of these 
parameters, and the advantages/disadvantages of the different systems. A 



summary of the characteristic values is given in table 7.1 the discussion 
below highlights the origins of some of these parameters. 



System 


Atom[U Q2] 


Atorngj 


SC qubit[3 


Exciton|5 




Optical 


Microwave 


Microwave 


Optical 


k/2tt 


1 MHz 


1 kHz a 


30 MHz 


30 GHz 


r/2vr 


3 MHz 


30 Hz 


3 MHz b 


0.1 GHz 


g/27T 


10 MHz 


50 kHz 


100 MHz 


100 GHz 


u/2tt 


350 THz 


50 GHz 


10 GHz 


400 THz 


X 3 /V 


IO" 5 


io- 1 


A/a = 1 


io- 1 


Q 


10 8 


10 8 


10 4 


10 4 


Other i max 


^flight 


~ 100/xs 


1/T 2 ~ 3MHz 


1/T 2 ~ 3GHz 



"Recent work reaches 10Hz|19j 



Not well known, dominated by cavity induced decay 



Table 7.1: Characteristic energy scales for different cavity QED 
realisations, as discussed in the text. 



Optical transitions of atoms 

We first consider real atoms, and cavities designed to be resonant for optical 
transitions of these atoms. In this case, the background atomic decay rate V 
is determined by the intrinsic properties of the atomic transition, reduced 
by the geometry of the cavity, and so this sets the scale that must be 
overcome for strong coupling. High quality mirrors at optical wavelengths 
can be made using dielectric Bragg mirrors. These consist of alternating 
layers of materials with different dielectric constant, with spatial period 
A/4, so that even if the dielectric contrast between adjacent layers is not 
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sufficient to cause strong reflection, the interference of multiple reflections 
will lead to strong overall reflection. 

To allow access for atoms to enter the cavity, the cavity length must 
be much larger than A ~ 0.8/im, and a typical sizepQ EE] is a cavity length 
L ~ 200/im, and beam waist w ~ 20/mi, which gives a ratio X 3 /V ~ 
1CP 5 . This small ratio means that a very high Q factor is needed to reach 
strong coupling, and Q ~ 10 8 is possible in current experiments. Typical 



values of the corresponding g,T,n are given in Table 7.1 In addition to 



the timescales included within Eq. (7.17), in some experimental systems, 
there is another timescale, that of the time for an atom to leave the cavity, 
however this is typically larger than all other timescales. 



Microwave transitions of atoms 

Remaining with atoms as the "matter" part of the Jaynes-Cummings Hamil- 
tonian, some of the problematic features above can be removed if one con- 
sidered instead microwave frequency transitions of the atoms. This in- 
creases the wavelength, allowing cavities comparable to the wavelength, 
and can also significantly improve the cavity quality. For microwave fre- 
quencies, superconducting cavities can be used (since the frequency is less 
than the BCS gap, Abcs ~THz, so the superconducting mirrors are dia- 
magnetic at these frequencies. The remaining limitation on cavity quality 
instead comes from scattering off the mirror surface, and any gaps required 
in the cavity for atom injection. For niobium cavities operating at 50GHz 
(A ~ 6mm), values of X 3 /V ~ 0.1 and Q ~ 10 10 are possible (19]. 

The relevant atomic transitions at these frequencies are transitions be- 
tween highly excited atomic states (Rydberg states). For states with n, I 3> 
1, the outermost electron remains far from the nucleus and so sees only the 
screened charge of +1, so the energy levels of such orbits are almost Hy- 
drogenic. For large n but I ~ 0, the orbit is perturbed by the enhanced 
potential near the nucleus, which can be incorporated by a quantum defect, 
Si to write: 

Ky 

(n - 5i 



with 5i —* for large I. It is thus clear that transitions between n, n — 1 
will have small energies, AE n tx 1/n 3 , and so for n ~ 50 one can have 
microwave frequency transitions. Atoms can be prepared in such states by 
using tuned laser pulses. 

As well as allowing better cavities at these frequencies, Rydberg atoms 
also have significantly reduced spontaneous decay rates without sacrificing 
coupling to the cavity mode. The total decay rate from atomic state nl is 
given by: 

r to t oc ^^iiy V \d n iy V \ 2 . (7.32) 

n'V 

Selection rules imply that V = I ± 1, which means that two different limits 
of decay rates exist, depending on whether / = n — 1 or / ^ n, as illustrated 
in Fig.[7T3 
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Figure 7.3: Cartoon of atomic level scheme, and transitions from 
a highly excited state with (a) I <C n, or (b) I = n — 1. Note that 
the actual values used are of the order n, I ~ 50, much larger than 
shown here. 



I <C n In this case, decay to almost all other values n' < n is possible (see 



arrows (a) in Fig. 7.3). We may consider the characteristic decay rate 



to nearby levels n! ~ n, or to the ground state n 1 = 1. 

• Decay to n' ~ n has u oc 1/n 3 , and the dipole matrix element 
depends on the characteristic size of the orbit for large n, i.e. 
d oc n 2 (this holds because both initial and final states have 
similar extensions). These combine to give T n i~ n oc n~ 5 . 

• Alternatively, for transitions to the ground state, oj is n inde- 
pendent, being more or less fixed at the the Rydberg energy. In 
this case the dipole matrix element is however much smaller, as 
it now involves the overlap between a large Rydberg state and 
a much smaller ground state orbital. Thus, the dipole matrix 
element depends on the overlap d oc giving T n /^i oc n~ 3 . 

As such, the latter process dominates for large n, and so the overall 
decay rate is a factor n~ 3 slower than for low excitation states. 

7i — 1 In this case, the only transition allowed by selection rules is to 
n' = n — 1. Thus, this has the decay rate in the first part of the 
above case, T oc n~ 5 and is yet even slower than for I <C n. These 
states with I = n — 1 are known as circular Rydberg states, as they 
correspond closely to classical circular atomic orbits, as expected from 
the correspondence principle. 



While the decay rates are significantly reduced by the reduced energy of 
nearby transitions, and reduced overlap of remote transitions, the coupling 
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to confined radiation for nearby transitions is not reduced in the same 
way. This is because the continuum density of states does not enter the 
calculation of g. Instead, one has g oc y/cud for resonant transitions, giving 
g oc y/n, which increases with n. However, compared to optical transitions, 
the value of g is significantly reduced because g oc 1/ y/V, and the cavity 



has V — A 3 . Thus, the value of g in table 7.1 is reduced compared to 
optical transitions, but is reduced by a smaller factor than the reduction of 
the atomic decay rate. 



Superconducting qubits in microwave resonators 

Continuing with microwave frequencies, one can also consider "artificial 
atoms" coupled to microwave cavities; these have the advantage of con- 
siderably increasing the coupling strength, by (in this case) increasing the 
number of electrons involved in the artificial atom. 




Figure 7.4: Schematic diagram of superconducting qubit capac- 
itively coupled to a stripline microwave resonator (left), and the 
equivalent circuit (right) , coloured for comparison of diagrams — 
after Koch et al. [20] . 



Figure |7.4| illustrates one schematic design of a superconducting island 
capacitively coupled to a stripline resonator. This particular design is 
known as the "transmon" qubit [20J, as the atomic states corresponds to 
quantised modes plasmon oscillations, modified by the transmission line. 
The equivalent circuit shows that the artificial atom consists of a pair of 
Josephson junctions, shunted by a large capacitance. The pair of junctions 
exists so that one may tune the effective Josephson coupling via a magnetic 
flux, giving Ej teS = Ej cos($ B /$ ), giving: 

g2 (fiQ Tobias )^ 

-ffatom = ^ — Ej jeS COs((p Q ), (7.33) 

ZL/ Q,eS 

where the number and phase operators obey canonical commutation re- 
lations. This qubit is capacitively coupled to the resonator mode, i.e. 
the charge difference across the resonator conductors voltage biases the 
qubit. This coupling can be engineered to be relatively large, of the order 
of 100MHz. In other designs of qubit, the coupling is reported to be able to 
be made yet larger, such that g > uo. This limit is referred to as ultra-strong 
coupling, in which the rotating wave approximation is invalid. 

Unlike the microwave system with Rydberg atoms, where atoms even- 
tually leave the cavity and can be measured, the observation of these super- 
conducting systems is typically via the emitted radiation. For this reason 
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K is chosen to be larger than physical constraints would require, so that 
sufficient photons escape to allow measurement. The origin of the artifi- 
cial atom decay rate in these experiments is not particularly clear, since 
coupling to bulk radiation modes should be negligible. However, since the 
system is in a solid state environment, other degrees of freedom exist that 
can lead to relaxation. These other degrees of freedom certainly lead to 
dephasing, so that in these systems, I/T2 is a significant rate, arising from 
charge and flux noise on the superconducting circuit. 



Quantum-dot excitons in semiconductor microstructures 

Finally, we consider artificial atoms at optical frequencies. One example 
of this concerns excitons in quantum dots, coupled to semiconductor mi- 
crostructures. The exciton states can be regarded as Hydrogen- like wave- 
functions of electrons and holes, with a reduced binding energy, and in- 
creased Bohr radius due to the dielectric screening £ re \ ~ 10 and reduced 
electronic mass. Although the excitonic binding energy is much less than 
optical frequency, the relevant transition is the creation of a bound exciton, 
which corresponds to the band gap less the binding energy. Typical semi- 
conductor band gaps are of the order of leV, leading to optical frequencies. 

Compared to real atoms, the notable improvement in table |7Tl~l is the sig- 
nificantly enhanced g; this arises from the increased Bohr radius, and hence 
larger dipole matrix element, as well as the much reduced mode volume. 
As the quantum dots are fixed inside the semiconductor, microstructures 
can be grown with characteristic sizes comparable to the wavelength of 
light. The mechanism of light confinement used for such cavities varies, 
combining one or more of dielectric contrast, Bragg mirrors, or photonic 
band gap materials. Dielectric contrast can produce reasonable confine- 
ment for whispering gallery modes (WGM) in circular resonators; in this 
case the WGM makes a shallow angle of incidence with the edge of the 
cavity, and so is confined by total internal reflection. Cavity decay rates 
are however typically much larger than cavities used for atoms, due to the 
lower quality of mirrors in integrated semiconductor heterostructures. The 
exciton decay rate is much larger than for real atoms due to the solid state 
environment. Just as for superconducting artificial atoms, this again leads 
both to reduced T\ and also T2 dephasing. 



7.4 Further reading 

For a discussion of the Purcell effect, and the crossover between weak and 
strong coupling, see e.g. Meystre and Sargent III [12 . The discussion 
of particular cavity QED systems in this chapter is based on various re- 
views: for Rydberg atoms, see Raimond et al. [4 or Gallagher |21j ; for 
superconducting qubits, see Blais et al. [3]; for excitons in semiconductor 
microcavities, see Khitrova et al. [5]. 
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Questions 

Question 7.1: Pseudo-mode decay rate. 

Consider a standing wave confined within the cavity, i.e. ipo(x) = 
y/2/a cos[(2n + l)irx/a]. By considering the decomposition of this wave 



onto the true eigenmodes of the extended system (as calculated in Sec. 7.1 ) 



show that the time evolution of the overlap |(V'qIV'q(*))I decays at rate re/2, 



as defined in Eq. (7.14). 



Lecture 8 



Collective effects of two-level 
atoms: open systems, 
superradiance 

So far in these lectures we have either considered behaviour of a single 
two-level atom, or if we have considered ensembles then we have assumed 
the atoms act independently. However, as first pointed out by Dicke[22 , 
this approach cannot be correct if the atoms are close together, as they 
see the same electromagnetic modes with sufficiently similar phases, and so 
emission is a collective process. This collective behaviour can significantly 
change how decay occurs, even when the electromagnetic modes are treated 
as a Markovian bath, with short memory times. 

In this lecture we will consider the case of many atoms coupled to a 
continuum of radiation modes (i.e. without a cavity); the next lecture will 
consider collective effects of many atoms with a cavity. 



8.1 Simple density matrix equation for collective 
emission 

We may begin by considering coupling between a collection of atoms and 
radiation modes: 

HsR = T (V4 ei(Wfc " e) *" ik ' ri + a+a k e i( - e -^ t+ik - r ^ . (8.1) 

k i 

Repeating the derivation of the system density matrix equation, we recover: 



d_ 
dt 



\9k 
I 2 



P (t)=-Y J j dfj: 

ij k 



1.2) 



Two extremes of this equation exist; if \n — rj\ Ao, then the atoms will 
act independently. In the other extreme, we may neglect all phase factors, 
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in which case we can introduce: 

^<Ti (8.3) 



and write: 

j t P{t) = ~\ [J + J-p ~ J-pJ + + PJ + J-} ■ (8.4) 

This second limit is what we will consider in this section. We must however 
note that when separations are small, one must also consider the effects 
of Coulomb interactions between the two-level systems. In the Coulomb 
gauge, this means the explicit Coulomb term will have an effect; in the 
dipole gauge, the problem is instead that long-wavelength modes (for which 
phase factors can not be neglected) describe the effects of Coulomb inter- 
action, leading to important energy shifts. We will discuss this effect in 



section 8.2 Putting aside these issues for the moment, let us study the 



consequences of Eq. (8.4) 



Dicke Enhancement of emission rate 

A simple understanding of how collective emission differs from independent 
emission can be found by considering eigenstates of | J\, J z . These are spin 
states and so obey: 

J 2 \J,M) = J(J + 1)\J,M) (8.5) 
J Z \J,M) = M\J,M) (8.6) 
J-\J,M) = yOVTT) - M(M- l)\J,M - 1) (8.7) 

We may then use these states as a basis for the density matrix equation, 
and writing equations for the diagonal elements Pmj = PMJ,MJ we find: 

jP m ,j = - r[j(j + 1) - m(m - i))p M ,j 

+ T[J(J + 1) - M(M + 1)]P M +V ( 8 - 8 ) 

The rate of emission is thus I = T[J(J + 1) - M(M - 1)] = T(J + M)( J - 
M + 1), whereas if incoherent, the total emission rate would be TN. 

If we consider cases where J = N/2, if J takes the maximum value 
possible consistent with the number of two-level systems, then the radiation 
rates at characteristic values of M are: 

M = f I = FN 

M =0 I (8.9) 

M =-f + l I = TN 

Thus, near M = 0, there is macroscopic enhancement of the emission rate 
compared to the incoherent case. Even when M = 1 — N/2, the emis- 
sion is much greater than it would be for a single excitation if considering 
independent atoms. 
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Solving the simple model 

The above estimates indicate the general idea of superradiance; as M de- 
creases, the collective effects increase in importance, leading to an increas- 
ing rate of transitions through the super-radiant cascade. Our aim in the 
following will be to describe this cascade. Before discussing a tractable ap- 
proximation scheme, we may note that the Laplace transformed equations 
can be straightforwardly solved to give: 

[s + r2 j]Pj,j(s) = i 

[s + T(J + M)(J-M + l)]P M ,j( 8 ) = T(J + M + 1)(J — M)P M+lyJ (s) 

The general result will thus be a sum of exponential decays with different 
time constants, and prefactors linear in t (degenerate decay rates exist, with 
M <-> 1 — M having the same rate, hence the matrix problem for eigenvalues 
is defective). While this approach may be appropriate numerically, it does 
not help to produce a qualitative understanding of the decay, so we will 
instead concentrate on approximately solving the density matrix equations. 

Semiclassical evolution after initial times 

As a first approach, we may write a semiclassical approximation for how M 
evolves. Since J does not change during the time evolution, we will drop 
the label J, and assume J = N/2 throughout. 

j t (M) = Y, M j t p M = -r Y^ M ~ ( M " + M )( J ~ M + l ) p M 

M M 

= -T((J + M)(J - M - 1)) (8.10) 

The semiclassical approximation is to assume that {(M— (M)) 2 ) <C (M) 2 so 
that expectations of products of operators can be represented as products 
of expectations, and one has: 

|(M)^-r(J 2 -(M) 2 ) (8.11) 

We have also assumed here that J, (M) 3> 1, such an approximation will 
be seen later to arise naturally whenever the semiclassical approximation 
is valid. 

These equations can be solved by substituting (M) = Jtanh(x), which 
lead to dx/dt = —TJ, hence the general solution is: 

(M) = - jtanh[rj(t - t D )}, i = -r^^ = rj 2 sech 2 [rj(t - t D )\ 

(8.12) 



Semiclassical results in Heisenberg picture The same results can 
be recovered in the Heisenberg picture by considering equations of motion 
for the operators J z ,J ± . This approach will be useful later on, when con- 
sidering an extended system. Using the commutation relations [J + , J~] = 
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Figure 8.1: Time evolution of (M) and the associated rate of 
radiation calculated semiclassically 



2J Z , [J z , J^] = ±J ± , we may write: 

j t ( ' r) = Tr ( J ~ I) = 4 Tr «- r J+ - J+J ^ = r(jv+> - 

(8-13) 

Similarly, one finds {d/ dt){J z ) = —T{J + J ). The semi classical approxima- 
tion corresponds again to factorising products of operators, thus yielding: 

j t (J z ) = -T\(J~)\\ ±(J~) = F(J*)(J-) (8.14) 

Clearly (J z ) = (M), and so the previous equations can be recovered by 
writing (J z ) = Jtanh(x) along with (J~) = Jsech(x)e' 1 ^ . The phase <p is 
constant, and is not determined by the equations of motion. While this 
substitution is most suitable to solve the equations, it is also worth noting 
(for future reference) that another substitution is more natural in order to 
understand the behaviour. By noting that the components J~ , J z represent 
parts of a collective spin, one is led to write: 

(J z ) = Jcos(#), (J~) = Jsin(0)e^, -> 6 = TJsm{6). (8.15) 

Hence, superradiance corresponds to the collective spin behaving as a damped 
pendulum, initially in its inverted state. 

Early time evolution 

In the above description, it is clear that at early times, when fully inverted, 
the classical equations will not start. However, at early times the semi clas- 
sical approximation fails. In order that semiclassics is valid, the distribution 
of M should not spread too much, and so the M dependence of the time 
evolution rate should be small compared to the mean evolution, i.e. 

2M-1|<J 2 -M 2 . (8.16) 

It is clear that this is true as long as \M\ C J, thus semiclassical evolution 
describes the late time evolution, but not the early time. At early times, 



dM 



(J + M)(J - M + 1) 
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M spreads, but at later times no significant spread occurs. This allows us 
to write: 

P M (t) = j daPj- a (to)S (M - (M(t))\ { M(t ))=j-s) • (8-17) 

i.e., if we find the density matrix at early times, then for each possible 
value of M arising from it, we may evolve forward semiclassically. The 
initial state, can be matched to the solution in Eq. (8.12), expanded for 
small s to give: 

(M(t )) = J-s~j\l- 2e 2rj (*°-' D )] (8.18) 
thus the "delay time" appearing in the semiclassical equation is given by: 

* D ( s ) = *° + ^j ln (v) ■ (8 - 19) 

This formula allows one to translate the result of the early time evolution 
at to m to the subsequent evolution of M. Small values of s, describing 
small deviations from an initially inverted state, correspond to longer delay 
times; because the inverted state is unstable the initial motion away from 
this point is exponential growth, hence the logarithmic dependence of delay 
time on the value of s. To complete the problem, we must find how the 
distribution of M spreads for early times. 

At early times, excitation numbers remain of order M ~ J, i.e. we may 
write Pj = m-s arid expand Eq. (8.8) for small s to give: 

~P S = — r(2J - s)(s + 1)P S + r(2J - a + l)sP s -i 
at 

~ -T2 J[(s + 1)P S - sP s -i] . (8.20) 

Solving the first two cases, with the initial condition Pq = 1, P s >o = 0, one 
finds: 

P = e - r2Jt , P 1 = e - r2J '(l-e- r2J *), ... (8.21) 
from which one may guess the general solution: 

Ps = e - r2Jt (l - e- V2Jt Y (8.22) 
This can be shown to satisfy the equation of motion, i.e.: 

j t P s = -T2J [P s + s(e- r2Jt - 1 + l)P s -i)] 

= -T2J[{l + s)P s -sP s - 1 }. (8.23) 

This is a Bose-Einstein distribution, P s oc z s with z = 1— exp(— T2 Jt) which 
has mean excitation (s) = z/(l — z) = e r2Jt — 1, which grows exponentially 
at early times. 

If we choose to such that 2rjio ^> 1, then (for large enough J) it is 
possible simultaneously to fulfil the condition for semiclassical evolution, 
and also to expand the above distribution for large t as giving: 

P s (t ) = exp (-r2Jt - se- r2Jto ) (8.24) 
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Putting together all of the above ingredients, one has the result: 



P M (t) = j dsexp (-T2Jt - se- T2Jt °) 



x 6 M+ Jtanh 



1 f 2J 



3.25) 



Since to was somewhat arbitrary, this result should not depend on to, and 



indeed (see question 8.2) it does not. Thus, one can set to = and find a 
full semi classical approximation of the time evolution of Pm- 



Early time evolution for Heisenberg picture The above analysis of 
Pm for short times has an equivalent manifestation in the Heisenberg pic- 
ture. For early times, one cannot factorise (J + J~) as |(J~)| 2 , but instead, 
one needs to worry about correlations of the operators. This can be treated 
approximately by assuming that the initial value of (J - ) should be drawn 
from a Gaussian distribution of mean 0, and variance (J J\J + J~\J J) = 2 J. 
These quantum fluctuations are then amplified by the semiclassical be- 
haviour, just as the semiclassical evolution of Pj amplifies small initial 
differences in M via the different classical trajectories they lead to. 



8.2 Beyond the simple model 

The above discussion assumed a contradictory set of conditions; it assumed 
that the atoms are close enough that phase differences between their cou- 
pling to light could be neglected, and that also assumed that the Coulomb 
interactions between the atoms could be neglected. As the following will 
show, one or other of these conditions can hold, but not both, and their 
violation will modify the superradiant behaviour. 



Coulomb interactions and dephasing 

If the atoms remain close together, the above treatment of coupling to 
propagating radiation modes is correct, but Coulomb interactions between 
atoms will break the indistinguishability that led to the collective enhance- 
ments. Indistinguishability is broken as each atom sees a different environ- 
ment of Coulomb interactions from its neighbours. The Coulomb interac- 
tion can be written as: 



-^coulomb = y^^ij [ a t (r j + <T t (T i ) • ( 8-26 ) 



i>j 

Treating atoms in the dipole approximation, the dipole-dipole coupling is 
given by: 

|2 



Before discussing the effect of this term, we may first notice that it de- 
pends on rather similar parameters to T, and in terms of the characteristic 
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strength, one may write: 

Mij \d a b\ 2 3vre c 3 



47re r3. \d ab \ 2 LU$ 32tt 3 



5.28) 



Thus, the same condition required in order to neglect phase differences also 
implies that Qij > T. 

Superradiance occurred because of constructive interference between 
the indistinguishable pathways of atomic de-excitation. The Coulomb term 
makes atoms distinguishable. Physically this means that while the coupling 



to light in Eq. (8.4) causes transitions between symmetric atomic states, 
Coulomb interaction make symmetric superpositions evolve into less sym- 
metric superpositions, thus reducing J. 

Coulomb interactions do not however change the excitation level M, 
it merely transfers excitation between atoms. As such, one may estimate 
the relative size of effect of the Coulomb term vs M by considering the 
number of possible states with a given M. When M = N/2, only one 
state is possible, the symmetric state with J = M = N/2, thus Coulomb 
interactions may shift the total energy, but not modify the state. For 
smaller values of M, one may use the number of representations of N spins 
with modulus J, given by: 

(2J+1)JV! 

Vn{j) = (N/2 + j + iy.(N/2-jy: (8 - 29) 

The number of possible states of given M is the number of representations 
having J > M, i.e. X]j=m u n(J)- This behaviour leads to the following 
scenario. For M ~ J, Coulomb terms do not significantly effect the decay, 
but as M approaches J (in fact, when M ~ ^/~J), Coulomb terms rapidly 
dephase the atoms, transferring the system to states with J < y/~N, and 
thus suppressing superradiance. 

If one were to work in the electric Dipole gauge, the equations may 
appear different, as no static Coulomb term exists in the Hamiltonian. 
However, in this case, the low k part of the photon mediated interaction 
describes the same physics. This means that even though a Markovian 



approximation can be made for the sum in Eq. (8.2) for nearly resonant 
modes, the low k modes must be treated carefully, and contribute an im- 
portant imaginary part (energy shift). In contrast, the different form of 
gfc in the Coulomb gauge means that the low k parts of this sum cause no 
particular problem, as the Coulomb term is already explicitly included. 



Extended systems 

To avoid the effects of Coulomb interactions, one may instead consider 
extended systems, where nj > A, and so Coulomb effects become weaker. 
However, in this case, we can no longer neglect the phase differences seen 
by different atoms, and so the evolution will be more complicated that 
Eq. djUb. 
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The simplest possible extended system is a one dimensional system 1 . In 
this case, we may replace the phase factors by exp[ik(zi — Zj)], and forget 
about summation over directions of k. Note that in this replacement, we 
have chosen to specifically consider forward propagating waves in the z 
direction. An equivalent description of backward propagating waves will 
exist, however for simplicity we neglect this complication. 

In this approximation, our equation of motion becomes: 



2 ik(zi-zj) 



>) = -£/ *E| 

ij k 

x [<rt<rje k P ~ °~ 3 P°t (& + Ck) + pvtvjtk] • (8-30) 
and we may note that: 

f dt'Ck = -n ~yT~n+ = -iP (—!—) + *8(u k - e) (8.31) 

In the previous discussion of decay, we have only considered the delta func- 
tion, and have neglected the principal value part. This can sometimes be 
valid, as the imaginary part gives an energy shift — i.e. it contributes 
a term that looks like [H e g, p], describing Hamiltonian dynamics. If this 
Lamb shift is site independent (or if we only have a single atom) then such 
a shift can be absorbed into renormalisation of the bare Hamiltonian. How- 
ever, it can also lead to effective interactions between the atoms. We have 
already mentioned one example of this; Coulomb interactions in the Dipole 
gauge originate from the low k contribution of these imaginary parts. In 
the current case, we work in the Coulomb gauge, so these low k terms are 
not important, however the imaginary parts of the resonant terms (i.e. for 
to k ~ e) will play an important role in the extended system. These terms 
lead to phase shifts that are responsible for causality emerging in the final 
result. 

To see this causality, we consider: 

El , 2 e ik(zi-Zj) 
- -A x r- (8-32) 



Concentrating on nearly resonant terms, we may approximate the density 
of states by T/2n, and extend the sum over k to ±oo, yielding: 

where kg = e/c. In the case that we may neglect phases, this becomes 
T/2, due to the step function. If we phase shift the system operators 



lr To be able to truly treat a ID system, one has two opposing constraints; the system 
must be narrow enough that there are not significant numbers of transverse modes. At 
the same time, diffraction must be small enough that propagation remains axial. These 
constraints are best satisfied by a tube of waist w ~ y/LX. 
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°i ~* °i e lk ° Zi then we find: 

Zi ^> Zj 

This describes causality, in that operators must always proceed from small 
z to larger z. Note that the previous factor of 1/2 is replaced by the 
restricted range of the sum. Because we have written a Markovian density 
matrix equation, we have neglected propagation time of electric fields, and 
so the above equation only applies for intermediate sized systems. For 
larger systems Heisenberg equations can be used to extend the range of 
validity. 

To describe the behaviour more concretely, let us consider the Heisen- 
berg equations within this Markovian approximation. Following the proce- 



dure of Sec. 8.1 we may write: 



Zi>Zj 

= -r £K + ^ +CJ X> ( 8 - 35 ) 

Zj<Z k 

j t {^ k )=2T ^KaT). (8.36) 

Zj <z k 

These equations cannot directly be factorised if they apply to single 
two-level systems, but a coarse grained version, in terms of: 

M{z) = - z^af), V{z) =J25(z- z^cr'), (8.37) 

i i 

can be considered instead. These give the semiclassical approximation 

-r dz' [V(z)V*(z') + V(z')V*(z)] (8.38) 



dN 
~dt 
dV 



dj 2T J dz'N(z)V(z') (8.39) 

If we define £ = f z dz'V(z') then this reduces to the simpler set of equa- 
tions: 

d ^L = -r[V£* + £V*), d ^ = 2TN£, d ^=V. (8.40) 
dt 1 J dt dz y ' 

Here, £ can be interpreted as proportional to the electric field strength, 
resulting from the integrated polarisation of the medium to the left of the 
current position. (Were we to have instead written Heisenberg equations 
from the begining, and were we to have avoided making the Markovian 
approximation, then a similar set of equations would have been found, but 
with the retarded time r = t — z/c appearing in place of the t.) 
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Equations (8.40) can then be understood by the same substitution as 



in Sec. 8.1, mapping them to a generaised damped pendulum problem. 
Writing M = Mq cos(9),V = Mo s'm(6)e 1 ^ one has the equations: 



d6_ 
dt' 



2T£, 



d£ 

dz 



giving the result: 



A/" sin(#)e^ 



dzdt 



2lWosin(0), 



/oc|£| 2 



dB_ 
dt 



$.42) 



Compared to the small system equation, the principal difference here is 
that this equation is now second order, i.e. it contains both damping and 
inertia. For sufficiently long samples, one thus has ringing oscillations, as 
each atom accumulates a different phase depending on where it is. This 
causes overshoot, and hence oscillations. 



8.3 Further reading 

A comprehensive review of superradiance may be found in Gross and Haroche 

Questions 

Question 8.1: Coherent states of spin 

Consider an alternative choice of coherent atomic states, known as co- 
herent states of spin: 

\^(a)) = Me aJ + | -J, J) . (8.43) 



8.1. (a) Show that normalisation of \ip(a)) implies M(l + la] 2 )^ = 1. 
8.1. (b) Evaluate (V>(«)| J z \ip{oi)) and (^(o)| J+J- \i/>{a)). 

8.1.(c) Thus, compare the rate of photon emission from a coherent state 
of spin and the Dicke states discussed earlier in this section. 

Question 8.2: to independence of superradiance P s (i) Show that 
equation (8.25) is independent of to, by evaluating the partial derivative 
with respect to to and showing it vanishes. 
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More collective effects: 
Cavity system and the Dicke 
model 

In lecture [8] collective atomic decay was discussed for coupling to a con- 
tinuum of modes. This lecture considers the case of collective atomic be- 
haviour in a cavity. Two particular cases are considered: the first section is 
on the possibility of spontaneous coherence in thermal equilibrium, and the 
second section discusses the general time evolution starting from a fully ex- 
cited atomic state. In the first section we ask whether the model in thermal 
equilibrium can undergo a phase transition in which the two-level systems 
polarise, and generate an expectation of the photon field. 

9.1 Phase transitions, spontaneous 
superradiance 

In this first section, we will consider the Dicke model without the rotating 
wave approximation. Written in terms of spin operators <7j for the two-level 
systems, this is: 

H = Y J eaf + + £ |( a t + a) {af + of) • (9.1) 

i i 

We wish to consider the thermodynamics of this model, by considering a 
mean-field theory, in which we assume a coherent mean-field for the photon 
state, and find the associated free energy. If the free energy is minimised 
by a non-zero photon field, (i.e. if such a state has a lower free energy 
than the vacuum state), then there is a phase transition to a spontaneous 
superradiant phase. For a coherent photon state exp(— |a| 2 /2 + aa^)|0), the 
associated partition function can be written as: 

Z[a\ = e-^H 2 JJ^exp (-(3 [2ga'af + eaf)) . (9.2) 

i <?% 

Here, we have written real and imaginary parts of a as a = a' + ia", and 
have used af + = 2af . One may then evaluate the sum over states of 
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the two-level system by transforming to the diagonal basis in the presence 
of field a'. The eigenvalues of the 2x2 matrices are then ±E/2, with 
E = \f e 2 + 4g 2 a' 2 , and so the free energy is given by: 



F(a) 



lo a 



N 



In 



cosh 



e 2 + 4g 2 a' 2 



(9.3) 



To determine where a non-zero value of a can minimise F, we should find 
the derivative of F with respect to a, to locate possible stationary points. 
Clearly, dF/da" = 2uia", and so for any stationary point we must have 



a 



0. For the real part of a we have instead 



da' 



4<?V 



(3 2 ^2 + 4 5 2 a /2 



tanh 



e 2 + 4g 2 a /2 



).4) 



Using the eigenvalue energies E defined above, this may be rewritten in the 
simpler form: 

dF 



da' 



2a 



Ng 2 ,(PE\ 
U -— tanhf — 1 



•5) 



A non-zero a' solution exists if the term in brackets can be made to vanish. 
In order that this vanishes at some temperature T > 0, we require that the 
function tanh(/3£' /2)/ E crosses the value uj/Ng 2 . As illustrated in Fig. 9.1 



at small temperatures, tanh(/3£'/2)/£' ~ 1/E, while at large temperatures 
it becomes 1/2T. Thus, for a crossing to occur, one requires 1/E > cu/Ng 2 , 
and since E > |e|, a transition may exist if e < Ng 2 /co. 




2 3 
Temperature 



Figure 9.1: Graphical representation of Eq 
for a solution to exist, it is necessary for l/|e 



(|9.5l. It is clear that 
> 1/E > Lu/Ng 2 . 



Now, recalling the definition of g in terms of microscopic parameters 



and matrix elements in Eq. (1.39), we have the condition: 



eui < 



V 



N (2\d ab \ 2 e 2 



e u 



(9.6) 
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Thus, apparently a transition can occur if the density of two-level systems 
is large enough. This phase transition was originally discussed by Hepp 
and Lieb |24| [25"] , and in the form presented here by Wang and Hioe [26 . 



9.2 No-go theorem: no vacuum instability 

If the density of two-level systems is large enough, the Dicke model is in- 
valid, as we have neglected the matter-radiation coupling originating from 
the A 2 terms. Let us restore these terms, and see at what density of two- 
level systems they become important. In the same single-mode approxima- 
tion in which we are treating the quantised radiation, the correction due 



these terms can be written following Eq. (1.35) as: 

(NB, the mass m appearing here for the contribution from a single atom is 
1/m = l/m e + 1/rrii, i.e. is the reduce electron-ion mass.) Including this 
term in the Hamiltonian, the function F{a) acquires an extra dependence 
on a', giving 



F(a) = uj\a\ 2 + (N4a' 2 - — In 



cosh I — \J e 2 + 4g 2 a' 2 



In terms of the condition for a phase transition to occur, this means one 
should replace u —* uj + 4£iV; i.e. the condition is now e(u> + 4£iV) < Ng 2 
Physically, this means that the A 2 terms in the Hamiltonian describe a 
dielectric response of the atoms, opposing large transverse fields; as a result, 
the energy cost of spontaneous polarisation has increased, and sponatneous 
effects can arise only if the two-level system susceptibility is larger than it 
was above. Combining the defintions of Q and g, we may write the new 
condition as: 

)<U— )■ < M > 

V \me w / V V Equj J 

As a function of density of two-level systems, this condition is clearly not 
satisfied at N/V = 0, but will be satisfied for some finite density of two-level 
systems if and only if: 

e ^<« £2 . (9.10) 

Thus, spontaneous polarisation requires that the energy and dipole matrix 
elements of an atom can satisfy 2|d a b| 2 e > q 2 /m. 



Thomas-Reiche-Kuhn sum rule 

The energy and dipole matrix elements obey a sum rule, restricting the 
relative values they may take. This sum rule can be derived from assuming 
that the internal atomic Hamiltonian depends on internal momentum like 
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Hq ~ p 2 /2m + V(r). (Such an assumption was already implicit in our 
derivation of g.) Then, as before, p/m = i[Ho,r], and so 

[r,[H ,r]] =-i-[ r ,p] = -. (9.11) 

Since this identity gives a c-number, the expectation of this identity for any 
state will match, thus: 

- = {il>\[qr,[H ,qr]]\il>)- (9-12) 
m 

Then, inserting a complete set of states \<f>), we have: 

I 
m 



qr 



J2{Wir\(f>) (<P\ [H ,qr] |V) - (VI [H ,qr] \<f> 
{d^^E^ — E^d^ — d < jy l / J (E^ — E^jd^} 
2^2\d^\\E^-E^). (9.13) 



This is the Thomas-Reiche-Kuhn sum rule, relating the total oscillator 
strength, weighted by energy, to the charge and mass of particles involved. 
There is an immediate corollary of this rule; suppose that ip is chosen to be 



the ground state. In this case, the right hand side of Eq. (9.13) is a sum of 
positive terms, and so is greater than any one of its terms. Taking the two 
lowest levels, one may then write: 

- > 2\d ab \ 2 e. (9.14) 
m 

The sum rule has therefore proved the opposite inequality to the one we 
require for a phase transition; the sum rule prevents such a phase transi- 
tion from occurring. This observation was first pointed out by Rzazewski, 
Wodkiewicz, and Zakowicz |27j . 

9.3 Radiation in a box; restoring the phase 
transition 

The above result shows that the phase transition by which the vacuum 
state becomes unstable, and spontaneously polarises is an artefact of that 
model, and not physical. This however is not the end of the story. There 
is a simple extension to the Dicke model, which is appropriate in a variety 
of recent works on confined quantum optical systems, which can restore 
the phase transition. That extension is to consider a closed system with a 
density of excitations. This means considering the at-first counterintuitive 
idea of a chemical potential for photons. 

The idea of a chemical potential for photons is not generally considered 
because photons cannot easily be confined, and so a fixed density of photons 
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S3 



is hard to achieve. Without confinement, photons can be exchanged with 
the bulk, which acts as a reservoir at zero chemical potential. However, in 
engineered cavity quantum electrodynamics systems, photon confinement is 
exactly what is being created, and so for such systems, including polaritons 
in semiconductor microcavities, Josephson junctions in microwave waveg- 
uides, and atoms in superconducting mirror cavities, a chemical potential 
for photons is a useful concept. 

Adding a chemical potential, we have 



H -> H = H - (j,M, 



M 



n + 



+ a) a 



(9.15) 



The net result is to replace e — ► e = e — fi and uj ^ Co = to — fi in the previ- 
ous Hamiltonian (although the factors of e and to in the coupling strengths 
remain unmodified). In this case, it is clear that one can satisfy the condi- 
tion e {Co + 4£iV) < Ng 2 by ensuring the chemical potential is close enough 
to the two-level system energy. Since this may be done at any density of 
two-level systems, we may assume a low density, and so neglect £ — i.e. 
neglect the A 2 terms in the Hamiltonian. 

Since a transition may occur, it is interesting to find the critical temper- 



ature of this transition. From the adapted form of Eq. (9.5), this condition 
is: 



(e — n) (u — n) = Ng 2 tanh 



li 



2T 



(9.16) 



This can be combined with the equilibrium expectation of the number of 
excitations (M) in the absence of an expectation of a, which gives. 



(M) 
N 



1 — tanh 



2T 



(9.17) 



Combining these allows one to plot the phase boundary, as shown in 



Fig. 9.2 See question 9.1 for further discussion 



9.4 Dynamic superradiance 

In this section, we consider the general time evolution of the Dicke model, 
starting from a fully inverted atomic state. This describes the analog of the 
superradiance discussed in Sec. [8] in the case where there is a good cavity. 

For illustration, we discuss here the simplest case of e = a>, and so our 
model is: 

H = ujJ z + | (a^r +aJ+) +uoa)a. (9.18) 
To study the semi classical dynamics, we first write the Heisenberg equa- 
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0.2 0.3 0.4 
Density: <M)/N 



co/gVN=0.0 
G)/gVN=0.5 
co/gVN=1.0 
o'.5 



Figure 9.2: Critical temperature vs excitation density, plotted for 
three different values of oj (measured with respect to e) , including 
the case of ui — e discussed in Q. 9.1 The inset shows 1/T vs 
density on a logarithmic scale, illustrating the asymptotic form at 
low densities. 



tions of motion: 



-ujjy + ^i(a-a^ J z , (9.19) 

+oj J x - |(a + J z , (9.20) 

-^i(a-a^ J x + ^(a + a^ J y , (9.21) 

9 -{J x - U y ) + wo. (9.22) 

The semiclassical approximation is then to replace these equations for non- 
commuting operators for equations for their commuting expectations, and 
to factorise the expectations, so ((a + a))J z ) — > (a + a)){J z ). Doing this, 
and writing (a) = a = a' + ia" the above equations can be written as: 



r 


= %[H, J x 


j y 


= i[H, Jy 


j z 


= %[H, J z 


ia 


= -[H,a 



gol 



h x J, h = ( -ga" J , ia = | ( J x - iJ y ) + coot. (9.23) 

In order to solve these equations, it is first convenient to transform to a 
rotating frame. This means substituting a — > ae~ lujt , and J~ — > J~e~ lujt , 
i.e. 

cos(u;t) — s'm(ojt) 

sin(cjt) cos(wt) | J. (9.24) 
1 




ia = | (J x - iJ y ) . (9.25) 



2 
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At this point, starting in the excited state J = Jqz, we can make the ansatz: 





(9.26) 



Then, substituting this into the equations of motion, we have: 



-a" cos(0) 
-a' cos(#) 
a' sin(0) 





' > 






cos(0) 











(9.27) 



which implies a" = and 9 = —ga'. Then, the equation for a = a' gives: 

a' = _^ si n(0). (9.28) 

Thus, the angle defining the Bloch vector obeys the equation of motion: 

9 2 Jo 



-ga' 



sin(0). 



(9.29) 



This is the equation of motion for an (initially) inverted pendulum, and so 
the Bloch vector makes traces out great circles passing through the entirely 
inverted, and entirely empty states. One should compare this second order 
equation for dynamical superradiance in a cavity to the first order equation 



in the decaying case discussed in Sec. |8.1| of the previous lecture. 

When non-resonant, i.e. e ^ lo, the comparable equations are more 



complicated (see question 9.2), however the general result is similar: af- 
ter transformation to a rotating frame, the Bloch vector traces out circles 



(see Fig. 9.3), and can always be mapped to the problem of an inverted 
pendulum. 




Figure 9.3: Bloch sphere, and paths that Bloch vector takes in 
the rotating frame, in the general (non-resonant) case. 
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9.5 Further Reading 

The original phase transition of the Dicke model is discussed in Refs. 2 1 
[26] . The "no-go" theorem is introduced in Ref. |27j and elaborated further 
by, e.g. Ref. |28j . The phase transition of a model with a chemical potential 
is discussed by Eastham and Littlewood [29, 30J. 

As hinted to in the discussion on dynamic superradiance, the Dicke 
model is integrable, and this integrability extends even to the case of a 
distribution of two-level system energies: 

H = Yl ei(j i + wata + Yl f ( ata -r + aa t) ■ ( 9 - 3 °) 

i i 

The integrability of this model can be proven by defining a vector function 
L(z), having the following properties: 



• The modulus of the vector L(z) is conserved. 



• The Taylor expansion in z contains as many powers of z as there are 
two-level systems. 



Together, these mean that L(z) contains sufficient conserved quantities that 
the system is integrable. This method is discussed in specific case of the 
Dicke model by Yuzbashyan, Kuznetsov, and Altshuler [31] . The history of 
this integrability, and its relation to a wider class of models is discussed in 
the review article by Dukelsky, Pittel, and Sierra [32 j . The time-dependent 
solution, starting from the excited state is given by the same method used 



in question 9.2 by Barankov and Levitov |33] 



Questions 

Question 9.1: Asymptotic form of phase boundary 



9.1. (a) By combining Eq. ( |9.16| ) with Eq. (9.17), show that at low den- 



sities, in the resonant case (e = oj), one may write: 

Tc ^ ln(<M)/iV)' (9 " 31) 
9.1.(b) Find the low density asymptote in the more general case e ^ oj. 

Question 9.2: Non resonant Bloch vector precession 

Consider the non-resonant case of the Dicke model: 

H = eJ z + - (a),r + aJ + ) + uja)a. (9.32) 



9.5. FURTHER READING 



87 



9.2. (a) Show that the semiclassical equations of motion may be written 
as: 

ij- = e J~ - 2gaJ z , (9.33) 
i j z = g (j+ a - J- a *) , (9.34) 
ia = gJ~ + u>a. (9.35) 



9.2. (b) By making a change to a rotating basis, but using a — ► ae~ ir,t , 
J~ — > J~e~ lTlt , show that this changes the above equations by replacing 
ijj — > Co = u — T) and e — > e = e — r\. 

The solution to these equations, analogous to the one discussed in 
Sec. 9.4 has the special property that a is real. Making the ansatz J~ = 
Aa + iBa, J z = Jq — Ca 2 , and using 9(a) = 0, show the following: 

9.2. (c) The imaginary part of the equation for J_ implies A = IB; 

9.2. (d) The equation for J z implies C = gB/2; 

9.2. (e) The equation for a implies gB /2 = 1, and u + e = 

The last of these conditions fixes rj, i.e. r/ = (to + e)/2. This leaves only 
the real part of the equation of motion for J_ . 

9.2. (f) Show that this equation has a first integral, and that this first 
integral is equivalent to the condition that \ J\ 2 = J + J~ +{J Z ) 2 is a constant 
of the motion. 

9.2. (g) Using the initial condition \J\ 2 = Jq, show that this final equa- 
tion may be written as: 

2 



d=|ayaQ — a 2 , Uq = 2Jq — 4(-) (9.36) 



9.2. (h) Thus, writing a = oto sin(0/2), show that the equation of motion 
has the form:: 



g- 



sm{9) 



(9.37) 



9.2. (i) Find the form of J z in terms of 9, and compare to the solution 
in the resonant case. 



Lecture 10 



Lasers and micromasers 



In this lecture, we move to discuss a somewhat more classical system, the 
laser. The laser is classical in the sense that far above the lasing thresh- 
old, it produces large coherent states with classical correlations, and it is 
based on stimulated emission of radiation, which can be described classi- 
cally. However, we will discuss not only simple lasers, for which reasonable 
classical descriptions may exist, but also more quantum mechanical sys- 
tems, such as the micromaser, in which two-level atoms fall, one at a time, 
through a cavity, and in which the quantum mechanical evolution of each 
atom matters. 

10.1 Density matrix equations for a micromaser 
and a laser 

The aim of this section will be to write down a density matrix equation 
describing a laser, and to find the steady state of this equation. The laser 
will consist of a cavity mode, coupled both to the continuum of modes 
outside the cavity, and to a gain medium that provides photons to the 
cavity mode. The modes outside the cavity lead to the standard decay 
term: 



For the gain medium, we will begin our discussion by discussing a formalism 
that can apply both to few atom lasing systems, and also to the typical 
incoherent gain medium, which will be introduced in the next section. 

Gain from two-level atoms 

To give a model of gain that encompasses both coherent dynamics of a 
single injected atom (for a micromaser), as well as describing incoherent 
evolution of an externally pumped atom (for a laser) , we will describe the 




which can be written in the Fock basis as: 



-rPnn< = (n + n')p nn , -2^J{n + l)(n' + l)p n +i,n'+i ■ (10.2) 

at decay 1 L J 
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Figure 10.1: Left: Micromaser scheme, atoms in their excited 
state are injected at a rate r, and fall through the cavity in time r. 
Right: Relation between a many- level lasing scheme, with decay 
via intermediate states to the ground state, and pumping to the 
excited state, and the quasi two-level description, with injection 
rate r, and two-level system lifetimes distributed according to the 
decoherence rate 7. 



dynamics of the gain medium, allowing for arbitrary rates of decoherence 
and pumping. We therefore consider a two-level system, which starts in 
the excited state, and interacts with the cavity for time r; new two-level 
systems are injected at a rate r. If we take r to be the same for all atoms, 
we expect to see periodic dependence on t, as one allows a given number 
of Rabi oscillations, transferring energy between the two-level system and 
the cavity. This model is appropriate for atoms falling through a cavity at 
a fixed speed. If instead r is drawn from a Poisson distribution, average 
rate 7, then one recovers the standard laser result; at a rate 7 the two-level 
system decays to some other states, and then rapidly decays to the ground 
state, and is then re-excited to the excited state of the two-level manifold 
at a rate r: the net result is loss of coherence at a rate 7, and injection of 



energy at a rate r. This relation is illustrated in Fig. 10.1 



To find the gain supplied by the two-level systems, we want to find 
the evolution of the field density matrix. If interaction between the two- 
level system and the cavity for time r causes the change Apg e id(r), and 
new excited two-level systems are created at rate r, then one can write the 
evolution of the field for either a fixed r, or a Poissonian distribution as 



dt 



Afield 



dt 



PReld 



/•oo 

■ / 7 (irA /0field (T)e-^. (10.3) 
J 



Thus, we require the change of the field state due to evolution for time r: 



ApfieldO) = Tr atom [p(t + t)] - /9fi e ld(*o) 
= Pfield.U (*0 + T) + Pficld,TT(*0 + 



•pfield,TT (*())• (10.4) 



In the second line we have used the fact that the two-level system is by 
definition excited at t = to- Now, using results of Sec. 3.2 we can follow 
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the evolution of a state |f, n) with n photons to time r: 



, n) 



cos i — - — r I + % cos(Zti) sin I — - — r 
+ isin(20)sinf%tir^ \i,n-l), 



IT,™) 



(10.5) 



with fl n = \/(e — to) 2 + g 2 n, tan(20 n ) = g^fnjie — u>) as in chapter [i] For 
simplicity, let us assume resonance, e = u>, so 9 = tt/A. Then, using this 
wavefunction evolution, one finds: 



Ap„„'(r) = p nn i cos ^|rv / n+ l) cos ^|r\/n' + 1^ - 1 

+ /3 n -l,n'-l sin (l 1 "^) sin (| tvV^ . (10.6) 

It is apparent this equation describes gain, since it describes transfer of 
probability from the state with n — 1 photons to the state with n photons. 



In distinction Eq. (10.2) describes evolution from n + 1 photons to n pho- 



tons. Note also that as anticipated, the fixed time r would lead to possible 
cancellation if integer numbers of Rabi oscillations occur. However, since 
the oscillation period is different for each number of photons, the general 
result is very complicated, and can lead to a very non-smooth probability 
distribution. 



Gain medium with decoherence — laser limit 

Now let us consider the case where r is a random variable, distributed ac- 
cording to exponential decay, rate 7. Thus, we should average Eq. (10.6) 
over such a distribution. For simplicity, we will present here only the cal- 
culations for the diagonal elements, p n = p nn ; for the general case, see 



Q. 10.1 



dp n 

dt 

where: 



gam 



sin 



2 (9 



rVn+T) ) p n + (sin 2 ( f r 



n ) Pn-l 



110.7) 



a 

sin" I — r 
2 



"fdre 



1 



a 

sin" I — r 
.2 

1 



. + 

7 + la 7 — la 



dre~ 1T (e iar + e~ iar - 2) 
a 2 /2 7 2 



1 + a 2 / 7 2 



(10.? 



Thus, the rate of gain due to a such a distribution of atomic coherence 
times is given by: 



dp n 
dt 



gam 



(ffV7 2 ) 

1 + (<?2 n / 7 2) 



Pn-l 



(g 2 (n + l)/ 7 2 ) 
l + ( 9 2( n + i)/ 7 2) 



Pn 



(10.9) 
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10.2 Laser rate equations 



Combining the gain in Eq. (|10.9|) with the diagonal part of the decay in 

Lotion for the 

( 5 2 (n + l)/ 7 2 ) 



Eq. ( 10.2 ), one finds the equation of motion for the probability distribution: 

(g 2 n/j 2 ) 



dp n 

dt 



1 + ( 5 2 n / 7 2) 



Pn-1 



l + ( 5 2( n + l)/ 7 2)^ 

- K [np n - (n + l)Pn+l] 



(10.10) 



From this equation we can now determine when lasing will occur, as well 
as the full probability distribution when either well above or well below 
threshold. 



Laser threshold condition 



Lasing occurs when the vacuum state becomes unstable, i.e. when dt(n) > 0. 
Thus, let us find: 

d t (n) = Y t nd tPn . (10.11) 

n 

Let us however assume that g 2 n/j 2 is small — i.e. that the average pop- 
ulation is small compared to (7/g) 2 , so that the denominators in the gain 
can be expanded. Then, 



t (n) = y^re < 

n ^ 



9 

9 n 



g 2 n\ g 2 (n + l) 

1 n" Pn-l 2 

T J T 



1 



g 2 (n + l) 



Pn 



- K[np n - (n + l)p n+1 ] 
Now, shifting the various sums, one can rewrite this as 



'10.12) 



d t (n) = 



Pn 



[(n+l) 2 -n(n + l)] 



(^) [(n + l) 3 -n(n+l) 2 ] 



[n 2 -n(n-l)] \, (10.13) 



which can be written as: 

-.2 



d t (n) 



rg^ 
2 7 2 



K 



, \,rg r 
{n) + 2Y 2 ~2 



k \ n 



^) 2 ((n + l) 2 ) 



(10.14) 



This can be understood as three terms; the first describes competition 
between loss and stimulated emission; the second describes spontaneous 
emission; the third describes the nonlinear susceptibility of the two-level 
systems, reducing the gain they provide when (n) is large. The first and 
third of these could be recovered from a semiclassical theory of lasing, 
however the existence of spontaneous emission requires a fully quantum 
theory. However, if our aim is the lasing condition, this can easily be 
identified from the first term: stimulated emission outstrips decay if g 2 > 
2«;(7 2 /r), and lasing then occurs. 
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Steady state probabilities 



The coefficients in Eq. (10.10) depend on n in such a way that one may 
write: 

= -B n p n + B n „ip n _i - C n p n + C n+1 p n+1 , (10.15) 

with excitation rate B n = rg 2 (n+l) /[y 2 +g 2 (n+l)] and decay rate C n = nn. 
To find a steady state that is valid for all n, we need to find a relation 
between terms on the right hand side such that this will always vanish. 
The only consistent equality is to choose: 

B n pn = C n+1 p n+1 <S=^ B n ^xVn-\ = C n p n . (10.16) 

The steady state distribution is then given by: 

nB m -i tt r(g 2 /-f 2 ) , 1(117 , 

^Po=U 2mh 2 )} Po, (10.17) 

m=l m=l L K iiii 

with po set by normalisation. This product can in general be expressed 
in terms of a hypergeometric function. Let us consider two limits where a 
simpler expressions result. 



Below threshold If we are at small powers, then small values of n will 
be most probable. In this case, one may neglect the term involving g 2 raj^ 2 
in the denominator. Thus, at low pumping rates, below threshold, one has 
p n oc (rg 2 /2^/ 2 k) u , which is a thermal distribution, with exp(— (3E e g) = 
(rg 2 /2j 2 k), from which one recovers 

1 r(q 2 h 2 ) 
<"> = ^K^l = 2 JVwW <10 ' 18) 



in agreement with the small pumping rate limit of Eq. (10.14). Clearly, 
such an approximation is only valid below threshold, otherwise p n cannot 
be normalised. In this low pump limit, there is an "effective pump rate", 
limited by the ratio of 5/7, indicating that emission competes with relax- 
ation. 



Far above threshold The opposite limit occurs when (n) is large enough 
that the term 1 in the denominator can be neglected. In this case one has 
p n oc (r /2k) u jn\. This is a Poissonian distribution, as one would expect for 
a coherent state, with average population r/2/t, set only by the balance of 
pump rate and decay rate, since strong fields mean that stimulated emission 
beats any other relaxation process, thus for every excited atom created, half 
a photon is emitted. 



10.3 Laser Linewidth 

Our analysis of the laser so far was based on the density matrix equation and 
its steady states. Just as in lecture[5j we can also calculate the spectrum, by 
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using two-time correlation functions via the quantum regression theorem. 
In the current context, we have a density matrix equation for the cavity 
field, and so we can proceed directly to evaluate: 

I{u) oc 25ft ( dte iut (a^{t)a{0)}\ (10.19) 



As in the case of resonance fluorescence, we can thus find this by determin- 
ing the expectation of (a'(t)}, when the initial state at t = is taken to 
result from acting with a on the equilibrium density matrix. 

If we have the time-dependent solution of the equation of motion for 
the density matrix, then we can write 

(ot(t)) = Tr [aV(t)] = \^+Tp n ,n+l(t). (10.20) 



We will define the shorthand p n ,n+k 

(t) = pfc(t), and so our task is to find 
Pn(t). Before writing its equation of motion, let us anticipate how it should 
behave. Regardless of the initial conditions, at long times this expectation 
should decay to zero, as it does not conserve the number of photons. How- 
ever, since the initial conditions are those set by ap e qbmj the expectation 
of (at(t)) will be non-zero; we are interested in how this function decays in 
time. 



Equation of motion 



Combining the off-diagonal gain found in Eq. (10.31) with the decay term 



from Eq. (10.2) gives an equation of motion for the off-diagonal density 



matrix that can be written in the form: 

( nPn ' n — lPn — 1 ' ti-\-XPti-\-X' 

Writing r\ = g 2 /j 2 for brevity, the coefficients are given by: 



dt 



(10.21) 



rjy/(n+ l)(ra + 1 + fc) 
1 + rj(n + 1 + fc/2) + r] 2 k 2 /16 



C\ = Ky/n(n + k), 

r)(n+ 1 + k/2) +r] 2 k 2 /8 



A 



1 + rj(n + 1 + k/2) + r] 2 k 2 /16 



+ k \ n + 



(see question 10.1 for derivation of the off-diagonal gain terms). Note that 

they 
0, the 



the terms B% , C„ have been written in such a way that for k ■■ 
reduce to the previous B n , C n discussed above. Whereas for k 



term A\ is simply related to the by A^ = + C°, no such simple 

relation exists in general. Since it was this simple relation that lead to the 
detailed balance condition for the probability distribution, no equivalent 
condition exists with k ^ 0. 
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Approximate solution 

Although the previous detailed balance solution does not apply to k ^ 0, 
it motivates a possible extension, which gives a solution decaying in time: 



P k n(t) 



n j^k 

II ~pk 1 P0> 
m=l m 



(10.22) 



where is a parameter to be determined self consistently. Substituting 
this ansatz for p k into Eq. (10.21) gives: 



- D k p k 



~^nPn ' D n-\ 



r>fc 

D n-\ 



+ Q 



k 

n+1 



r>k 1 n 



(10.23) 



The terms in brackets are the expressions for p\±\ in terms of p k from 
Eq. (10.22). If D k ~ -D^-tu then the exponential terms may be dropped, 



giving: 



r\k Ak /-ik rtk 

1 ] -<^T) - ' > I 1 



(10.24) 



Thus, if we solve Eq. ( 10.24 ) , and find it varies slowly with n, then D ~ 
defines the decay rate of the off-diagonal density matrix. We therefore 
want to find D 1 to solve for the linewidth. If we assume we are far above 
threshold, so (n) 3> 1, then we can expand the square roots in C\ and B^, 
to give: 



D 



rj({n) + 1 + 1/2) + r/78 - r/((n) + 1 + 1/2) 



7](n) + 1 + 7/(1 + 1/2) + 7/2/16 



+ k I (n) + - - (n) 



1 + 



2(n> 



1 

5(77)2 



+ o 



1 

(77)2 



1 



+ K 



(10.25) 



8(n) V2 7 2 

Thus, using this time-dependence in the two-time correlation, we have: 



9i(t) 



-Dt 



(10.26) 



We have been working in an interaction picture, neglecting free evolution of 
the photon fields at frequency to. Restoring this, the frequency dependence 
of the intensity is then given by: 



D 



(v-lu) 2 + D 2 

Hence the FWHM of this Lorentzian lineshape is given by: 



2D 



1 



4(n) \2t 2 



r 5 



(10.27) 



(10.28) 
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Interpretation of linewidth 



This result for the lineshape has a simple interpretation in terms of phase 
diffusion. If we consider the radiation field of a laser to be described by an 
ensemble of coherent states, with effectively fixed magnitude |or| = \J (n), 
but different phases, then decay of correlations occurs because of growth of 
phase differences. Each time a photon is incoherently emitted emitted into, 



or lost from the cavity the phase changes. This is illustrated in Fig. 10.2 




Figure 10.2: Cartoon of coherent state undergoing a random walk 
due to spontaneous emission events, leading to phase diffusion for 
the coherent state. 



To turn this cartoon into a physical estimate, we start from the rate 
of such changes per unit time, which is d = k + rg 2 /2'y 2 , by adding rates 
of photon decay and spontaneous emission (hence the appearance of the 
effective pumping rate, as appearing in the threshold condition, rather than 
the total pump rate, as in the probability distribution far above threshold). 
If we assume a random walk, with d events per unit time, the change to the 
coherent state is governed by the random walk probability distribution: 



P(Aa) 



exp 



(Ao) 
td 



(10.29) 



However, the variation is entirely phase, so Aa = \J (n)A6. Using this, the 
decay of correlations after time r can be found by averaging the state: 

gi (t) = < qt ffl, Q i 0) > = ( e -iie+m)] e i<> 




dW-exp(-(A*)-!A_^ 



exp 



td 



, (10.30) 



which reproduces the previous result, D = d/4(n). 



10.4 Further reading 

The laser rate equations presented here follow closely the presentation in 
Scully and Zubairy [ID]. For a different presentation of laser theory, the 
review article by Haken |34j . The discussion of laser linewidth in this 
Chapter is based on Scully and Zubairy [ID] . 
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Questions 



Question 10.1: Off-diagonal matrix elements 

Evaluate the general off-diagonal matrix elements describing gain, de- 
scribed in Eq. (10.6), and thus show that, writing p n ^ n +k = Pni one nas: 



dpn 
dt 



gam 



(n + l + fc/2)( g V) + fcV/8 
7 4 + (n + 1 + k/2)g 2 ^ 2 + k 2 g 4 /16 

y/n(n + k)(g 2 j 2 ) 



+ 



7 4 + (n + k/2)g 2 ~f 2 + k 2 g A /\Q 



P k n-l 



(10.31) 



Question 10.2: Second-order correlation functions 

Using the quantum regression theorem, find an expression for the second 
order correlation function: 



92 (r) 



1 



(10.32) 



in terms of the equations of motion and initial conditions for the variables 
Pn = (n\apa)\nj . 

Solve these equations and initial conditions explicitly in the two limits 
g S> 7 and 5^7 (i.e. far above, and far below threshold), and comment 
on the form of gi (r) in the two cases. 

[For values near the threshold, one can solve these equations numerically 
to find the full evolution of the second order coherence function.! 



Lecture 11 

More on lasers 



In the previous lecture, we developed a quantum theory of the laser and 
micromaser, in terms of the density matrix of the photon field. The aim 
of this lectures is firstly to repeat that derivation, this time starting from 
a more familiar density matrix equation for the coupled system. This will 
then make clear the approximations required in the results discussed last 
time. It will also allow one to see the relation between the previous results, 
and the semiclassical (Maxwell-Bloch) equations. 

In discussing the relation of the semiclassical theory, and the results 
of the full quantum theory, we will see that there is a parameter, f3, that 
controls the ratio of emission into the cavity mode vs all other decay chan- 
nels of the excited gain medium. When /3 is small, there is a well defined 
threshold, and the semiclassical theory can describe this adequately. When 
(3 is large (which will occur for strong coupling) the threshold becomes less 
well defined. To further study the sharpness of threshold, we also consider 
the uncertainty in photon number, which peaks at threshold. In terms of 
this quantity, we can also re- analyse the micromaser, looking at its steady 
states and noise levels, seeing that there are peaks of noise near transi- 
tions, but sub-Poissonian noise between them. Armed with these concepts, 
we will then briefly discuss recent cavity QED realisations of single atom 
lasers and micromasers, noting the novel features these experiments show. 



11.1 Density matrix equation 

As in the previous lecture, we will consider an effectively three-level de- 
scription of the laser gain medium. The description effectively involves 
only three levels, as we assume recycling to the ground state is fast fol- 
lowing any transitions out of the lasing levels, thus only the ground and 
excited states of the lasing transition, g, e and the true ground state are 
important. (See Fig. 

The coherent dynamics in this system is therefore controlled by the 
simple Hamiltonian: 

#=§EH 9 +°M e )- ( 1L1 ) 

i 
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e 




y -, 



lasing 
transition 



Figure 11.1: Labelling of three level scheme 



The notation erf is the operator that takes atom i from state g to state e. 
Since we have three active levels, we need this extended notation, rather 
than the two-level system definitions of a^,a z . In addition to the Hamilto- 
nian, there are incoherent processes describing pumping, decay and dephas- 
ing. The description of pumping as purely incoherent would correspond to 
coupling the — > e transition to an inverted reservoir. Alternatively, it may 
be seen as an approximation description of coherent pumping to a fourth 
level — this is discussed further below. 

In describing pumping and decay, we assume these processes act on 
each atom in the gain medium separately. As such, we separate the rate 
of pumping of an individual atom T from the total pumping rate r that 
appeared in the previous lecture — we will see later what quantity we 
should identify with r. Also, in distinction to the last lecture, we include 
both decay of the two lasing levels, and in addition pure dephasing, with a 
rate 711. Putting these ingredients together we have: 



In writing this, we has used the fact that each term acts on a single atom 
to collapse adjacent atomic operators, e.g. a 0e a e0 = a 00 , hence the sim- 
pler form of the pumping and decay terms. In the following, we will drop 
the subscript % on the atomic degrees of freedom, assuming all atoms are 
equivalent, and so the sum can be replaced by the number of atoms Na- 

Reduction to photon density matrix 




+ I (of p - 2of pa? + P°t e + °?P ~ 2<rJVf + P°f ) 
- I [« £ " of 8 ) 2 P ~ 2 « e " of) P (cr! e ~ <?l 9 ) + P (of " < 



f) 2 ]}- 
(11.2) 



We will now reduce Eq. (8.4), which contains both field and atomic degrees 
of freedom, to an equation for just the field degrees of freedom. This means 
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we want to introduce: 

p,p = TV atoms (p) = Pee + Pgg + Poo (H-3) 

The extra notation pij indicates an operator in the photon space, but a 
c-number in the atomic degrees of freedom; i.e. we define: 

Pij = Tr atoms (a J V) (11.4) 

where n, m are photon number states. With this notation, taking a trace 
of Eq. (11.2) over atom states gives: 

■9 N A ( , + t 



dtPip = —i— ^— ( apge + a' p e g - p ge a - p eg a} 

[a^ap^ — 2ap^a) + p^a^aj . (11.5) 



K 

2 



Thus, to proceed, we must also find equations for p eg etc. Taking a trace 
over atoms with a factor of a ge then gives: 

dtPeg = -y (ap gg - p ee a) - ^ (a)ap eg - 2ap eg a) + p eg a)aj 

-l{l + l)p eg - 1 l(l + 2+l)p eg . (11.6) 

At this point, we start to make use of the same approximations as made 
in the previous lecture; i.e. k <C 7. This firstly means we can drop the 
k term in the above equation. It secondly means that, for the purpose of 
substituting into Eq. (11.5), we may find the steady state of the atomic 
degrees of freedom. Defining jt = 7 + 7|| , this gives the condition: 

= dtPeg = ~— (ap gg - p ee a) - ^tPeg, (H- 7 ) 

giving the result: 

Peg = -^(ap gg - peed). (11.8) 

27t 

We can now repeat the above logic to find find the relation of the di- 
agonal elements p ee , p gg , poo- Dropping the k terms, and looking for steady 
states, one find: 

= -y (a ] p eg ~ pgeOj ~ ~fPgg 

= -|— (atapgg + p gg a)a - 2a) p ee aj - jp gg (H.9) 



"y (aPge ~ PegO^ - 7p ee + T p 0Q 

- — — I nn^ n J- n n n t — Onn rjt 



, y aa 1 p ee + peeaa 1 -2ap gg a l j -JPee + Tpoo (11.10) 

"EpoO + 7 (Pee + Pgg) ■ P- 1 - 11 ) 
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Equations (11.9) and (11.10) have made use of the result in Eq. (11.8) 



Furthermore, Eq. (11.11), along with the definition in Eq. (11.3) implies 
that: 



Poo 



7 



r + 7 



Pee + Pgg 



r + 7 



(11.12) 



We then have a closed set of three equations for p^, p ee , p gg . However, 
because these equations involve photon operators, it is rather involved to 
invert them to get a single equation for p^. To do this inversion, it is easiest 
to work in the photon number basis. Restricting to the diagonal elements, 
we have equations for p n ,p e ,n and p g . n , which become: 



9 

,2 



{2np g)n - 2np e!n - 1 ) - 7p fl) , 



= ~j— (2(n + l)p e ,n - 2(n + l)p g . n+ i) - -yp e , n + J~?—Pr, 
47t 1+7 



(11.13) 
(11.14) 



Shifting n — > n — 1 in Eq. (11.14), it is clear one can combine the two 



equations, by taking their difference to give: 

,2 







9 



-n { 



r 7 



; (j>g,n ~ Pe,n-l) ~ 7 (j>g,n ~ 7Pe,n-l) ~ TTT Pn-1- 

It 1+7 
This equation can be solved for — p e , n -i to give: 

1 r 7 



(11.15) 



Pg,n Pe,n—l 



7 + 5 2 n/7 t T + 7 



Pn-l- 



(11.16) 



Then, noticing that the atom-photon coupling terms in Eq. (11.9) and 



Eq. (11.10) add up to give the required term in Eq. (11.5), one has the field 



equation in the form: 



dtPn 



9 2 N A 
47i 



[2n(p fl)Tl - p e ,n-l) - 2 (ra + l)(p g ,n+l -Pe 



nnpn + «(n + i)p n +i (11-17) 



and substituting Eq. (11.16) recovers the previous field equation: 



_ r ff 2 ra/7 f r g 2 {n+l)/^ t 

(JtPn — n i 9 / Pn-1 n 9/ ,w Pn 

- Knp n + n(n + l)p n +i (11-18) 

where r = Na^T /(7 + T). Clearly, if 7 T, then r ~ TNa, so T is the 
pumping strength, but if F becomes too large, then the pumping rate is 
limited by the recycling rate 7, as the true ground state becomes depleted. 
If 7$ = 7, (i.e. if 7|| vanishes), the above reduces exactly to the results of 
the previous lecture. 
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Four-level scheme 

As a brief digression, let us note how a four level scheme, as illustrated in 



Fig. 11.2 reduces to the above problem in a limit of fast relaxation. In this 
system, there is also a level 1, with a coherent coupling H = (0/2) (a 01 + 
cr 10 ), and the density matrix equation has an additional term: 

dtp ^d t p-^ {a u p - 2a el pa le + pa 11 ) . (11.19) 

The pumping term for p ee , which previously was Tpoo now becomes: 

d tPee\ pump = Ippll (11.20) 

If we assume 7 P is very rapid, then in the same spirit as the above derivation 
we have: 

in 

dtpu = — — (poi - Pw) ~ IpPn, (11.21) 
dtPio = — 2" (/ooo - Pu) - y/oio, (11.22) 
and then setting these time derivatives to zero one has: 

in . . n 2 . . . 

Pio = — (Pu - Poo) , = (pu - poo) - TpPli (11.23) 

1p 1p 

we can then find pn,pio in terms of poOi which allows us to derive an 
effective incoherent pump rate of: T e ff = 7 P 2 / (O 2 + 7p). If O is small, then 
r c g oc n? /7 p , but if pumping is too large, there is an additional constraint 
on the pump rate, set by the relaxation to e. If 7 P is not large compared 
to other timescales (i.e. if g were to become very large), then it would no 
longer be possible to eliminate the dynamics of this level. 

Maxwell Bloch equations 



Having shown that Eq. (11.2) corresponds to the same physics as discussed 



in the previous lecture, let us now discuss a different approximate descrip- 



tion of Eq. (11.2). This is the semiclassical approximation, which should 
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be valid when numbers of photons and atoms are large. We will thus use 
the density matrix equations to write equations for expectations of three 
operators: 



a 



(a), V = N A (an, M = N A (a e 



a 



(11.24) 



corresponding to the electric field, polarisation and inversion respectively. 
The approximation scheme will be that we replace expectations of prod- 
ucts of operators by products of expectations. From the density matrix 
equations of motion, one immediately finds: 



d t a 



.g k 

-i-V a, 

2 2 ' 



dtP = +i^aM - jtP. 



(11.25) 



For the equation for M we must write: 

■ 9 



i| {aV* - a*V) + r(a 00 ) - 7 (a ee ) 



d t (<J ee ) 
dt{a B0 ) 

d t (a 00 ) = -r(a uu ) + 7 (0 +i{ct") 



+i^(aV* -a*V) - -/(a") 

,00 \ 



We will make an additional assumption here, that T <C 7, in which case 
{a 00 ) ~ 1, and we can identity r = TN A - Thus, in this case, the difference 
of the above equations becomes: 



dtM 



-ig (aV* - a*V) + r - 7JV 



(11.26) 



We may now consider Eq. (11.26) and Eq. (11.24), and look at their 



steady states, and normal mode spectrum, to understand their behaviour. 



Steady states of Maxwell Bloch equations 



For a steady state, Eq. (11.25) immediately implies that: 



V 



K 

i—a 
9 



2jt 



Ma. 



(11.27) 



This implies either that a = or that M = 2^t^/g 2 - The latter solution 
(corresponding to lasing) implies that inversion locks at its threshold value; 
this is known as gain clamping. For the last equation, we have either that: 

Below threshold, a = 0. In this case, dtM = implies N = r/7 (hence 
threshold occurs at r = 2/-t(77 t /g 2 ), as found previously). 

Above threshold, J\f = A4at = ^tn/g 2 . In this case, we may use the gain 



equation as a condition for a, using Eq. (11.27) to eliminate V, and 
then writing: 







-2n\a\ 2 + r - 7A4 



sat. 



a 



2k 



sat 



[11.28) 



Thus, above threshold \a\ 2 increases linearly with pump rate in this mean 



field theory. This is illustrated in Fig. 11.3 
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yN. 



Figure 11.3: Left: Field strength, and inversion arising from 
mean field theory. Right: Decay rates of fluctuations of Maxwell- 
Bloch equations. 



Fluctuations of Maxwell-Bloch equations 

The instability of the non-lasing state at threshold can be understood 
by considering fluctuations, and linearising the Maxwell-Bloch equations 
around a = V = 0,N = r/~/. One then finds: 






d t \ bV = igr/2 7 - 7t \ \ SV \ (11.29) 



Thus, inversion fluctuations are clearly damped. The photon field and 
polarisation have collective oscillations at a frequency to given by: 




- = -Mt + ^)±m/(t-|) 2 + C7- (1O0) 

This frequency is purely imaginary, so the modes are either growing or 
decaying 1 . If r > r t h, one of these roots becomes positive, and such fluctu- 



ations will grow. These modes are shown in Fig. 11.3 

If 74 3> ft, the two modes are always well separated, and the slow dy- 
namics can be found by eliminating the matter degrees of freedom, writing: 

<p = ^-m r= ( 7+ ^!| a |2W (n.31) 
^it v it j 

If we then define n = \a\ 2 , the remaining equation becomes: 

d t n = - — kti (11.32) 

2 77t + 9 n 

This obviously has the same steady states as found above, the fluctuations 
are not identical, but become so asymptotically if 74 ^> k. 

1 NB, because of the pumping and decay scheme chosen, inversion is always positive 
in this system, Af = r/7 — a more realistic model, with thermal transitions between e 
and g might allow for negative inversion. If inversion is large enough and negative, then 
Rabi oscillations are possible, as in the cavity QED discussion earlier. However, the Rabi 
splitting collapses when the system is inverted, so these are not relevant in the current 
model of pumping. 
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11.2 Spontaneous emission, noise, and j3 
parameter 

In the previous section we found a semiclassical equation for n, the mean- 
field photon number. This equation is valid as long as numbers of photons 
are large. When numbers of photons are small, spontaneous emission be- 
comes important, and effectively n — > n + 1 in the previous equation. This 
equation, with n — > n + 1 can also be seen as the result of writing: 

d t (n) = (B n - C n ) (11.33) 

with B n , C n as derived in the previous lecture, and re-derived above. If we 
make the mean-field assumption (B n ) ~ B/ n \, and we write n for (n), then 
we have: 

r g 2 (n+l) 

2 llt + 5 2 (n + l) V ; 

We will investigate here the effects of this +1 — we should be aware however 
that whenever the effects of this term become significant, so do fluctuations 
(as we will investigate more fully below), in which case no mean-field theory 
can be applicable — i.e. it is necessary to study the full probability dis- 
tribution, not just (n). However, if we can understand when the +1 term 
will matter, then we can see the conditions under which semiclassics is a 
reasonable approximation. 

The important question regarding the effects of the +1 is to understand 
whether single photons will affect the dynamics. This in turn means the 
question of whether one photon has a significant effect on saturating the 
transition, i.e. what is the ratio of g 2 /"f"ff It is standard to define the 
parameter j3 slightly differently, as: 

P = 9 \ 2 - (H-35) 

lit + g 2 

If g 2 <C then single photons have a negligible effect, and this correspond 
to p <C 1. If g 2 3> 77t, single photons are important, and — > 1. In terms 



of P, we may rewrite Eq. (11.34) as 



rP n + 1 

dm = nn 11.36) 

2 1+pn v ' 

Based on the semiclassical threshold discussed above (valid if p <S 1, and 
so P ~ 9 1 Hit) w e will define rth = 2k/ p. In terms of this threshold, we 
can then solve the steady state number of photons to be given by r{n + 1) = 
r t h«(l + Pn), with the solution: 



-- 1± \/(-- 1 ) 2+4 ' 3 - 



(11.37) 



This equation (normalised by plotting Pn) is shown in Fig. 11.4. For small 
P it is clear that one recovers the sharp step like feature at r = rth- However 
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as (3 increases, the step becomes more smeared, and eventually, at ft = 1, 
one has n/3 = r/r^ for all r, i.e. no step survives at all. This indicates 
that in the limit of strong coupling, there is no clear "laser threshold" , but 
rather a thresholdless lasing behaviour. 
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Figure 11.4: Numbers of photons vs pump rate (normalised by 
7"th, and Fano factor, for different values of /3 




Fano factor 

As well as the sharpness of (n) vs pump power, further insight into the role 
of (5 can be found by considering the full probability distribution, and the 
relative size of fluctuations. We choose to measure the Fano factor, defined 
by: 

f =<"•>-<">'. (11.38) 
(n) 

This parameter measures deviations from Poissonian statistics, as F = 1 
for a Poisson distribution. The fluctuations can also equivalently be defined 
by 52(0), where: 

{o)a^aa) (n 2 ) — (n) 
(ata) 2 (n) 1 

Another parameter occasionally quoted is the Mandel Q parameter. The 
three parameters are related by Q = F — 1 = (n)(<?2(0) — 1). We will 
consider only F below. 

Let us now determine F from the full probability distribution as derived 
in the last lecture, and assuming [3 is not too large. Far above threshold, 
the distribution is Poissonian, and so F — > 1. Below the threshold, we may 
use the Bose-Einstein distribution p n oc (r/Vth) n to give: 

(»> = (» 2 > = r/ ;; h(1 ! r/ S h) (n.40) 

l-r/r th (l-r/r th ) 2 

and so: 

r/r th 



52(0) = ,„ + ^ 2 = JZX2 • (11.39) 



F = l+ , ' (11-41) 
1 - r/r th 
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which again becomes 1 is far from threshold. 

Exactly at threshold, one can make an approximation for small (3 that: 

Pn oc JJ(1 + (3m)~ l ~ exp I -f3^m ) ~ exp (-(3n 2 /2) . (11.42) 

m \ m / 

This is valid if /3 <C 1, in which case, p n has decayed to a small value well 
before (5n ~ 1. With this one-sided Gaussian distribution, one then finds: 

<»>-\/5 F -$a{f*-£)- (11 ' 43) 

Thus, if /3 -C 1, there is a sharp spike of Fano factor at the transition — the 
uncertainty in photon number becomes very large here in the semiclassical 
limit. On the other hand, everywhere else, the Fano factor is small. The 



full behaviour of the Fano factor vs (5 is shown in Fig. 11.4 It is clear that 
for (5 ~ 1, there is no clear spike, and instead one has a broad region near 
the ill-defined threshold where fluctuations are large. 



Micromaser and noise 




The above results for the laser have F > 1 in all cases, meaning that the 
noise is always at least Poissonian, and is sometimes worse. On the other 
hand, the micromaser that we discussed previously can have sub-Poissonian 
noise. Recalling that in that case, the probability distribution is given by: 

n/r\ /sin 2 flT-v/m/2\ 
(«)( m )• (1L44) 

There are two important features of this distribution; the first is the irregu- 
lar dependence on gr, which was commented on before. This is illustrated 
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in Fig. 11.5 More importantly, there is the dependence on n, which is what 
can lead to number squeezing (i.e. sub-Poissonian noise). 

The number dependence of the gain means that the distribution is at- 
tracted to particular number states. This is illustrated schematically in 
Fig. 



11.6 



which plots rsm 2 {gT\Jn + 1/2) — kit, as a function of n. Because 
the width of the attraction regions varies in a non-smooth way, the width 
of the distribution is not simply connected to the mean number of photons. 
Thus, the variance can be far less than it would be for a Poisson distribution 
with the same mean number of photons. 




Figure 11.6: Gain vs n for micromaser equation 



The number dependence of the gain is also responsible for the discon- 
tinuous jumps that can occur on increasing r. As r increases, new at- 
tractors appear, i.e. new peaks of sin 2 gr\Jn + 1/2 fall within the range 
n > r/n. The photon number and Fano factor vs pump rate is shown in 



Fig. 11.7 This kind of mean-field-like argument can describe the possible 
steady states reasonably well (under the standard n>l validity condition 
of mean field theory). It cannot however describe the points of transition 
between the different possible steady states. This can be clearly seen in 



Fig. 11.7 



11.3 Single atom lasers 

Armed with the above ideas of the Fano factor, describing noise, and the 
validity of semiclassics, and the (3 factor describing how saturation effects 
occur, we can now consider cavity QED lasing systems, where small num- 
bers of atoms are involved. 



Firstly, we should note that the density matrix equation in Sec. 11.1 
could reasonably describe a single atom. However, it does rely on the 
separation of timescales, and if very strong coupling is required, then the 
approximations used in deriving it will become invalid. To see why this 
might occur, let us rewrite the threshold condition as: 

N A Tj 2k 

r =rT7 = 7 (1L45) 



If N A = 1, and 7 cannot be significantly tuned (as it describes intrinsic 
decay rates of atoms), then to reach threshold with a single atom, it is 
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Figure 11.7: Numbers of photons and Fano factor vs coupling gr 
for a micromaser. Blue squares show the most probable number 
of photons. Dotted lines correspond to the mean-field results, 
r sin 2 (gr^n + 1) = nn. (gr is fixed at 2.5) 

generally necessary to increase /3, so that one will no longer have (3 <C 1. In 
this case, fluctuations will matter more, the behaviour will be thresholdless, 
and if g becomes very strong then no separation of timescales is possible. 
Thus, while single atom lasing need not automatically imply a breakdown 
of the above description, reaching threshold with a small number of atoms 
generally requires strong coupling, and this in term means large /3, hence 
thresholdless lasing. 

Optical cavity, real atom 




Trap 

Figure 11.8: Schematic of single atom optical lasing scheme, and 
level scheme 



Lasing from a single atom, using an optical transition and a regu- 
lar cavity has been achieved by [35 . The system and level scheme are 
shown schematically in Fig. |11.8| The parameters for this experiment were 
(g,K,j) = 2tt x (16,4.2,2.6)MHz. This implies ~ 0.97, and so the lasing 
will be thresholdless. In addition, because of the strong coupling, the adia- 
batic elimination of the non-lasing levels is not a good approximation in this 
system. This allows the possibility of sub-Poissonian behaviour, which was 
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observed experimentally, and can be described theoretically by modelling 
the full four level system illustrated in Fig. 11. 8| 



Josephson qubit micromaser 

The other example we will discuss is a rather different system and rather 
different method of lasing. This system is a superconducting qubit as the 
artificial atom, coupled to a mechanical resonator, thus this is phonon lasing 
rather than photon lasing|36] [37]. The circuit is shown schematically in 
Fig. 



11.9 



Resonator 



SC Island 



Vg+ v ds 



IE, 



2 



Figure 11.9: Circuit diagram of single artificial atom laser 

The two-level "atom" here consists of two charge states of the island, 
having either zero or two extra electrons on the island (two electrons be- 
cause of Cooper pairing, and the consequent odd/even energy gap). Thus 
one has cr + = a 20 etc. This gives the Hamiltonian: 



H = [E c + A(a + a ] )]a z + Eja x + ua)a 



(11.46) 



where Ec depends on the capacitance of the island, and A describes the 
shift due to the location of the mechanical resonator. An equivalent model 
would also describe a capacitively coupled stripline resonator. As well as 
this Hamiltonian behaviour, the bias Vd s leads to a single-electron current 
through the island, with individual electrons hopping in a two step pro- 
cess. Eliminating the singly charged state (which decays fast), one has the 
effective density matrix equation: 



dtp = —i[H, p] — — ( a)ap — 2apa) + pa) a 



(a + a p — 2a pa + + pa + a 



(11.47) 



Lasing can occur if Ec + uj — 0, in which case, there is a cycle of 
incoherent transition from 2 — ► 0, followed by coherent transition from 
— > 2 with emission of a phonon to conserve energy. However, because the 
photon number affects the energy of the state Ec, this leads to micromaser 
like behaviour. In the theoretical work |37|I38]. the parameters are however 
in a regime where mean field theory can describe some of the features seen 
by full density matrix evolution. 
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11.4 Further reading 

Discussion of density matrix equations and Maxwell-Bloch equations for 
our model laser can be found in many textbooks, e.g. Meystre and Sargent 
III [12], Yamamoto and Imamoglu |llj . The effect of the (3 parameter, and 
the nature of the threshold is discussed in Rice and Carmichael p9j. A 
review on cavity QED particularly focusing on micromaser behaviour can 
be found in Walther et al. [40] . 



Lecture 12 



Three levels, 

electromagnetically induced 
transparency, and coherent 
control 

Most of the lectures so far has been restricted to two-level atoms. While 
there many are obvious extensions to the case of more than two levels, this 
section focuses only on one phenomena which requires three levels: That is 
the coherent interference of pathways for atomic transitions, and the possi- 
bility of suppressed transitions. By using an applied electric field, one can 
suppress transitions between levels, and thus prevent absorption, leading 
to Electromagnetically Induced Transparency. This lecture starts by dis- 
cussing the phenomena in terms of the density matrix evolution, allowing 
incoherent decay of atomic energy levels, and considering radiation semi- 
classically. In this way, we will see that conditions under which the group 
velocity of light is significantly reduced while absorption of light remains 
small is possible in this system — a combination of features that does not 
occur for a two level system. Then, by considering the atomic wavefunction 
evolution, it is possible to understand the underlying interference phenom- 
ena which lead to the cancellation that prevents absorption. Finally, we 
return to a quantum optics description of the states that are free from 
absorption — dark state polaritons — discussing how the nature of these 
coherent quantum states can be controlled. 



12.1 Semiclassical introduction 



We consider the level scheme shown in Fig. 12.1 Dipole matrix elements 
exist between the ground state b and the excited state a, and also between 
the metastable state c and excited state a. The transition between b and 
c is however dipole forbidden. The states a and c are coupled by a strong 
"pumping" radiation field, described only by the effective Rabi splitting it 
induces, Q p . The transition between a and b is then probed by a weaker 
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field. We allow the possibility of incoherent decay from a to b, and also a 
(much weaker) decay from c to b. The coherent part of the evolution can 




Figure 12.1: Scheme of atomic energy levels and allowed, and 
driven, transitions between them. 



be described by the Hamiltonian: 



H = | - E d* ab e iuJt oob | . (12.1) 



Susceptibility; slow light. 

Let us consider the absorption and emission coefficients, by finding the 
susceptibility of the three level system; i.e. finding how the polarisation, 
P = (N SLtonis /V)d a bPabi depends on the field E applied resonantly with the 
a <-> b transition. We assume the field E is weak, so we consider only linear 
response, and can assume the populations of the a and c levels remain small. 
We are interested in the evolution of the density matrix, p = —i[H, p] +Cp, 
where C represented a superoperator, describing the incoherent evolution, 
described in lecture [4] Since Trp = 1, small populations of a, c imply 
Pbb ~ 1 and p ac ~ 0, and Hermiticity means the only relevant components 
are p a b and p c b- Thus: 

Pab = — i [{H aa • Pab + H a b • 1 + H ac ■ p c b) 

- (0 • H ab + Pab ■ Hbb-0- Hcb)] - ^abPab 

= - (iujab + j a b) pab + Xlpe**** p A + iEd ah (T^\ (12.2) 
Pcb = -i [{H ca ■ Pab + • 1 + H cc ■ p cb ) 

- (0 • H ab + Pcb ■ Hbb - • Hcb)] - IcbPcb 

= - {ioj cb + jcb) Pcb + i^y^pab. (12.3) 

Let us now assume that the pump is resonant, so that uj p = oo ac = w a — u c . 
Then, rewriting p a b = Pab&~ luli and p c b = Pd>^ ~ U3p ' t i and defining the 
probe detuning A = uj a b — w, we have: 




Pab \ _ ( -(iA + Jab) i% \( Pab \ ( iEd ah \ 

Pcb ) V in * P -{i^ + lcb) J \ Pcb J \ )■ 



(12.4) 
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This has the form X = MX + A, and we want the steady state solution (af- 
ter decay of any initial transients), which is clearly given by X = —M~ 1 A. 
Thus we can extract the polarisation by finding p ab : 



P = 



N atoms iVatoms, (lA + J cb )lEd ab _ iwt , 10 _s 
Tj. dabPab = T~ r aabTTj— w- A , n , i n i 2 e ■ l 12 - 5 J 

Then, noting the time dependence of E{t) = Ee~ tujt , one can easily identify 
the susceptibility, x by P = X £ oE. i.e. 

= ^atoms \d ab \ 2 j(zA + ~f cb ) 

X ~ V e {iA + lab ){iA + lcb ) + \£l p \z- 1 j 

Let us now discuss the properties of this susceptibility. It is clear that 
x(— A) = — x(A)*, thus if x = x' + ^x" i it is clear that x' is an °dd 
function of A (and vanishes at A = 0), and x" an even function. However, 
x(A = 0) oc 7 C 6, which we assumed earlier to be small, as this transition is 
dipole forbidden. To be precise, x"(A = 0) oc r y c b/(lablcb + l^p| 2 ); thus the 
absorption part is small if *y cb <C £l p ; the coherent driving field has induced 
transparency. The full real and imaginary parts of susceptibility are plotted 



in Fig. 12.2 both with and without a driving field S7 f 




Figure 12.2: Real and imaginary parts of susceptibility, plotted 
for 71 = 0.5,73 = 0.01, and two values of fl p as indicated. 



Since both real and imaginary parts of x are zero, or very small at 
A = 0, let us consider linear expansion; clearly only the real part exists to 
linear order. Setting ^ cb = 0, one trivially finds: 

X = - iVatoms|4b|2 A + 0(A 2 ). (12.7) 

V e \n p \^ 



The linear part of susceptibility relates to the dispersion of the probe field; 
n(A) = ^/l + x(A)) which can be used with u = uj ab — A, and n{uj)uj = ck. 
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One thus has phase velocity v p = c/re(A) and group velocity v g = c/[n 
(uj a b — A)dn/dA]. Assuming A remains small, one can expand: 



n 



1 



N atoms 

V 



to 



Thus, to leading order: 



A 



dn ^ -/Vatoms l^abl 2 1_ 

dA ~ V e 2\n p 



(12. 



c, 



i , N u ab \d ab \ 2 ' 
1 T |n p p 2e V 



(12.9) 



where we have rearranged factors to identify the atom-radiation coupling 
strength g 2 = u a ;,|(f a b| 2 /2£o^ [cf. the definition in Eq. (1.39), in the case 
e = u ab \. Thus, if Ng 2 j 



ojk = e = 0J a b\. Tnus, it l\g" /\\l p \ is large, one can have small absorption, 
yet large enough dispersion to significantly reduce the group velocity, and 
thus have non- absorbed slow light. 



Decay length of slow light 

The previous analysis shows that near zero detuning, one can have a small 
group velocity, but small absorption. Let us now be more specific about 
how small the absorption is. What matters in practice is the attenuation 
of the signal as it propag ates, following 1(1, oj) = /(0)e -fc,x "^. Since the 
attenuation is frequency dependent, this may also lead to distortion of the 
propagating signal. The frequency dependence of the transmission intensity 



is shown in Fig. 12.3 




-5 -4 -3 -2 -1 1 2 3 4 5 
Detuning from resonance 

Figure 12.3: Transmission after L = (50/7r)A, for the same pa- 
rameters as in Fig. |12.2| 



Considering the case 7 C & = 0, for small A, the second order expansion 
yields: 



X « 2i — —A 



N 9 2 A 2 lab 

ck I i^rt I ^ 



(12.10) 
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where we have written uj a b = ck. One can thus identify a frequency depen- 
dent decay length: 

1 1 l^l 4 

L = c——,—, 1 -^-. (12.11) 

12.2 Coherent evolution alone; why does EIT 
occur 

Let us now ignore the decay rates, j a b and 7ft c , and look at the evolution 
of atomic states alone. Thus, since we are not interested in the dielectric 
response, we may denote Ed a b = £l s as the signal field. Considering the 
fully resonant case, to = uj a b and u> p = Lu ac , one can write: 

c a e~ iuJat \a) + c b e-™ bt \b) + c c e~ iuJct |c) (12.12) 

(12.13) 









id t 


Cb 


H 










l \b) + c c e 


-iu) c t 




C a 


T ) ( 


Cb 




Cc 



One can easily find the eigenvalues and vectors of this problem; in terms of 
VLq = y^^sl 2 + \Q p \ 2 one has: 





A-0: —I n f , ]: X = ±n: — | <>.; j. (12.11} 

Thus, if one starts in the state (£l p \b) — £l s |c)) /0,q, then the evolution never 
reaches the excited state \a). Since the only substantial decay constant was 
from the state |a), this then means one has no absorption: The evolution 
is trapped in this dark state, and shows no decay. One may also describe 
this process as destructive interference between the transitions \b) — > \a) 
and \b) — > \a) — > |c) — ► \a) etc. 

Since such a dark state exists, from which no further excitation is pos- 
sible, the full evolution (with decay rates) can be understood as follows: 
Regardless of the initial state, decay out of states \a) allows one to reach 
various superpositions B \ b) + C |c). Eventually, one will reach dark state, 



after which no further evolution can occur. [See Q. 12.1 . 



12.3 Coherent control of quantised radiation: 
Dark state polaritons 

The previous sections show that the existence of the classical driving field, 
resonant with the c <-> a transition modify the response seen by a probe 
field at the b «-> a transition: One can have vanishing absorption, while 
still having strong dispersion, allowing propagation of "slow light". In this 
section, we consider the probe field quantum mechanically, and ask how the 
classical driving changes the evolution of quantised radiation. As illustrated 



in Fig. 12.4, we consider the driving as still classical, but the probe field is 



quantum mechanical. 
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Quantum 
Field 




Classical 
Driving 



Figure 12.4: Illustration of which transitions are classically 
driven, and which described quantum mechanically 



In order to consider propagation through the EIT medium, we consider 
a continuum of radiation modes, as well as explicitly including a sum over 
atomic states. Note that in the earlier treatment, the sum over atoms 
occurred implicitly via the factor iV atoms in Eq. (12.6) Assuming resonance 
for the driving field, and writing a^ y for the three by three matrix describing 
transitions from level y to level x of atom i, one can write: 



H = Yl cka W + Y 



k>0 



n* p af a + H.c. 



(12.15) 

Then, making a gauge transform to a rotating frame one can remove the 
explicit e~ luJabt time dependence, in favour of replacing ck — > ck — uj a b. 
Then, one may make a Fourier transform of the photon field, defining: 



a{x) = y^ctfce 



i(k-u) ab /c)x 



(12.16) 



By writing k — ui a b/ c rather than k, one recovers the correct energy from 
the derivative w.r.t. x. Thus, one can write 

H = J dxa\x)(-icd x )a(x) + \-go\ a a ] { Xl )e~ iwM/c - n* p af a + H.c. 

(12.17) 

Note the factor of e~ lulabXi ^ c appearing in the photon-atom coupling term. 
Let us introduce collective variables: 

X{x) = e?5{x - Xi) 



P(x) = £ o?8{ 



-iuj ab Xi/c 



S( x ) = y £ j a\ c 5{x-x i )e- iuJ - 



b%i/c 



;i2.i8) 



The definitions here are not the most symmetric, but are physically moti- 
vated: P represents a polarisation wave, while S represents a "spin wave" 
- this language is appropriate to cases where the low lying levels b, c arise 
from opposite spin states. The definitions of P(x) and X(x) are easily 
understood, as they allow one to rewrite Eq. ( 12.17[ ) as: 



H= / dxa) \x){-icd x )a{x) - gP{x)a\x) + %X(x) + H.c. . (12.19) 
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To understand the definition of S(x), we must consider commutation 
relations. Writing commutators using [ct*- 5 , a kl ] = 5jk<7 il — buo k i , the cross 
commutators are given by: 



\P{x),X\x')\ = S(x)S(x-x'), 
[S{x),X{x')} = P{x)5{x-x'), 
[P\x),S(x')\ = X\x)5{x-x'). 



(12.20) 



Other than these, and their conjugate forms, all other cross commutators 
vanish. Thus, it is clear that the phase factors in the definition of S{x) in 



Eq. (12.18) is required to give these simple forms. 



We now have a tractable form of the quantum Hamiltonian, which allows 
us to consider how an initial state will evolve. Let us solve the simpler 
problem, of finding the quantum states that are trapped in the absorption- 
free subspace. It is clear that the atomic ground state, with no photons is 
an eigenstate, which we will denote |0). Let us consider excited states of 
the form: 



$) = $t | ) = f dxcf>(x) (aa)(x) + fiS^xj) |0) , 



(12.21) 



and ask for the condition under which no excitations to the excited atomic 



state, a, occur. Note the distinction to the discussion in Sec. 12.2; in that 



case we wanted superpositions of b and c states so that in the presence of 
a given probe electric field there were no excitations; in the current case 
we want a superposition of radiation and spin- wave excitations. Since the 



trial state in Eq. (12.21) contains no excitations, we require only that the 



time derivative of number of excitations vanishes. Thus, the condition we 
require can clearly be written as: (0| P \H, $t ] |0) = 0. Let us consider the 



commutator, using the results in Eq. (12.20), along with [a(x), cv(x')} 
S(x — x'): 



H, (aatfx) + (3S\x) 

a (icd x a)(x) - gP j (x)^j + (3 (ga\x)X(x) - £l p P ] (x)} . (12.22) 



The term proportional to w(x)X{x) is second order, involving both atomic 
and radiation excitations, for low intensities it can be neglected. The num- 
ber of excitations then vanishes under the condition: ag + (3£l v = 0. If we 
wish to find a normalised wavefunction, we require: 

S(x - x') = [$(x), &(x')\ = \a\ 2 5{x - x') + \(3\ 2 ^ S(x - x'). (12.23) 



Thus, one has 1 = \a\ + in addition to the previous condition, thus: 



a 



(3 



y^WTg^N' 



(12.24) 
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From Eq. (12.22), we find this state evolves according to the Heisenberg 



equation of motion, 

irQ 

- id t &(x) = [H, *t] = p 2A A aHx). (12.25) 

To eliminate cv(x), we wish to write this in terms of the field <J?t(x). The 
field & is the linear combination of a) and which is dark — i.e. never 
undergoes a transition to the excited atomic state; another orthogonal lin- 
ear combination of at and exists, let us call that W. Then, one can 
write: 

* f ) = ( 3% ) ( JfSl ) ■ (12 - 26) 



Here we have included a factor of 1/yN required to achieve bosonic com- 
mutation relations into the definition of a bosonic spin operator Sft /y/N. 
Inverting this to find a) in terms of <J?t and , and assuming an initially 



dark state so that ^ — > 0, Eq. (12.25) becomes: 



id t ^(x) = - iC ^ p d x - S &{xY (12.27) 



i.e. 



Thus, we have found a darA; staie polariton [5T], a state which is never 
excited to the radiative atomic state o, and is half light, half spin-wave 
excitation. It propagates at a reduced velocity, which is the group velocity 



found in Eq. (12.9). Importantly, one can control the nature of the state 
by the ratio of f2 p to gyN. In a strong driving field O p 3> gV~N, the 
excitation is almost pure photon; for a vanishing field O p -C gV^N, the 
excitation becomes pure spin wave. Thus, starting with a large Q p , one can 
inject an arbitrary quantum state of photons. By adiabatically reducing f2 p , 
the quantum state is exactly transfered to the spin waves. These propagate 
slowly, until O p is increased again, and the state transfered back to photons. 
Hence, coherent control allows the trapping of quantum states of light. 



12.4 Further reading 

The general topics of EIT, dark states, and coherent control in three level 
systems are discussed in a review by Fleischhauer et al. [42 . The question of 
coherently transferring the quantum state of light to "dark state polaritons" 
is introduced in Ref.[41j, and discussed further in Ref . |43] . An interesting 
application of this idea is described in Ref. |44j . where a quantum state of 
light is transported by transferring atoms between two condensates, and 
then releasing the state of light from the second condensate. 
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Questions 



Question 12.1: Evolution to dark state 



As mentioned in Sec. 12.2 decay should be responsible for projecting 



into the dark state. Considering Eq. (12.4); find the evolution to the dark 



state if one starts in the ground state, i.e. find the general time dependent 
solution starting from p a i = p c b = 0. Sketch the time dependence of 3?[/? C 6] • 
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